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1  INTRODUCTION 

This  report  describes  development  of  a  software  package  that  simulates  tissue  damage  from 
small  arms  projectiles  and  fragments  penetrating  the  lower  extremities.  The  software  can  be 
interfaced  with  software  that  performs  a  detailed  biomechanical  assessment  of  residual  physical 
capability,  specifically  moment  generating  capability  of  joints  affected.  Finally,  a  graphical 
interface  is  provided  for  medical  simulation  and  training. 

Two  subcontractors  assisted  Mission  Research  Corporation  (MRC),  MusculoGraphics,  Inc. 
(MGI)  of  Evanston,  Illinois  and  GB  Labs  (GBL)  co-located  with  MRC  in  Fountain  Valley, 
California.  MGI  developed  the  interface  allowing  visualization  of  the  wound  tract.  This 
visualization  is  incorporated  in  the  virtual  anatomy  developed  by  MGI  as  part  of  the 
DARP A/MGI  Surgical  Simulation  for  Limb  Trauma  Management  Technology  Reinvestment 
Program  (TRP).  GB  Labs  performed  ballistic  testing  on  biosimulant  materials  to  develop  a 
phenomenological  database  as  well  as  help  validate  MRC  modeling.  This  testing  involved 
innovative  biosimulant  target  configurations  with  in  situ  instrumentation. 

Major  accomplishments  of  MRC  included  the  development  of  models  describing  soft  tissue 
damage  from  a  penetrating  wound.  This  description  includes  the  permanent  wound  cavity,  tissue 
transient  response,  and  hemorrhagic  zone  of  extravasation.  MRC  also  developed  models  that 
describe  the  interaction  of  bone  tissue  with  the  missile  including  mode  of  fracture  produced  and 
the  resulting  kinematics  of  the  projectile  in  three  spatial  dimensions  and  time. 

The  models,  technical  insights,  and/or  biosimulant  configurations  are  already  being  employed  in 
five  other  Government  projects.  In  the  DARPA/DLA/NRaD  Sensate  Liner  Development  for 
Combat  Casualty  Care  (NRaD/MRC  N66001-96-C-8641),  MRC  developed  non-invasive 
sensors  and  signal  processing  algorithms  that  can  acquire  the  acoustic  signature  of  projectiles 
penetrating  soft  tissue  and  determine,  in  near  real-time,  the  trajectory  of  the  projectile  through 
the  human  anatomy.  This  project  exploited  the  .soft  tissue  modeling,  biosimulant  technology,  and 
wound  profile  analysis  developed  in  the  DARPA/USAMRMC  project. 

In  a  project  for  the  Army’s  Soldier  System  Command  (USASSCOM/MRC  DAAH02-98-C- 
4012)  Development  of  a  Character  Simulator  for  Battlefield  Virtual  Environments,  MRC  is 
developing  a  virtual  human  that  incorporates  the  DARPA/USAMRMC/MRC  wound  ballistic 
penetration  algorithms.  As  a  subcontractor  to  Scientific  Applications  and  Research  Associates 
(SARA),  MRC  is  developing  models  and  performing  experiments  using  biosimulant  human 
phantoms  developed  in  the  DARPA/USAMRM/MRC  project  to  describe  the  biological  effects  of 
non-lethal,  high-energy,  toroidal  vortices.  As  a  subcontractor  to  Simulation  Technologies, 
Incorporated  (STI),  MRC  incorporated  for  the  USASSCOM,  new  wound  ballistic  and  casualty 
assessment  algorithms  in  the  Army’s  CASRED  code.  The  algorithms  were  based  on  the 
penetration  algorithms  developed  in  the  DARPA/USAMRMC/MRC  project. 

Finally,  as  previously  mentioned  the  DARPA/USAMRM/MRC  software  is  being  directly 
integrated  into  the  Lower  Extremity  Battlefield  Trauma  Simulator  being  developed  by  MGI  for 
DARP  A.  In  this  DARP  A/MGI  project,  the  MRC  software  is  being  used  to  provide  an  analytical 
description  of  tissue  damage  from  a  penetrating  wound.  This  is  necessary  because  although  an 
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entrance  and  exit  wound  may  be  empirically  observed  and  digitized,  intervening  tissue  damage 
may  only  be  surmised.  The  DARPA/USAMRM/MRC  project  thus  filled  this  void  by  providing 
physically  based  analytical  models  that  explicitly  describes  the  interaction  of  a  penetrating 
missile  with  tissue. 

There  are  two  major  shortcomings  of  the  DARPA/USAMRM/MRC  models  however.  First, 
while  the  distribution  and  sizes  of  ejected  bone  fi-agments  in  soft  tissue  is  modeled,  the  finhatirpH 
tearing  of  soft  tissue  resulting  from  the  ejected  bone  fi’agments  acting  as  stress  concentrations  in 
the  projectile  wake  is  not  modeled.  Second,  the  correlation  of  the  models  with  autopsy  data  was 
not  as  extensive  as  what  we  would  have  liked.  These  shortcomings  were  not  due  to  any 
fundamental  technical  obstacles  but  rather  due  to  a  lack  of  program  resources. 

This  report  is  divided  into  five  major  sections,  as  well  a.5  two  Appendices.  The  second  section. 
Section  2,  describes  development  of  the  analytical  models  and  the  resulting  BONEGEL  code. 
The  BONEGEL  code  describes  the  physical  interaction  between  projectiles,  soft  tissue,  and 
bone.  Sections  2.1  through  2.4  describe  development  of  the  analytical  methodologies  and 
insights  that  are  incorporated  into  the  BONGEL  software.  Section  2.5  describes  the  detailed 
structure  of  BONEGEL.  The  third  section.  Section  3,  describes  the  experimental  effort  that 
supported  analytical  model  development.  Section  4  describes  development  of  the  graphical 
interface  and  incorporation  of  the  BONEGEL  code  into  a  software  package  that  enables  reading 
tissue  geometries  and  properties  from  a  virtual  anatomy  based  on  the  visible  human.  Section  5 
details  conclusions  and  recommendations  from  the  entire  effort.  Appendix  A  includes  a  paper 
originally  presented  st  Medicine  Meets  Virtual  Reality  IV.  This  paper  gives  an  overview  of  the 
phenomenology,  modeling,  experimental  efforts,  and  accomplishments  of  this  program. 
References  to  other  presentations  made  on  behalf  of  this  program  are  listed  below.  Appendix  B 
lists  the  MRC,  MGI,  and  GBL  persormel  involved  in  this  effort. 

*  ft-  D.,  Eisler  and  A.  K.  Chatteijee,  Analytical  Simulation  of  Wound  Tracts  from  Missile  Penetration, 
Proceedings  of  Seventh  International  Symposium  of  Weapon  Traumatology  and 
Wound  Ballistics,  St.  Petersburg,  Russia,  20-23  September  1994  (Unclassified). 

*  Eisler  and  A.  K.  Chatteijee,  Simulation  and  Assessment  of  Musculoskeletal  Trauma  from 
Ballistic  Penetration,  MEDICINE  MEETS  VIRTUAL  REALITY,  m.  Sponsored  by  the  University  of 
California  Medical  School  at  San  Diego,  the  Commission  on  the  European  Union,  and  the  U.S. 
Advanced  Research  Projects  Agency,  San  Diego,  21  January  1995  (Unclassified). 

*  A.  K.  Chatteijee  and  R.  D.  Eisler,  Ballistic  Penetration  into  Gelatin,  IMPACT,  WAVES,  AND 
Fracture,  Proceedings  Of  The  Wemer  Goldsmith  Symposium,  sponsored  by  the  /^plied 
Mechanics  Division  of  the  American  Society  of  Mechanical  Engineers  and  the  University  of 
California  at  Los  Angeles,  Los  Angeles,  28  -  30  June  1995,  ASME  .^plied  Mechanics  Division, 
Volume  205,  pp.  9-20, 1995  (Unclassified). 


•  R.  D.,  Eisler,  A.  K.  Chatteijee,  and  G.  H.  Burghart,  Simulation  and  Modeling  of  Penetrating  Injuries 
from  Small  Arms,  published  in:  HEALTH  CARE  IN  THE  INFORMATION  AGE:  FUTURE  TOOLS  FOR 
Transforming  medicine,  presented  at  the  Medicine  Meets  Virtual  Reality  4  International 
Synqjosium  sponsored  by  The  University  of  California  School  of  Medicine,  the  Advanced  Research 
Projects  Agency,  American  Psychiatric  Association,  Institute  for  Telemedicine,  Society  for 
Minimally  Invasive  Suigery,  Society  of  Gastrointestinal  Endoscopic  Surgeons,  and  Society  of 
Cardiovascular  and  Interventional  Radiology,  San  Diego,  California,  17-20  January  1996 
(Unclassified). 
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2  Analytical  Development 

During  the  course  of  this  program  various  formulations  of  ordnance  gelatin  (250A  from  Knox 
and  Kind)  were  used  during  ballistic  experiments  as  surrogates  for  human  tissue.  These 
experiments  are  elaborated  in  Section  3.  Ordnance  Gelatin  was  selected  based  on  two  &ctors. 
First,  ordnance  gelatin  demonstrates  phenomenology  that  is  similar  to  penetration  in  human  soft 
tissue.  Second,  ordnance  gelatin  is  a  non-linear,  viscoelastic  material  that  has  many  of  the  same 
unknowns,  in  terms  of  material  properties,  that  in~vivo  human  tissue  has.  That  is,  a  fundamental 
technical  obstacle  that  thwarts  analysis  of  projectile-human  tissue  interaction  is  the  lack  of  in- 

vivo  dynamic  properties  of  living  tissue.  Further,  it  is  unlikely  that  these  properties  will  ever  be 
acquired. 

Another  complication  is  that  human  tissue  undergoes  phase  transitions  during  the  penetration 
process.  At  high  penetration  velocities  tissue  tends  to  act  more  like  a  viscous  fluid.  At  low 
penetration  velocities  tissue  tends  to  behave  as  a  viscoplastic  soUd.  Most  of  the  material  adjacent 
to  the  projectile  however,  particularly  at  intermediate  velocities  (which  is  where  most  of  the 
wound  tract  is  established)  is  a  multiphase  medium  with  even  less  defined  material  properties. 
Thus,  any  technique  that  relies  on  the  input  of  these  properties  for  human  tissue,  such  as  finite 
element  or  hydrocode  techmques,  at  least  for  the  analysis  of  “dynamic  events”  involving  human 
tissue,  is  of  questionable  utility.  Rather,  a  technique  is  required  that  can  analyze  the  dynamic 
interaction  between  a  penetrating  missile  and  human  tissue  in  terms  of  data  that  exists  and  is 
appropriate  for  this  type  of  interaction. 


An  assumption  in  this  program  has  been  that  if  ballistic  penetration  into  ordnance  gelatin  can  be 
analyzed  using  the  material  properties  available  or  measurable  for  ordnance  gelatin  targets,  then 
the  same  methodology  can  be  applied  to  human  tissue  since  the  same  material  property 
unknowns  and  phenomenology  applies  to  human  tissue.  Note  this  assumption  is  far  less 
restrictive  thM  what  the  lethal  effects  community  normally  assumes  in  using  ordnance  gelatin 
targets.  That  is,  the  lethal  effects  community  normally  assumes  that  penetration  in  human  tissue 
can  be  linearly  related  (i.e.,  scaled)  to  penetration  in  either  a  10%  or  20%  ordnance  gelatin 
target.  Without  addressing  the  veracity  of  this  controversial  issue,  we  make  no  such  assumption. 
We  merely  assume  that  the  phenomenology  and  that  the  material  property  unknowns  are  similar. 

Using  this  set  of  assumptions,  our  approach  has  been  to  parametrically  vary  the  gelatin 
formulation  to  document  the  change  in  gelatin  mechanical  properties  and  effect  on  projectile 
terminal  ballistics.  A  model  is  then  developed  that  describes  the  observed  trends  in  terms  of  the 
parameters  and  mechanical  properties  that  were  varied. 


This  model,  based  on  the  previous  assumptions,  should  have  the  same  parametric  form  as  the 
analogous  model  describing  the  same  phenomenology  in  human  tissue.  If  the  model  parameters 
are  formulated  in  terms  of  stable,  independently  measured  material  properties  that  are  available 
for  gelatin  as  well  as  human  tissue,  with  the  appropriate  substitution  of  parameters,  the  response 
in  human  tissue  can  be  described.  Validation  of  this  approach  is  achieved  by  comparing 
predicted  wound  tracts  to  autopsy  data.  It  should  be  noted  that  the  analytically  predicted  wound 
tracts  are  based  on  a  sequence  of  highly  non-linear  interacting  events.  These  events  include  the 
transient  expansion  and  recovery  of  tissue  in  the  projectile  wake.  This  expansion  and  recovery 
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results  in  a  hysteresis  effect  and  the  accumulation  of  inelastic  strain  which  results  in  the 
permanent  wound  tract.  Given  the  highly  non-linear  nature  of  the  projectile-tissue  and  projectile- 
gelatin  interaction,  any  agreement  of  analytically  predicted  wound  tract  characteristics  and 
observed  wound  tracts  would  greatly  exceed  random  probability  and  therefore  must  be  the  result 
of  correctly  modeling  the  physics  of  the  interaction. 

In  general,  analytical  predictions  successfully  correlated  with  the  experimental  results  on  gelatin 
targets  and  the  more  limited  autopsy  data  (which  we  would  have  liked  to  pursue  if  program 
resources  were  available). 

2.1  Resistance  Characteristics  of  Projectile  Motion  in  Gelatin  and 
Tissue  Simulants 

During  penetration  of  a  projectile  inside  a  target,  the  retardation  force  that  decelerates  the 
projectile  arises  due  to  its  contact  with  surrounding  target  material.  The  nature  of  these  contact 
forces  depends  on  the  kinematic  state  of  the  projectile  and  possibly  on  its  kinematic  history  as 
well.  For  example,  in  the  case  of  a  penetration  velocity  that  is  much  slower  than  the  speed  of 
sound,  Cp^  in  the  material,  the  retardation  force  is  approximately  equal  to  the  compressive 
fracture  tfeeshold  of  the  material,  and  hence  can  be  regarded  as  a  material  constant,  cf. 

Alternatively,  when  the  penetration  velocity  is  comparable  to  Cp ,  the  target  material  behaves  like 
a  viscous  fluid  medium,  and  hence  the  force  of  resistance,  F,  is  no  longer  a  material  constant  but 
depends  on  the  local  kinematics  and  possibly  on  kinematic  history  of  the  projectile  penetration. 
Most  experimental  studies  done  in  the  fluidized  state  of  the  material  assume  that  the  character  of 
F  depends  on  the  steady  state  fluid  drag  on  the  projectile  even  when  the  problem  is  temporally 
unsteady. 

Assuming  that  this  is  a  valid  approach,  F  is  proportional  to  in  this  regime  where  v  is  the  steady 
state  velocity  of  the  flow.  But  since  F  is  the  surface  integral  of  the  contact  forces  over  the  contact 
area,  it  is  physically  evident  that  for  unsteady  problems,  the  contact  area  may  not  only  depend  on 
the  instantaneous  velocity  of  the  projectile  but  may  also  depend  on  the  history  of  loading  prior  to 
the  current  state.  The  contact  area  also  depends  on  the  locations  of  boundary  layer  separations. 
Determination  of  boundary  layer  separations  in  unsteady  flow  is  a  very  complex  mathematical 
problem  that  has  not  been  solved  to  date  except  by  purely  numerical  algorithms  implemented  on 
super  computers. 

Here  we  have  taken  a  different  approach.  Our  goals  are  to  bypass  the  complex  step  of  numerical 
integration  of  nonlinear  differential  equations  and  to  develop  a  technique  based  on  a  semi- 
analytical  approach,  that  can  be  used  to  estimate  the  retardation  force,  F,  for  a  given  target  over  a 
given  range  of  entry  velocities. 

For  our  analysis,  we  assume  that  there  exists  an  experimental  database  that  contains  an  array  of 
penetration  depths  and  corresponding  entry  velocities  for  the  target  material  under  study.  Our 
analysis  also  deals  with  the  relationships  between  the  model  parameters  and  physical  parameters 
of  the  problem  as  much  as  practicable,  so  that  an  extension  of  our  analysis  to  unknown  materials 
(where  no  penetration  data  are  available)  can  be  made.  With  these  goals  in  mind,  we  start  with 
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the  assumption  that,  even  in  an  unsteady  problem,  F  depends  only  on  the  instantaneous  velocity 
of  the  projectile.  We  also  conceive  that  our  projectile  is  made  up  of  a  rigid  collection  of  particles 
so  that  F  can  be  determined  by  integrating  the  forces  on  such  particles  over  the  surface  of  the 
projectile.  Thus,  we  need  a  basis  of  determining  the  force  on  one  such  spherical  particle.  For  our 
purpose,  we  used  small  spherical  projectiles  to  simulate  such  particles.  Questions  may  arise  as  to 
whether  these  projectiles  are  suitable  for  simulating  particles  that  constitute  a  particular 
projectile.  The  success  or  failure  of  our  approach  will  determine  whether  the  use  of  other  particle 
simulants  is  warranted. 

Thus,  we  have  reduced  the  problem  of  determining  the  retardation  force  on  a  projectile 
penetrating  a  given  target  to  the  problem  of  determining  the  retardation  force,  F(v)  per  unit  mass, 
on  a  spherical  projectile  where  F  is  an  unknown  function  of  the  instantaneous  velocity  of  the 
projectile  inside  a  target.  The  penetration  problem  of  such  a  particle  inside  a  target  is 
schematically  shown  in  Figure  1. 

We  also  assume  that  we  have  an  experimental  database  relating  the  distribution  of  the 
penetration  depth  function,  S{vq)  (Figure  1)  as  a  function  of  entry  velocity,  v^.  Furthermore,  the 
low  and  high  velocity  behaviors  are  assumed  to  be  known  and  these  are  given  by 

lim  F{v)  =  Cf  ,  lim  F{v)  -  av^  ^ 

i;->0 


In  Equation  1,  both  ^and  a  depend  on  target  material  properties. 


Figure  1.  Schematic  of  projectile-target  interaction 


Let  us  first  establish  an  exact  relation  between  the  two-function  jfv)  and  F(v)  introduced  above. 
The  equation  of  motion  of  the  center  of  mass  (CM)  of  the  projectile  is  given  by 


v—  =  -F{v) 
dx 


(2) 
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Integrating  (2)  from  x=0  to  x=S,  we  have 


Since  (3)  is  valid  for  all  we  have,  after  differentiation  of  (3)  with  respect  to  v^, 

dv 


(3) 


(4) 


Equation  (4)  is  valid  for  all  values  of  the  local  velocity  v  of  the  projectile  inside  the  target,  and 
establishes  the  exact  relation  between  the  retardation  force  F(v)  per  unit  mass  and  the 
penetration  depth  6.  Equation  4  can  be  used  to  estimate  F(v)  from  the  experimental  database  of 
d(v)  and  V.  Below  we  present  a  method  of  estimating  F(v)  and  v  from  known  6(v)  and  v. 

2. 1. 1  Determination  of  the  Projectile  Retardation  Function,  F(v) 

2. 1.1.1  Low  Velocity  of  Penetration 

Form  the  first  equation  in  (1)  and  (4),  it  is  clear  that 

lim  5{v) »  (5) 

r->0 


and  hence  in  the  low  velocity  range  S(v)  can  be  expanded  in  a  Taylor's  expansion  as 

d(v)  =  v^(a  +  bv +  ciP'+...)  (6) 


We  can  now  estimate  the  constants  a,  h,  c  ...  in  (6)  by  first  dividing  the  experimental  data  by 
and  then  fitting  a  polynomial  to  the  reduced  data.  The  relation  between  the  constant  a  in  (6)  and 
the  physical  constant  <^can  be  found  by  using  (1),  (4)  and  (6).  These  give 


V 


Cf  =  lim  F{v)  =  lim 

v-^O  i;->0 

dv 


_1_ 

2a 


(7) 


Thus,  the  leading  term  in  the  expansion  of  S(v)  in  (6)  can  be  found  from  physical  data  of  ^that 
can  be  interpreted  as  the  ultimate  static  strength  of  the  material. 

2. 1 . 1 .2  High  Velocity  of  Penetration: 

From  the  second  equation  in  (1)  and  (4),  it  is  clear  that  S(v)  behaves  like  the  natural  logarithm  of 
V  when  the  velocity  v  is  comparable  to  the  sonic  velocity,  Cp^  of  the  target  medium.  Hence  in  this 

region,  S(y)  can  be  expanded  as 
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^(y)  =  a'ln(i;)+— +-^+....  (8) 

The  constants  a ,  h ,  c  etc.  can  now  be  calculated  from  the  experimental  data  for  large  v  using  a 
fit  of  the  form  of  (8)  to  d(v).  Since  a  in  (1)  is  related  to  the  drag  coefficient  for  the  target 
medium,  a  is  a  physical  parameter  which,  using  Equations  (1)  and  (8),  can  be  related  to  the 

coefficient  of  a'  of  the  leading  term  in  (8).  This  yields  a  =  — . 

a' 

2.11.3  Solution  in  the  Intermediate  Velocity  Region 

In  the  intermediate  region  of  the  penetration  velocity  v,  the  solution  for  F(v)  and  the  analytical  fit 
for  S(v)  can  be  found  by  using  matched  asymptotic  expansions  of  these  functions  in  the  low 
and  high  velocity  range.  This  procedure  is  discussed  in  more  detail  in  Section  2. 1.1.9.  The 
matched  asymptotic  technique  is  based  on  the  simple  rule  that  there  is  an  overlapping  region  in 
the  intermediate  velocity  range  where  both  the  low  velocity  (Equation  6)  and  high  velocity 
(Equation  8)  expansions  yield  the  same  result.  This  procedure  is  discussed  and  is  illustrated  by  the 
example  discussed  in  Section  2. 1.1. 5  below. 

2. 1.1. 4  Final  Polynomial  F\i\oF(v) 

After  the  determination  oiF(v)  is  complete  using  the  procedure  given  above,  it  is  possible  to  fit  a 
polynomial  in  v  to  fit  the  resulting  The  virtue  of  such  fit  can  be  realized  in  the  theoretical 

prediction  of  the  penetration  depths  for  velocities  in  the  interior  of  or  outside  the  experimental 
database.  Since  the  calculation  of  6(v)  from  a  fitted  F(v)  involves  an  integration  (which  is  a 
smoothing  process)  as  in  Equation  3,  estimates  of  6(v)  are  usually  very  good.  These  expansions 
are  however  not  valid  for  characterization  of  S(v)  using  the  differential  relation  in  (4)  as  they 

may  predict  the  incorrect  behavior  of  these  functions  at  both  low  and  high  velocity  domains  For 
example,  if  the  polynomial  fit 


(9) 

is  used  to  predict  the  behavior  of  d(v)  for  high  velocity  ranges  using  (4),  we  get 

S{v) »  ^  and  F{v) «  (10) 

Equations  10  contradict  Equation  8  and  the  second  equation  in  (1)  respectively. 

2. 1.1. 5  An  Example 

For  the  case  of  a  target  made  of  20%  gelatin,  we  have  gathered  experimental  data  on  penetration 
depth  and  entry  velocity  for  small  spherical  projectiles.  These  projectiles  are  composed  of  steel 
with  a  mass  of  about  0.344  grams  and  a  diameter  of  1.0  cm  and  density  of  8.46  grams/cc, 
approximately.  The  entry  velocity  has  been  varied  from  340  ft/sec  to  6,000  ft/sec  and  the  total 
penetration  depth  was  measured  in  inches.  A  schematic  of  these  experiments  is  shown  in  Figure 
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1  and  the  experimental  data,  with  supplementary  experimental  data  from  the  Edgewood  Arsenal, 
are  shown  in  Figure  2. * 

2. 1 . 1 .6  Data  Processing 

2. 1 . 1 .6. 1  Low  Velocity  Region 

As  indicated  previously,  in  this  range  we  first  compute  6(v)/v^  versus  v  from  the  known 
experimental  data  prior  to  a  polynomial  fit  to  obtain  S(v)  in  the  form 

=  (11) 

i=0 

where  n  is  the  degree  of  the  polynomial  and  a's  are  unknown  constants  to  be  determined  from  the 
polynomial  fit  to  the  experimental  data.  For  n=2,  we  find 

Oq  =  5.13  X 10"^,  Oi  =  -2.844  x  10"^,  02  =  5.0  x  10"'^ . 

2. 1 . 1 .6.2  High  Velocity  Region 

For  high  velocities,  d(v)  behaves  like  ln(v)  and  hence  can  be  expanded  as 

n 

S{v)  =  o'  ln(t>)  +  b{v~^  (12) 

i=l 

The  a'  and  b'.  's  can  be  calculated  from  fitting  (12)  to  the  experimental  database  for  high 

velocities.  This  yields  o'  =  1261,  bl  =  -7.67  x  10^ .  Comparisons  between  the  low  and  high 
velocity  asymptotes  and  the  experimental  data  are  shown  in  Figure  3.  As  indicated  previously, 
there  is  an  intermediate  region  of  velocity  where  both  of  these  asymptotes  yield  the  same 
penetration  depth.  For  this  experimental  database,  this  region  is  seen  to  lie  between  2000  ft/sec  to 
2200  ft/sec  (Figure  3). 


’  R.  D.  Eisler,  et  al,  Phase  11  Development  of  a  Soldier  Protective  Ensemble  Computer  Aided  Design  (SPE/CAD) 
System,  Final  Report  for  U.S.  Army  Natick  Contract  DAAK60-C-92-0(X)8,  Mission  Research  Report  MRC-COM- 
94-R-0455, 4  July  1997,  To  Be  PuWished  as  U.S.  Army  Soldier  System  Command  Technical  Report. 

^  BRL:  Report  No.  2423, 1974. 

^  Huelke,  D.  F.  et.  al,An  experimental  study  in  bio-ballistics:  Femoral  fractures  produced  by  projectiles-II,  Shaft 
Impacts,  J.  BIOMECHANICS,  Vol.  1,  pp.  313-321,  1%8 

^  R.  D.  Eisler  and  A.  K.  Chatteijee,  Effect  of  Coatings  and  Processing  on  Failure  of  Polycarbonate  Substrates 
During  Ballistic  Impact,  Final  Report  for  Natick/MRC  contract  DAAK60-91-C-0087,  MRC  report  MRC-COM-R- 
97-500,  To  Be  Published  as  Natick  Teclmical  Report,  3 1  December  1997. 

^  MRC  conducted  an  extensive  literature  search  to  acquire  data  relative  to  bone  trauma  produced  by  ballistic  events. 
Less  than  twelve  sources  were  identified  since  1900.  The  most  recent  data  is  from  the  People’s  Republic  of  China 
during  the  inid-1980’s.  In  general,  the  data  from  these  experiments  is  compromised  by  the  relatively  imsofdiisticated 
instrumentation  available  at  the  time.  We  submitted  a  proposal  to  DARPA  in  conjunction  with  Maj.  Paul  Dougherty, 
M.D.,  USAMC  to  conduct  ballistic  tests  on  femurs  embedded  in  gelatin  matrices  with  in-situ  instrumentatioa  In 
addition  to  the  in-situ  instrumentation,  the  bones  would  be  tested  in  a  fixture  that  would  impose  tendon  and  gravity 
loads  representative  of  a  standing  soldier.  Although,  this  proposal  was  selected  for  fimding,  the  overall  Government 
program  fimding  these  activities  ended  and  hence  this  proposal  remained  in  abeyance. 
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Figure  2.  MRC  cmd  Edgewood  Arsenal  Data  for  Penetration  of  Spherical  Projectiles  into  20% 
Ordnance  Gelatin  (Reference  1) 

2. 1.1. 7  Determination  of  Retardation  Force 

Using  the  low  and  high  velocity  asymptotes  of  ^(v)  given  in  Equations  1 1  and  12,  respectively, 
and  the  relation  between  F(v)  and  dip)  given  in  Equation  4,  we  can  determine  the  retardation 
force  F(v).  For  the  current  database,  the  computation  of  F(v)  is  shown  in  Figure  4  by  the  solid 
line  curve.  We  fit  polynomials  of  degrees  2, 3  and  4  to  the  theoretical  F(v)  distribution  and  these 
results  are  shown  in  Figure  5.  The  coefficients  of  these  polynomial  fits  are  tabulated  in  Table  1. 
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Figure  3.  Comparison  Between  Experimental  Penetration  Data  and  High  and  Low 

Velocity  Asymptotes 


Table  1.  Polynomial  Fit  to  the  Retardation  Force,  F(v)  Function 


Degree  of 
Polynomial 
Fit  to  F(v) 

Coefficients  in  ascending  order  of  the  variable 

2 

2.45x10^ 

-570.63 

0.044 

3 

0.95x10® 

106.46 

2.287x10"® 

7.04x10"® 

4 

0.88x10® 

195.33 

-0.16 

1.48x10"^ 

1.07x10"® 
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Figure  5.  Comparison  of  Experimental  Penetration  Data  and 
Theoretical  Results 


Once  we  have  determined  the  theoretical  distribution  of  the  retardation  function  F(v),  we  can 
now  calculate  the  theoretical  penetration  depth  function,  6(v),  so  that  the  accuracy  of  the 
theoretical  calculations  can  be  estimated.  Using  (3)  and  the  polynomial  fits  given  in  Table  1,  we 
have  compared  the  experimental  data  with  the  theoretical  prediction  of  the  penetration  function. 
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The  results  are  shown  in  Figure  5.  From  Figure  5,  we  see  that  the  cubic  polynomial  fit  of  the 
retardation  function,  F(v),  gives  the  best  result  wWIe  the  quadratic  fit  yields  the  worst  result.  To 
improve  the  fit  of  the  quadratic  form  of  the  retardation  function  we  have  adjusted  the  coefficients 
of  the  second  row  of  Table  1  so  as  to  yield  the  same  experimental  penetration  depth  at  the 
maximum  velocity  of 3400  ft/sec  in  the  current  database  from  spherical  projectiles  into  20% 

gelatin.  The  new  coefficients  to  the  quadratic  fit  to  F(v)  are  4.046x10^,  -941.54, 0. 72.  The 
resulting  <5(V;  from  this  fit  is  compared  with  the  experimental  data  in  Figure  6. 


2. 1.1. 8  Application  to  Flechette  Penetration  Database 

We  have  experimental  results  on  the  penetration  of  flechettes  into  20%  gelatins  from 
Natick/MRC  contract  DAAK-92-C-0003.  As  part  of  this  effort  MRC  developed  numerical 
algorithms  which  can  describe  the  two-dimensional  motion  of  a  flechette  when  the  quadratic 
form  of  the  retardation  force  is  prescribed.  To  estimate  the  accuracy  of  these  algorithms,  which  is 
based  on  modeling  the  projectile  as  a  rigid  ensemble  of  particles  where  the  resistance  to  particle 
penetration  is  derived  from  the  retardation  functions  derived  from  small  spherical  projectiles,  we 
have  applied  the  above  estimated  polynomial  fit  to  the  retardation  force  function  as  input  to  our 
flechette  code.  Unfortunately,  at  the  time  of  development  of  the  flechette  code,  we  assumed  that 
a  good  quadratic  fit  is  possible,  and  hence  a  quadratic  retardation  force  and  local  velocity  is 
assumed  in  the  code.  Since  it  is  resource  intensive  to  modify  our  flechette  code  to  incorporate 
higher  order  polynomial  fits,  we  decided  to  use  the  modified  coefficients  for  the  quadratic,  F(\) 
and  V  (Figure  6)  in  our  flechette  code  to  obtain  and  compare  the  theoretical  predictions  to  the 
experimental  flechette  database.  Figure  7  compares  results  from  a  ballistic  test  for  a  flechette 
versus  predicted  results.  The  results  agree  reasonably  well  considering  the  fact  that  we  have  not 
used  a  good  fit  to  the  retardation  function  F(v).  If  the  current  flechette  code  were  modified  to 
accept  higher  order  F(v)  and  v  polynomials,  we  would  expect  much  better  agreement  between 
experimental  data  and  the  theoretically  predicted  results. 
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Figure  7.  Comparison  of Ballistic  Data  versus  Predicted  Results  for  a  19. 6-grain  flechette  striking 
20%  ordnance  gelatin  at  1200 ft/sec  with  a  6.6  degree  yaw 


2.1 .1 .9  Calculation  of  Retardation  Coefficients  from  Penetration  Data 

As  concluding  remarks  to  the  above  analysis,  we  would  like  to  point  out  some  important  aspects 
of  the  theoretical  modeling  of  penetration  problems.  First,  in  the  absence  of  basic  understandings 
of  the  physics  for  low  and  high  entry  velocities,  it  is  not  possible  to  model  the  nature  the 
retardation  force  from  experimental  data.  In  this  regard.  Equation  4  should  be  used  as  a  guiding 
equation  to  connect  the  characters  of  the  two  fundamental  functions;  d(v),  the  penetration  and 
entry  velocity  function  and  F(v),  the  retardation  force  and  local  velocity  function.  Special  care 
should  be  taken  when  we  use  the  derivative  of  6(v)  to  obtain  F(v),  and  large  errors  may  be 
introduced  if  the  proper  functional  forms.  Equations  7  and  8,  are  not  used. 

Secondly,  the  usual  quadratic,  polynomial  form  of  the  retardation  force  function,  F(v)  and  v,  does 
not  yield  acceptable  solutions  when  compared  with  the  experimental  data,  and  correct  forms 
should  be  carefully  incorporated  in  the  analysis.  It  is  interesting  to  mention  here  that  a  quadratic 
form  is  tempting  because  it  predicts  correct  behaviors  at  both  high-  and  low-velocity  regions.  But 
our  problem  here  is  to  come  up  with  a  functional  form  that  can  be  used  for  intermediate 
velocities  as  well.  This  is  exactly  where  such  a  functional  form  fails  to  provide  the  correct 
results.  Since  6(v)  depends  on  the  integral  of  F(v),  error  accumulation  coming  from  the 
intermediate  local  velocities  are  manifested  in  the  incorrect  prediction  of  the  final  penetration 
depth,  6(v),  for  a  given  entry  velocity  v. 

In  our  analysis  we  have  used  the  assumption  that  the  force  of  retardation  per  unit  mass,  F(v), 
during  penetration  of  gelatin,  can  be  expressed  as  a  quadratic  function  of  the  projectile 
instantaneous  velocity,  v,  i.e.,  F(v)  =  av2  +  y,  where  a,  p,  and  y  are  undetermined 
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coefficients.  From  previous  reduction  of  experimental  data  it  can  be  shovm  that  a  is  related  to 
the  high  velocity  behavior  of  the  penetration  depth  function,  Sfv), 

<5('')w=a'ln(v)+-+... 

V 

where  v  is  the  entry  velocity  while  y  is  related  to  S(v)  at  low  velocities. 


w  =v\a  +  bv  +  cv^  +. . .) 

High  velocity  refers  to  about  90%  of  the  maximum  velocity  of  interest,  Vmax  (~  4,500  ft/sec  for 
MRC  small  spherical  projectile  penetration  data)  where  S(v)  behaves  like  ln(v).  Although,  ln(v) 
is  the  first  order  and  dominant  term  of  the  series  at  high  striking  velocities,  S(y)  high,  there  are 
second  order  terms  that  include  contributions  from  the  linear  and  constant  terms  in  the  quadratic 
retardation  force  function.  This  is  not  the  case  for  the  low  velocity  behavior  of  S(y)  however. 
These  results  were  established  earlier  from  the  exact  functional  relationship  between  F(v)  and 
S(v),  i.e.,  S'(v)=v/F(v)  where  prime  indicates  differentiation  with  respect  to  the  argument  v. 

Using  previously  generated  experimental  data,  the  coefficients  above  have  the  following  values 
for  small  spherical  projectile  penetration  into  20%  ordnance  gelatin. 

Low  Velocity  Coefficients 
a  =  5.13E-6 
b  =  -2.84E-9 
c=  5E-13  . 

High  Velocity  Coefficients 
a' =  1.2.61 
b'  =  -7.67E3 

Depending  on  the  choice  of  the  offset.  Sc,  is  constant  and  be  estimated  by  writing  the  high 

velocity  asymptote,  Sfvjia^,  as  d(v)=ln(v)/a  -  Sc.  In  the  following  calculations,  v  is  measured  in 
ft/sec  while  lengths  are  measured  in  inches.  From  the  penetration  data  on  small  diameter 
spherical  projectile  data  into  gelatins.  Sc  is  approximately  2  inches. 

The  P  term  cannot  be  directly  inferred  since  it  represents  contributions  from  the  various  phases 
of  the  gelatin,  which  predominate  at  intermediate  velocities.  The  following  procedure,  which  can 
only  be  applied  to  the  high  velocity  asymptote  which  include  secondary  contributions  from  other 
retardation  terms,  is  therefore  used  to  estimate  p.  For  any  given  velocity,  v  =  vh,  in  the  high 
velocity  region  of  the  parameter  space,  calculate: 


a’\n(y)-S(ytj)\^=Sc 

This  forces  the  parameter  to  fit;  a'  ln(v)  -  S^  for  v  =  vh. 

The  error,  E(P),  for  any  given  vh  can  therefore  be  estimated  from 
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=  I  ^  )  -  ^c] 

Our  objective  is,  to  find  P  for  a  known  a  and  y  such  that  the  root  of  E(3)  =  0.  The  coefficients  of 
the  retardation  function  have  been  estimated  fi'om  the  experimental  data  using  two  approaches. 

Estimated  values  for  a  and  y  using  fits  to  the  parametric  forms  of  the  high  and  low  velocity 
a  =  -^  =  0.793  r  =  —  =  0.974E5 

asymptotes  are  and  2a  Coefficients  estimated  fi'om  fitting  a 

quadratic  retardation  function  to  the  entire  set  of  experimental  data  (which  as  discussed  in 
previous  progress  reports  yields  poor  estimates  of  penetration  depth)  results  in  the  following 
coefficients  (distinguished  by  the  subscript  “c”). 

ac  =  0.8637  pc  =  1408.68  Yc  =  0.687E6 

In  both  cases  v^ax  =  4500  ^s. 

The  first  term  in  the  expression  for  E(P)  above  is  identically  equal  to  5(v)  as  described  in 
previous  progress  reports.  Assuming  a  quadratic  form  of  the  retardation  function.  8(v)  can  be 
integrated  as  follows. 


^(v)  =  . 


ccv^  -  Pv  +  Y 


la 


A{P)  + 


P 


ccMP) 


arct 


lav-  p 


+  arctanj 


y^p) 


The  first  term  in  the  expression  above  can  be  simplified  to 


Knowing  that  E{P)  =  ^(v)  -  +  5^  and  using  5c  =  2  inches  which  is  the  empirically 

a 

observed  offset  fi'om  the  logarithmic  curve  at  high  striking  velocities,  we  can  determine  E(P)  =  P 
as  shown  Figure  8. 


The  optimum  P  which  minimizes  the  error  term,  E(P)  =  0,  \s  p  =  -478.54. 
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P 

Figure  8.  E(beta)  versus  beta 


Figure  9.  Comparison  of  Penetration  Depth  versus  Striking  Velocity  using 
Empirical  Fit  to  Quadratic  Function  versus  Semi-Analytical  Approach 
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Previous  experimental  data  fit  to  a  quadratic  retardation  function  yields.  Pc  =  1408.68.  A 
comparison  between  both  approaches  and  estimated  penetration  depth  is  shown  in  27  where 
discrepancies  accrue  for  a  significant  portion  of  the  intermediate  velocities. 

At  high  velocities  both  approaches  agree  with  the  data  and  each  other,  6(4500 Jps)  =  8.608 
inches  and  Sc(4500 Jps)  =8.517  inches.  To  check  the  estimation  procedure  for  P,  the  process  was 
repeated  to  estimate  Pc.  For  this  case  when  a  =  Oc  and  y  =  yc  the  results  in  Figure  10  are  produced 
for  E(Pc)  versus  Pc,  see  Figure  10. 

The  actual  root  of  {E(P),  P)  =  0  corresponds  to  P=  1.415E3.  This  compares  to  the  empirically 
determined  Pc  =  1.409E3.  This  appears  to  corroborate  the  procedure  to  estimate  coefficients  of 
terms  in  the  retardation  function  that  predominate  at  intermediate  velocities. 

The  next  step  in  the  analysis  was  to  determine  if  there  was  continuity  in  the  slope  of  the  high  and 
low  velocity  asymptotes  and  determine  sensitivity  of  these  slopes  to  the  specific  value  of  the 
coefficients.  This  is  a  more  stringent  requirement  on  the  fidelity  of  the  retardation  function  than 
simply  ascertaining  the  coefficient  of  an  assumed  polynomial  coefficient.  In  order  to  have 
continuity  of  the  first  derivative,  the  functional  representation  must  be  much  more  highly 
resolved.  For  low  velocities  where  we  assume  a  threshold  velocity  of  v/  =  0.  Iv^ax'. 

c+bv^  =  av^  +pvi+y  (a) 

b  =  2av,  +p  (b) 

and  for  high  velocities  where  we  assume  a  threshold  velocity  of  V2  =  0.9vmax 

avl  =  a\>l  +Pv^+Y  (c) 

2av2=2av2+>9  (d) 

If  we  normalize  the  striking  velocities  to  Vmax  ie.,  V;  =0.1  and  vj  =  0.9  and  assume  an  order  of 
magnitude  difference  between  a  (=  10.0)  and  c  (=  1.0)  then  the  following  applies. 


.2 

vl 

vl 

1 

vl 

1.25 

-0.125-1.25 

1.125 

2vl 

1 

0 

-1 

A-»  = 

-2.25 

0.225  2.25 

-1.025 

A  = 

v2^ 

v2 

1 

0 

1.013 

-0.101-0.012 

0.011 

2v2 

1 

0 

0 

-2 

-0.8  2 

-0.8 

a 

c 

p 

:=A-1- 

0 

Y 

a-v2^ 

b 

2av2 

24 


Mission  Research  Corporation 


DAMD17-94-C-4099 


Figure  10.  E(betac)  versus  betac 


Figure  11.  Retardation  Function,  F(V),  with  Low  and  High 
Velocity  Asymptotes 


form  which  we  can  obtain: 


a=  11.375,  p  = -2.475,  y=  1.114,  and  b  = -0.2. 


Re-doing  this  problem  with  threshold  velocities  assumed  to  be  v;  =  0.01  and  V2  =  0.99 
result  in: 

a  =  1 1 .021,  P  =  -2.022,  y  =  1 .001,  and  b  =  -1 .802. 
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A  plot  of  this  retardation  function  with  the  low  and  high  velocity  asymptotes  is  shown  in  Figure 
11. 

2.1.1.10  Summary 

Predicting  the  geometry  of  wound  tracts  in  humans  requires  a  retardation  function  that  describes 
the  velocity  history  of  the  projectile  during  penetration  in  human  tissue.  This  velocity  history  is 
used  as  an  initial  condition  to  other  models  that  describe  tissue  displacement  normal  to  the 
wound  tract.  Since  ballistic  testing  caimot  be  done  on  humans,  ballistic  testing  is  implemented  on 
surrogate  materials,  which  phenomenologically  manifest  behavior  similar  to  what  would  be 
expected  in  human  tissue.  A  model  is  then  developed  for  the  surrogate  and  assumed  to  have  a 
functional  form  similar  to  human  tissue.  The  model  parameters  however  need  to  have  rigorous 
physical  interpretations  so  that  they  can  be  ascertained  independent  of  a  ballistic  test.  In  our  case 
we  need  to  be  able  to  quantify  the  model  parameters  based  on  tissue  properties  that  are  available. 

The  procedure  discussed  in  this  section  enables  the  determination  of  p  for  known  a  and  f  where 
a,  p,  and  fare  coefficients  in  the  polynomial  representation  of  the  retardation  function.  The 
retardation  function  describes  the  velocity  history  of  the  projectile  during  penetration  and  is  an 
essential  element  of  the  methodology  to  predict  the  wound  tract  in  humans.  The  quadratic  term 
of  the  retardation  function  has  been  shown  to  be  associated  with  the  dynamic  pressure  produced 
by  passage  of  the  projectile  and  is  significant  at  high  velocities  (see  Figure  16).  Given  this 
interpretation,  the  only  tissue  property  that  is  relevant  for  calculating  this  term  in  the  retardation 
function  is  tissue  density  or  specific  gravity.  Other  properties  that  are  relevant  are  associated 
with  the  projectile  geometry  and  striking  conditions.  There  is  evidence  to  suggest  that  second 
order  terms  effects  can  be  important,  which  is  the  reason  for  the  offset,  8c  from  the  logarithmic 
term  in  the  plot  of  penetration  depth  versus  striking  velocity.  These  are  all  properties  that  can,  for 
the  most  part,  be  determined  independent  of  ballistic  testing.  This  means  we  can  fairly  reliably 
predict  the  value  of  this  model  parameter  applicable  to  human  tissue. 

Similarly,  the  constant  term,  y,  can  be  shown  to  be  equal  to  [ Uc  ApJ/nip,  where  ctc  is  a 
characteristic  stress  in  the  tissue,  Ap  is  the  projected  area  of  the  projectile,  and  nip  is  the  projectile 
mass.  If  we  assume  that  the  characteristics  stress  is  the  ultimate  tensile  strength  of  the  target 
material,  which  in  the  case  of  gelatin  is  between  10  and  100  psi  and  low  strain  rates,  for  the 
ballistic  testing  that  has  been  accomplished  employing  a  341  mg,  0.171  diameter  steel  ball,  we 
obtain  an  estimate  for  y  that  is  somewhere  between  0.4E5  and  0.4E6.  Recall  that  this  compares 
with  y  =  0.974E5  based  on  calculations  using  S(y)  and  yc  =  0. 687E6  using  a  quadratic  fit. 

The  detailed  quasistatic  stress  strain  curve  for  the  different  formulations  of  gelatin  can  be 
compared  with  the  differences  in  penetration  depth  observed  at  the  low  velocities. 

The  detailed  stress  strain  curve  is  also  required  for  the  transient  response  model.  This  model 
describes  tissue  displacement,  strain,  and  stress  normal  to  the  wound  tract.  Fortunately, 
quasistatic  stress  strain  data  is  available  for  many  human  tissues  from  fresh,  unembalmed, 
cadavers.  This  would  allow  model  parameters  to  be  calculated  for  penetrating  injuries  in  humans. 
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2.2  Rotational  Kinematics  of  a  Projectile  inside  a  Viscoelastic  Fluid 

Bodies  penetrating  a  highly  viscous  material  suffer  retardation  forces  mainly  originating  from 
three  fundamental  sources.  These  are:  (a)  a  dynamic  pressure  drag  force  which  is  usually 
proportional  to  the  second  power  of  the  instantaneous  velocity,  v,  of  the  contact  point  between 
the  projectile  and  the  target,  (b)  a  viscous  drag  force  proportional  to  the  magnitude  of  v,  and  (c)  a 

drag  independent  of  v.  The  last  term  is  basically  the  force  needed  to  fracture  the  material  in  the 
appropriate  mode. 

From  a  dimensional  point  of  view,  independent  dimensions  of  both  ^P^  and  v  can  be  found  from 
the  phy  sical  properties  of  the  target  and  the  projectile  and  hence  any  power  of  such  non- 
dimensional  groups  can  be  deemed  to  influence  the  aoxb  retardation  terms.  Obviously,  the 
simplest  function  of  these  non-dimensional  groups  is  a  linear  function  of  the  form  a\^+bv+c 
where  a,  b,  and  c  are  constants  depending  on  the  size,  shape  and  other  material  properties  of  the 
projectile  arid  the  target.  Determination  of  these  constants  for  a  specific  case  is  however 
e.^tremely  difficult  since  the  general  motion  of  the  projectile  is  at  least  two-dimensional  where 
the  projectile  center  of  mass  moves  with  a  superposed  motion  about  its  center  of  mass. 

For  some  applications,  extensive  experiments  have  been  conducted  to  determine  the  depth  of 
penetration  as  a  function  of  the  incident  velocity  of  the  projectile  for  specific  combinations  of 
targes  and  projectiles.  These  experiments  are  useful  since  such  data  can  be  directly  inverted  to 
obtain  the  functional  form  of  the  force  of  retardation  on  the  motion  of  the  center  of  mass.  Since 
the  motion  of  the  center  of  mass  depends  only  on  the  vector  sum  of  all  the  forces  acting  on  a 
rigid  body,  we  can  use  this  information  to  estimate  the  force  on  a  point  mass  moving  in  the  same 
media  when  the  instantaneous  velocity  of  the  point  mass  is  known.  Using  this  concept,  a 
mathematical  model  was  developed  which  describes  the  two-dimensional  motion  of  projectile 
inside  the  target  when  the  initial  conditions  are  prescribed.  This  motion  is  divided  into  two  parts; 
motion  of  the  center  of  mass  and  the  motion  about  the  center  of  mass.  Since  the  order  of  these 
two  motions  is  immaterial  for  a  rigid  body,  these  two  motions  completely  determine  the  location 
of  the  projectile  inside  the  target  as  a  function  of  time. 

The  motion  of  the  center  of  mass  {CM)  of  a  projectile  entering  a  target  can  be  calculated  by 
integrating  the  differential  equation  of  motion  in  the  velocity-distance  space  under  the 
assumption  that  the  coefficients  a,  b,  c  in  the  retardation  force-velocity  relation  are  known  from 
experimental  or  theoretical  calculations.  These  coefficients  represent  average  values  that  are 
capable  of  predicting  the  motion  of  the  CM  only.  To  extend  the  results  to  foil  two-dimensional 
prediction  of  projectile  motion  inside  a  target,  we  need  to  make  some  assumptions. 

One  of  the  most  important  assumptions  is  that  in  the  current  mathematical  formulation,  the 
associated  coefficients  are  also  capable  of  estimating  the  force  of  retardation  on  a  point  mass 
mowng  inside  a  target.  If  this  is  true,  then  a  rigid  body  may  be  conceived  as  an  ensemble  of 
particles  where  the  ret^dation  force  on  each  particle  is  known,  and  hence,  the  equations  of 
motion  governing  motion  about  the  CM  can  be  written  down.  If  we  ignore  the  fiictional  drag  on 
these  particles,  then  the  resistive  force  can  be  calculated  by  the  same  functional  form  shown 
above  except  that  v  should  be  replaced  by  the  normal  component  of  the  velocity  of  a  particle 
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Figure  12.  Projectile  location  and  loading  at  time  t 


mass  on  the  projectile.  Here  the  normal  direction  is  defined  by  the  outward  normal  direction  at 
the  projectile  surface. 

Assuming  that  the  force  of  retardation  R(v)  per  unit  mass  on  each  of  these  collection  of 
particles  constituting  the  projectile,  we  can  write  R(v)  as 

R{v)  =  av^  +  bv  +  c  (13) 

where  the  v  is  the  normal  component  of  the  velocity  as  described  above.  The  distribution  of  these 
forces  on  the  body  and  fin  area  of  a  flechette  is  shown  in  Figure  12. 

At  any  given  time  t,  let  [xQ(t),  yG(t)J  and  Q(t)  be  the  coordinates  of  the  center  of  mass  G  of  the 
projectile  and  the  angular  coordinate  defining  the  inclination  of  the  projectile  body  with  they- 
axis.  These  coordinates  completely  specify  the  location  of  the  projectile  at  time  t.  In  order  to 
determine  these  coordinates;  we  need  to  use  the  equations  of  motion  of  and  about  the  CM  Since 
the  motion  of  the  CM  is  governed  by  a  corresponding  particle  mass  problem  where  the  mass  of 
the  particle  is  equal  to  the  mass  of  the  flechette  and  all  the  forces  on  the  flechette  are  assumed  to 
pass  through  the  CM,  we  have  the  following  equations  of  motion  for  the  projectile  CM. 

I 

mxa  =  -J/?f?(|up;i|)sgn(i;^)cos^  -myi?(|ye„|)sgn(i;e„)sin^ 

-t 
l 

mya  =  +  jA^^(|ypn|)sgn(i;p„)sin^dr  -  myi?(|i;e„|)sgn(i;e„)cos^ 

-l 


(14) 
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The  equation  of  motion  about  the  CM  can  be  obtained  by  taking  moments  of  the  equations  of 
motion  of  the  ensemble  of  particles  constituting  the  projectile,  about  the  CM.  This  yields 

/ 

lG^=-\pR^Vpn\sQr\{Vpri)rdr  (15) 

-l 

In  Equations  14  and  15,  p  is  the  mass  of  the  projectile  per  unit  length,  r  is  measured  from  the  CM 
(Figure  12),  mj  is  the  mass  of  the  base  of  the  projectile,  21  is  the  projectile  and  Iq  is  the  moment 

of  inertia  of  the  projectile  about  its  CM.  R  is  the  force  of  retardation  per  unit  mass  for  the  particle 
mass  problem. 

The  velocity  of  a  particle  P  at  a  distance  r  from  the  CM  in  the  direction  normal  to  the  flechette 
body  can  be  expressed  in  terms  of  the  velocity  components  of  the  CM  in  the  x-  and^'-directions 
and  angular  velocity  6 ,  as 

Vp,i=r0+kQCOs9-yQS\n0  (16) 

Using  Equations  (13)  and  (16),  P  can  be  written  as  in  terms  of  the  retardation  coefficients  a,  b,  c 
as 


■^(|^pn|)  =  +  {2aA ± b)9r  +  {afi}  ±bA  +  c) 

(17a) 

where  A  -  XQCOs9-yQS\r\9 

(17b) 

Similarly,  i?(|t'e„|)  is  given  by 

R<\ven\l  =  oB^  ±bB  +  c 

(18a) 

where  B  =  XQS\r\9+yQCXis9 

(18b) 

Substituting  Equations  17  and  18  into  Equations  14  and  15,  and  integrating,  it  follows  that 


{mi,+mf)xQ 


-  El  cos  0  -  mjr2  sin  ^ 


{mh+mf)yQ  = 


21 


Fj  sin  0  -  m;r2  cos  9 


Ig 


■9=-^r, 
21  ^ 


(19) 


In  (19),  r 's  are  defined  as 
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Ti  =  +  2A0I^}]s  +  (oA^  +  +  2lbA 

r2  =  {oEp-  +  c)  sgn(B)  +  bB 

Ts  =  +  (0^2  +0)75)5 

(20) 


1  ,3 


For  the  projectile,  Iq  is  given  by  Ig  =  -mi^l  +  m fl 


In  (20),  75b  ,n  =  0,l,2,3axe  integrals  that  keep  track  of  load  reversal  due  to  the  change  in  the 

location  of  the  instantaneous  center  of  rotation.  In  order  to  explain  the  cause  of  this  load  reversal, 
the  expression  that  calculates  the  normal  component,  Vp^ ,  of  velocity  (see  Figure  12)  needs  to 

be  considered.  As  shown  earlier,  Vp^  is  given  by  (1 6).  To  find  the  regions  of  positive  the 

following  four  cases  may  occur  depending  on  the  relative  values  and  signs  of  A  and  9 .  With 
reference  to  Figure  12,  Vp^  is  positive  in  the  following  regions. 


L  ril  =  -jl<r<I,  forA<0,  ^>0 


2.  -r2^—\<r<l,forA>0,d>0 


3.  -  Z  <  r  <  rJ  =  — -\  <r  <1,  for  A  >0,  0<O 


4.  -l<r< -r2|^=  — J <r  <1,  forA< 0,  0<O 


With  the  above  definitions  of  rj  and  r2 ,  the  integrals  I  can  be  expressed  as 


75*)  =  mi  r^dr  +  m2ffdr 

tn+l  ^n+1 


-(m2  +(-!)” mi) +  ^^(-l)”'*'Hmi  -m2) 


7^*)  =  mi  J  ‘  r^dr  +  m2  J  r^dr 


tn+l  n+1 

- -(m2  +(-!)” mi)  +  -!—-(-l)”+*(mi  -m2) 

n  +  1  n  +  1 


In  Equation  22,  mi  and  m2  are  either  +7  or  -7  and  defined  as  follows.  For  the  75*)  integral  (i.e., 
when  Ti  is  feasible  which  corresponds  to  cases  2  and  4  of  Equation  21),  mi  =  -sgn(A)and 
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m2  -  sgn(^) .  Similarly,  for  the  integral  (cases  1  and  3  of  Equation  21),  m|  =  sgn(A)and 

/Wj  —  —  sgn(^) .  These  integrals  keep  track  of  the  reversal  of  the  retardation  forces  when  the 
location  of  the  instantaneous  center  of  rotation  changes  dynamically  with  the  relative  values  of  A 
and  d. 

Equations  19  through  22  completely  determine  the  two-dimensional  motion  of  the  projectile. 
These  equations  are  coupled,  nonlinear,  second  order  differential  equations  that  have  no  closed 
form  solutions. 

2.2.1  Initial  Conditions 

To  solve  the  system  of  differential  equations  above,  initial  conditions  specifying  the  location  of 
the  center  of  mass,  G,  and  the  inclination  of  the  projectile  at  time,  /  =  0,  need  to  be  defined.  The 

initial  velocity  {xQ,yQ)  of  the  CM  and  the  angular  velocity,  9 ,  of  the  projectile  also  need  to  be 
specified. 

When  the  initial  conditions  above  are  specified,  the  differential  equations  of  motion.  Equations 
21  and  22,  can  be  solved  using  a  fourth  order  Runge-Kutta  numerical  integration  routine. 

2.2.2  Motion  of  a  Partially  Embedded  Flechette  inside  a  Target 

Equations  14  through  22  describe  the  motion  of  a  projectile  that  is  fully  embedded  in  a  target  at 
all  times.  However,  prior  to  full  penetration,  the  projectile  undergoes  a  phase  of  partial 
penetration.  During  this  phase,  the  limits  of  the  integrals  in  Equation  14  should  be  replaced  by 
proper  limits  that  define  the  part  of  the  projectile  inside  the  target.  Other  than  these  changes,  the 
methodology  used  above  are  applicable  for  the  initial  phase  of  partial  penetration  prior  to  the 
problem  of  the  fully  embedded  projectile  moving  in  a  target  discussed  above. 
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2.3  Analysis  of  projectile  interaction  with  Soft  Tissue 


2.3.1  Determination  of  Temporary  and  Permanent  Cavities  and  Nonlinear  Spring-Mass 
Model 

We  have  developed  a  computer  code  that  calculates  the  transient  response  of  an  elastic  medium 
due  to  a  moving  projectile.  This  problem  is  solved  by  simulating  the  actual  physical  problem  of 
determining  the  transient  response  of  a  medium  when  ballistic  penetration  takes  place.  This 
problem  is  solved  by  the  superposed  response  of  properly  distributed  transient  point  sources.  We 
have  used  the  fundamental  Green's  function  solution  for  a  transient,  stationary  point  source  with 
an  arbitrary  temporal  distribution  as  the  building  block  for  the  desired  solution.  The  fundamental 
solution  has  three  terms  each  signifying  the  response  for  various  time  periods  at  a  given  point 
inside  the  medium.  Physically,  the  mathematical  solution  represents  the  sum  of  the  near  field 
solution,  far  field  P-wave  and  far  field  S-wave  solutions.  The  desired  solution  of  the  above 
problem  is  then  constructed  fi'om  this  fundamental  solution  by  properly  distributing  a  series  of 
sources  along  the  penetration  path  of  the  projectile.  The  temporal  distribution  for  each  of  these 
transient  sources  has  a  support  function  that  is  equal  to  the  travel  time  of  the  moving  projectile 
fi'om  one  source  location  to  the  next  one.  Adjusting  the  spacing  between  the  point  sources  tests 
solution  convergence.  From  our  earlier  results  on  the  nature  and  magnitude  of  the  drag  forces 
when  the  instantaneous  velocity  of  the  projectile  is  known,  the  intensity  of  a  particular  source  is 
calculated  as  a  fimction  of  location  and  time.  The  intensity  depends  on  the  velocity  of  the 
projectile  when  it  is  passing  through  the  transient  point  source  location. 

2.3.2  Calculation  of  Transient  Response  of  a  Target  and  Temporary  Cavity  due  to  a 
Penetrating  Object  using  a  Moving  Force  of  Variable  Intensity  Formulation 

When  a  projectile  of  small  size  e  g.,  a  bullet  penetrates  a  gelatin  or  similar  highly  viscous 
material,  the  trailing  path  of  the  projectile  may  partially  or  almost  fully  close  due  to  material 
settling.  The  response  at  a  given  point  Q  inside  the  target  is  very  complicated  due  to  the 
nonlinear  field  equations  and  temporal  boundary  conditions.  The  current  state  of  the  art  is 
represented  by  a  BRL  model^  that  assumes  that  such  moving  source  produces  spherically 
symmetric  waves  emanating  fi'om  its  instantaneous  location  so  that  the  total  impulse  any  given 
point  Q  (Figure  13)  is  given  by  the  integral 


R 

Zp  “X 

Impulse  at  Q  =  [  P{z)- — dz 

0  ^  (23) 

From  our  analysis  using  a  moving  point  source  formulation,  it  is  evident  that  the  BRL  model  is 
simply  a  far  field  point  source  solution  where  a  viscoelastic  damping  term  is  added  in  the  form  of 

_R 

the  exponential  product,  e  ,  in  the  integrand  in  Equation  23.  This  model  suffers  from  the 
following  deficiencies: 
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Figure  13.  Schematics  of  BRL  Calculation  of  Pressure  Field 


•  A  spherical  radiation  pattern  assumed, 

•  The  Force  intensity  is  not  an  explicit  function  of  projectile  velocity,  and 

•  Only  far  field  P-wave  response  is  considered. 

2.3.3  Spherical  Radiation  Pattern 

The  loading  on  the  target  due  to  a  penetrating  projectile  in  a  specified  direction  is  opposite  to  the 
retardation  force  that  the  projectile  experiences  during  its  motion.  For  an  elastic  medium,  the 
radiation  pattern  from  such  a  directive  force  can  be  calculated  from  the  corresponding  Green's 
function  in  the  medium  and  is  shown  in  Figure  14.  This  radiation  pattern  is  significantly  different 
from  the  spherical  wave  radiation  pattern  used  in  the  above  BRL  model,  especially  behind  the 
projectile  since  the  displacement  has  opposite  signs  in  these  two  models. 

2.3.4  Force  Intensity  Dependence  on  Instantaneous  Projectile  Velocity  and  Target 
Properties 

In  the  development  of  our  models  to  analyze  penetration  history,  we  concluded  that  the  force  of 


A 

^  (Field  Point) 


Y 

,.7 


Force 


A 


Figure  14.  Radiation  Pattern  due  a  Point  Force 
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retardation  depends  on  the  instantaneous  speed,  v,  of  the  projectile.  The  magnitude  of  the 
retardation  force  can  be  approximately  calculated  from  a  second  order  polynomial  in  v  of  the 
form  a\^+bv+c  where  a,  b  and  c  are  dependent  on  the  properties  of  the  target  and  the  projectile. 

2.3.5  Far  Field  P-wave  Response 

This  follows  from  our  moving  force  model  described  below.  With  reference  to  Figure  15,  let  the 
location  of  the  projectile  at  any  given  time  t  be  described  by  the  function:  x=l(t).  Behind  the 
projectile,  the  penetration  path  is  composed  of  partially  open  and  fully  closed  material.  The 
effect  of  the  open  penetration  path  introduces  additional  perturbation  in  the  medium  due  to  a 
vanishing  stress  field  at  these  locations.  We  do  not  know  at  this  point  how  much  influence  such  a 
perturbed  field  produces  at  a  given  field  point.  To  estimate  that,  we  first  assume  that  the 
penetration  path  closes  completely  after  its  passage.  If  we  ignore  any  effect  of  any  finite 
boundaries  the  response  of  the  medium  can  now  be  calculated  from  the  Green's  function  due  to  a 
time  varying  point  force  in  the  direction  of  the  projectile  located  at  its  instantaneous  position. 

The  physical  model  is  described  in  Figure  15. 

The  fundamental  assumptions  of  this  model  are  described  below: 

1.  A  moving  load  of  known  intensity  is  assumed  along  the  x-axis 

2.  The  force  intensity  is  a  known  function  of  the  velocity  v  in  the  form  av^+bv+c  where  the 
constants  a,  b,  c  depend  on  the  target  and  projectile  properties. 

3.  The  force  location  x=l(t)  is  a  known  function  of  time  t. 

Assumptions  in  2  and  3  follow  from  the  retardation  models  developed  in  previous  sections. 

2.3.6  Calculation  of  Transient  Response  due  to  an  Applied  Force 

The  displacement  at  a  point  x  at  time  t  due  to  a  point  force  at  the  origin  at  time  r  =  0  in  an  infinite 
medium  is  given  by  the  sum  of  the  following  three  terms: 


Figure  15.  Schematic  of MRC  Model  of  Moving  Point  Source 
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3  -S-  ^ 

1 .  Near  Field  Response:  Ui(x,  f)  =  f  tX{t  -  T)dr 

Anpr  * 

a 

2.  Far  Field  P-wave  Response:  ui{x,t)  =  -  — ) 

ATipa^r  « 


3.  Far  Field  S-wave  Response  Ui{x,t)  =  y  x{t  -  — )  (24) 

Anp^r  oc 

Closure  of  penetration  path  immediately  after  projectile  pass  is  assumed.  Effects  of  non-closure 
will  be  later  added  as  a  perturbation  to  the  above  solution. 

2.3.7  Calculation  of  the  Transient  Cavity  due  to  Projectile  Penetration 

The  mechanics  of  projectile  penetration  can  be  described  by  considering  the  target  as  a  three- 
dimensional  network  of  molecular  bonds  between  atoms.  When  a  projectile  penetrates  a  target,  it 
severs  molecular  bonds  between  particles  behind  the  projectile  while  the  particles  ahead  of  the 
projectile  experience  high  pressures.  The  radiation  pattern  of  these  pressure  waves  is  not 
spherical  in  shape,  but  rather  is  highly  directive  due  to  the  penetration-induced  directive  loading 
on  the  target.  The  shape  of  the  near-  field  and  far-field  radiation  pattern  for  this  directive  loading 
is  very  similar  those  given  by  the  Equations  24  even  though  these  equations  are  in  theory  valid  in 
an  infinite  medium  with  no  fracture  induced  free  surface  present. 

The  transient  cavity  is  formed  when  these  pressure  waves  try  to  displace  target  material  away 
from  the  fracture  region  while  no  molecular  resistance  to  such  motion  is  offered  from  the 
fractured  region  through  which  the  projectile  has  passed.  Thus,  at  any  instant  in  time,  the  target 
material  manifests  various  stress  distributions  along  with  new  stress-free  surfaces.  The  boundary 
of  the  transient  cavity  can  therefore  be  calculated  from  the  locus  of  zero-pressure  near  the 
projectile  path.  We  denote  this  surface  by  transient  cavity  boundary  1  or  tcbl. 

According  to  the  model  used  in  the  BRL  analysis  employing  a  spherical  radiation  pattern,  the 
boundary  of  the  transient  cavity  is  a  locus  of  critical  pressure  points  in  the  target  at  any  given 
point  in  time.  We  call  this  surface  as  the  transient  cavity  boundary  2  or  tcb2i.  The  difference  in 
the  transient  cavity  prediction  can  be  schematically  shown  in  Figure  16.  Since  the  pressure  is 
zero  on  tcbl,  by  definition,  and  is  again  zero  at  distances  far  from  the  fracture  region,  the 
variation  of  pressure  will  have  distribution  similar  to  that  shown  in  Figure  16. 

In  the  comparison  with  the  BRL  predicted  tcb2  and  MRC's  tcbl,  we  have  used  three  of  the  cases 
for  which  experimental  data  are  available.  These  are  solid  steel  spheres  penetrating  20% 
ordnance  gelatin  at  speeds  of  614,  783  and  984 ft/sec  respectively.  These  spheres  have  a  0.635- 
cm  diameter.  Since  we  have  assumed  in  our  analysis  that  the  path  of  penetration  closes  after  the 
passage  of  the  spheres,  the  lateral  displacement  is  zero  along  the  penetration  path  of  the 
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Figure  16.  Pressure  Distribution  during  Projectile  Penetration 

projectile.  For  the  calculation  of  tcbl,  we  need  to  determine  the  displacement  of  the  fractured 
particles  that  are  situated  along  the  trajectory  of  the  spherical  projectiles. 

We  know  that  the  lateral  displacement  is  zero  on  the  penetration  path  as  well  as  at  distances  far 
from  the  penetration  path  (radiation  condition).  Thus,  it  should  assume  a  local  maximum  at  some 
distance  from  the  penetration  path  (Figure  16).  We  can  then  calculate  the  maximum  lateral 
displacement  for  a  given  location  along  the  projectile  path  for  various  values  of  the  lateral 
coordinate.  As  expected,  this  maximum  occurs  very  close  to  the  penetration  path.  To  reduce  the 
singularity-induced  instability  (Equation  24  for  small  r)  in  our  calculation,  we  exclude  a  small 
region  near  the  point  source  for  which  the  contribution  at  a  given  field  point  is  calculated.  The 
pulse  width  for  each  point  source  is  taken  to  At  =  while  the  radius  of  the  excluded  zone  is 
taken  to  be  /?Af  where  p  is  the  shear  wave  speed  in  the  target.  The  calculated  results  are  then 
scaled  to  correspond  to  the  projectile  size  since  the  formulation  used  in  the  analytical  model 
corresponds  to  a  moving  point  force  in  the  target  medium.  The  retardation  coefficients  used  in 
the  MRC  developed  penetration  model,  corresponds  to  the  small  sphere  data  obtained  from 
correlation  with  MRC  generated  experimental  data  on  the  penetration  depths  of  these  spheres  in 
20%  ordnance  gelatin  targets.  These  coefficients  are,  in  cgs  units, 

a  =  0.227  (cm”^) 
b  =  -1126x10'*  (sec"*) 
c  =  1675x10^  (cm /sec^) 

The  problem  of  penetration  of  a  target  by  high  velocity  projectiles  can  be  divided  into  two 
problems.  First,  we  need  to  know  the  characteristics  of  the  retarding  force  acting  on  the  projectile 
and  then  we  also  need  to  calculate  the  target  response  due  to  projectile  penetration.  Experimental 
data  from  ballistic  testing  described  in  Section  3  has  been  implemented  for  different  formulations 
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of  ordnance  gelatin  targets.  These  data  show  the  distribution  of  projectile  penetration  depth  as  a 
function  of  striking  velocity.  An  analytical  inversion  process  has  converted  these  data  to 
determine  the  distribution  of  the  retarding  force  as  a  fonction  of  instantaneous  projectile  velocity 
inside  the  target.  The  distribution  of  the  retarding  force  yields  excellent  agreement  with  the 
experimentally  observed  projectile  velocity  history  inside  the  target. 

Let  us  now  look  at  the  problem  of  the  target  response  due  to  the  projectile  penetration. 

Physically,  this  involves  the  determination  of  the  temporary  transient  cavity  formed  behind  the 
projectile  as  well  as  the  final  long  term  cavity  known  as  permanent  cavity  along  the  penetration 
path  of  the  projectile  after  the  target  comes  to  rest.  The  second  problem  requires  the  knowledge 
of  the  viscoelastic  characteristics  of  the  target. 

To  solve  the  first  problem  by  direct  mathematical  modeling  of  the  actual  problem,  we  should 
consider  the  response  of  a  viscoelastic  medium  due  the  penetration  of  a  projectile  with  variable 
velocity  of  penetration.  The  temporal  character  of  the  velocity  of  penetration  is  assumed  to  be 
known  a  pr/or/.  Thus,  the  full  form  mathematical  model  is  a  three-dimensional  model  with  axial 
symmetry  along  the  penetration  path.  The  penetration  path  should  be  a  stress  free  surface  for  the 
model.  The  only  forcing  fiinction  in  the  target  medium  is  the  force  exerted  by  the  projectile  at  its 
instantaneous  location.  This  driving  force  is  also  assumed  to  be  known  since  this  force  is  exactly 
equal  and  opposite  to  the  force  of  retardation  experienced  by  the  projectile.  The  resulting 
problem  can  now  be  approximately  solved,  at  least  numerically,  by  using  the  full  three 
dimensional  form  of  the  representation  theorem  of  linear  elastic  theory  for  the  transient  part  of 
the  response  and  then  using  viscoelastic  theory  to  calculate  the  permanent  cavity  diameter. 

Some  first  and  second  order  estimates  of  the  transient  target  response  can  be  obtained  by 
simplifying  the  problem  by  eliminating  some  of  the  aspects  of  the  full  three-dimensional 
problem.  These  models  along  with  their  limitations  are  discussed  below: 

2.3.8  Three  dimensional  Target  Response  under  the  Assumption  of  Crack  closure 

If  we  assume  that  the  path  behind  the  projectile  closes  right  after  the  passage  of  the  projectile, 
then  we  can  get  the  target  response  by  superposing  the  solution  of  the  appropriate  Green's 
functions  for  point  sources  of  known  intensities  along  the  path  of  the  projectile.  These  Green's 
function  have  two  parts;  one  describes  the  response  at  a  point  inside  the  target  when  the  P-wave 
radiated  from  the  source  arrives  at  the  location  and  the  other  one  comes  from  the  S-wave  passage 
through  this  point.  Each  of  these  responses  has  their  near  and  far  field  representations.  These 
Green's  function  solutions  need  to  superpose  with  proper  time  shifts  to  incorporate  the  moving 
source  associated  with  the  projectile  motion.  The  intensities  also  are  time  dependent  since  they 
represent  the  intensity  of  the  force  that  is  a  function  of  the  instantaneous  velocity  of  the  projectile 
at  the  source  location.  The  shortcoming  of  such  a  simplified  model  response  is  that  the  lateral 
displacement  will  be  zero  on  the  projectile  path  (due  to  the  symmetry  of  the  problem  with  respect 
to  the  line  of  penetration)  as  well  as  zero  at  large  distances  from  the  sources  due  to  the  decay  of 
the  radiated  fields  with  distance.  For  the  original  problem,  the  displacement  should  be  maximum 
along  the  penetration  path.  Since  the  reduced  problem  has  a  local  maximum  at  a  point  away  from 
the  penetration  path,  we  can  extrapolate  the  transient  cavity  radius  from  such  simple  models. 
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2.3.9  Response  of  the  Target  due  to  a  Two-dimensional  Crack  Moving  with  Known, 
Variable  Velocity: 

If  only  the  early  response  of  the  target  is  of  importance,  we  may  assume  that  the  crack  formed  by 
the  projectile  penetration  is  a  two  dimensional  surface  because  the  radiation  from  the  other  stress 
free  surface  away  from  the  path  of  the  projectile  will  contribute  to  the  late  time  response  of  the 
target.  Here  again  we  assume  that  the  stress  ahead  of  the  crack  tip  is  known  i.e.,  the  stress 
intensity  factor  ahead  of  the  crack  tip  is  assumed  to  be  known  as  a  function  of  the  instantaneous 
velocity  of  the  projectile.  The  response  of  the  target  can  now  be  obtained  from  the  descritized 
version  of  the  representation  theorem  for  such  a  problem.  The  shortcoming  of  this  two- 
dimensional  crack  model  is  that  the  late  time  response  of  the  modified  function  has  a  long  tail 
while  the  original  problem  has  a  much  faster  decay  in  time. 

Once  the  analysis  of  the  results  from  these  two  simplified  models  are  complete,  we  will  be  able 
to  understand  the  physics  of  the  full  problem  more  clearly,  and  then  solve  the  fiill  three- 
dimensional  problem  described  above  if  such  solutions  are  warranted  for  better  correlation  of  the 
experimental  data  and  theoretical  analysis. 

2.3.10  Reduction  of  three-dimensional  inhomogeneous  material  distribution  to  an 
equivalent  non-linear  mass-spring  system 

In  the  determination  of  mechanical  response  of  parts  of  bodies  due  to  impact  or  other  similar 
loading,  it  is  currently  conceived  that  it  is  possible  to  replace  the  relevant  part  of  the  body  by  an 
equivalent  mass  spring  system.  There  are  obviously  some  questions  that  need  to  be  addressed 
before  such  idealization  is  useful.  For  a  given  body  part  the  main  issues  with  this  approach  are 
discretization  of  the  mass  particles  that  preserves  proper  mass  distribution  inside  the  body  as  well 
as  the  conservation  of  mass  for  the  given  body,  and  the  determination  of  the  equivalent  spring 
constants. 

Since  the  two  critical  parameters  in  the  mass-spring  model  of  a  given  body  part  are  the  particle 
mass  distribution  and  the  equivalent  spring  constants  connecting  these  mass  particles,  we  looked 
at  some  of  the  fundamental  issues  associated  with  their  description.  Discretizing  the  given  bodies 
mass  into  particle  masses  is  not  a  problem.  In  a  three-dimensional  scenario,  if  we  divide  a  given 
body  into  voxels,  any  mass  element  inside  a  cube  can  be  decomposed  into  eight  equal  particle 
masses  at  the  comers  of  the  elements.  Summing  all  the  particle  masses  at  any  cube  comer  from 
all  the  adjacent  elements,  we  get  the  equivalent  particle  mass  for  that  point. 

The  calculation  of  the  spring  constant  for  the  equivalent  mass-spring  system  is  not  that  simple. 
The  complexity  comes  from  the  ways  the  inhomogeneities  inside  a  cube  element  are  distributed. 
If  we  assume  that  all  components  of  a  body  part  have  known  properties,  each  component  can  be 
assigned  to  a  known  nonlinear  description  of  force-displacements  relations.  If  two  of  these 
elements  are  stacked  in  parallel,  then  we  can  calculate  the  equivalent  spring  constants  in  the 
following  way.  Let  us  assume  that  two  different  elements  are  situated  in  parallel  inside  a  cube 
element.  Then,  under  a  given  force  F  if  w  is  the  displacement  of  the  element,  then 
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2 

Fi  =  k^u  +  cju 

F2  =  k2U  +  C2U^  (25) 

F-ku  +  cu^  =  Fi+F2 

This  gives  the  equivalent  spring  constants  fc,c  in  terms  of  the  individual  spring  constants 
^,Cj  and  k2,C2  as  k  =  k\  +k\,c  =  q 


For  nonlinear  springs  in  series,  such  a  description  is  not  possible  simply  because  F  can  not  be 
written  in  terms  of  u=uj+u2  alone.  This  necessitates  the  imposition  of  critical  constraints  so  that 
a  given  cube  can  not  contain  two  different  material  types  in  series.  If  this  is  not  possible,  then  the 
inverse  relations  between  the  equivalent  spring  properties  and  individual  spring  properties  should 
be  used  even  though  they  do  not  exactly  satisfy  the  required  kinematic  relations.  If  this  is  done, 
then  an  error  estimate  should  be  obtained  to  give  some  credibility  to  the  solution. 

As  previously  mentioned,  the  penetration  problem  by  high  velocity  projectiles  can  be  divided 
into  two  sub-problems.  First,  we  need  to  know  the  characteristics  of  the  retarding  force  acting  on 
the  projectile  and  then  we  also  need  to  calculate  the  target  response  due  to  projectile  penetration. 
We  can  analytically  invert  ballistic  data  relative  to  projectile  penetration  depth  as  a  function  of 
striking  velocity  to  determine  the  retarding  force  acting  on  the  projectile  as  a  function  of  the 
instantaneous  projectile  velocity  inside  the  target. 

Let  us  now  look  at  the  problem  of  the  target  response  due  to  the  projectile  penetration. 

Physically,  this  involves  the  determination  of  the  temporary  transient  cavity  formed  behind  the 
projectile  as  well  as  the  final  long  term  cavity  known  as  permanent  cavity  along  the  penetration 
path  of  the  projectile  after  the  target  comes  to  rest.  The  second  problem  requires  the  knowledge 
of  the  viscoelastic  characteristics  of  the  target. 

To  solve  the  first  problem  by  direct  mathematical  modeling  of  the  actual  problem,  we  should 
consider  the  response  of  a  viscoelastic  medium  due  the  penetration  of  a  projectile  with  variable 
velocity  of  penetration.  The  temporal  character  of  the  velocity  of  penetration  is  assumed  to  be 
known  a  priori  since  we  have  already  shown  the  method  of  solution  for  such  velocity  histories. 
Thus,  the  full  form  mathematical  model  is  a  three-dimensional  model  with  axial  symmetry  along 
the  penetration  path.  Besides,  the  penetration  path  should  be  a  stress  free  surface  for  the  model. 
The  only  forcing  fiinction  in  the  target  medium  is  the  force  exerted  by  the  projectile  at  its 
instantaneous  location.  This  driving  force  is  also  assumed  to  be  known  since  this  force  is  exactly 
equal  and  opposite  to  the  force  of  retardation  experienced  by  the  projectile.  The  resulting 
problem  can  now  be  approximately  solved,  at  least  numerically,  by  using  the  full  three 
dimensional  form  of  the  representation  theorem  of  linear  elastic  theory  for  the  transient  part  of 
the  response  and  then  using  viscoelastic  theory  to  calculate  the  permanent  cavity  diameter. 

Some  first  and  second  order  estimates  of  the  transient  target  response  can  be  obtained  by 
simplifying  the  problem  by  eliminating  some  of  the  aspects  of  the  full  three-dimensional 
problem. 


39 


Mission  Research  Corporation 


DAMD17-94-C-4099 


2.3.11  Selected  Model  Results 


Figure  17  shows  predicted  scaled  displacement  versus  scaled  time  for  permanent  and  temporary 
cavities  using  a  quadratic  model  for  the  spring  stiffiiess.  The  governing  equation  of  motion  is 
given  below. 


m- 


<rx 

IF 


—x+—x^ 
m  m 


The  displacement  is  given  by  the  parameter,  x;  m,  t,  and  Ki  represent  the  lumped  mass,  time,  and 
stiffiiess  of  spring  /,  respectively.  We  assume: 

K.  X  —  ^  — 

X  =  =  and  t=  —  where  vo  is  initial  velocity;  then,  x  =  —  and  t  =  —  where  x  and  t  are  the 

Vo  Xs  K 

non-dimensional  scaled  parameters. 


Figure  1 7.  Effect  of  Energy  Absorption  on  Temporary  and  Permanent  Cavity 
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Figure  18.  Cavities  at  Various  Penetration  Depths  for  Non-Linear  Spring 
Model  with  Three  Particles  (K1  =  1.0  and K2=  0.35) 


Figure  19.  Cavities  at  Various  Penetration  Depths  for  Non-Linear  Spring 
Model  with  Five  Particles  (K1  =  1.0  andK2=  0.35) 
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Figure  18  and  Figure  19  shows  transient  and  permanent  cavity  displacement  versus  penetration 
depth  for  identical  non-linear  spring  models  but  where  the  lumped  mass  has  been  discretized  into 
3  and  5  particles.  Figure  20  shows  the  “affected  zone”  superimposed  on  Figure  19.  The  affected 
zone  is  established  by  identifying  those  areas  beyond  the  permanent  cavity  where  a  hysteresis  or 
inelastic  strain  accumulates.  A  hysteresis  effect  is  illustrated  in  Figure  21. 

A  typical  tissue  stress-strain  curve  is  shown  in  Figure  21.  The  projectile  impact  results  in  the 
material  being  loaded  up  to  point  A  in  Figure  21.  When  elastically  unloaded  to  point  B,  a  residual 
strain  or  hysteresis  effect  is  evident.  The  affected  zone  in  Figure  20  is  established  by  examining 
each  spring  in  the  lattice  (see  Figure  35  and  associated  discussion  regarding  non-linear  spring- 
mass  lattice,  all  Appendix  A)  to  determine  where  a  non-zero  residual  strain  exists.  If  the  initial 
part  of  the  stress-strain  curve  were  linear,  as  is  usually  the  case,  then  for  sufficiently  small- 
imposed  loads,  unloading  would  result  in  return  a  zero  residual  or  initial  strain  condition.  It  is 
only  when  the  load  extends  into  the  concave  portion  of  the  stress-strain  curve  and  then  elastically 
unloaded  that  a  hysteresis  effect  is  evident. 

Figure  22  shows,  as  a  function  of  penetration  depth,  the  size  of  the  temporary  cavity  for  a  344 
mg,  0. 171  inch  diameter  steel  sphere  striking  20%  ordnance  gelatin  at  1667  feet  per  second.  The 
plot  shows  predictions  using  different  discretization  schemes  where,  particularly  at  low 
velocities,  it  is  seen  that  convergence  is  very  fast. 


0  1  234567 


Scaled  Penetration  Depth 

Figure  20.  Affected  Zones  and  Cavities  as  a  Function  of  Penetration  Depth  for 
Non-Linear  Spring  Model  with  Five  Particles  (K1  =  1.0  andK2  =  0. 35) 
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STRESS 


Figure  21.  Schematic  of  Hysteresis  Effect 


Notice  in  Figure  22  that 
convergence  (i.e.,  the  offset 
between  the  various  curves)  is 
slower  at  high  penetration 
velocities  than  at  low  penetration 
velocities  (where  the  offset  between 
the  various  curves  decrease  very 
fast  with  increasing  numbers  of 
particles). 

Analytical  models  were  also 
developed  that  calculate  the  wound 
tract  geometry  for  motion  of  a 
bullet  or  other  similar  bodies  of 
high  aspect  ratio  in  a  bone-tissue 
environment.  We  assumed  that  the 
bullet  suffered  no  rotation 
until  it  hits  the  bone  so 
that  its  center  of  mass 
moves  in  a  straight  line.  It 
is  possible  to  relax  this 
condition  on  the  projectile 
kinematics  but  the 
resulting  problem 
becomes  three- 
dimensional.  Besides,  for 
most  problems  of  interest, 
this  condition  does  not 
impose  any  limitations  in 
applications  of  these 
models. 


The  problem  was 
formulated  as  follows.  We 
consider  a  target 
composed  of  tissue  matrix 
where  a  cylindrical  shape 
bone  is  embedded.  Here 
the  tissue  material  is 

replaced  with  a  20%  gelatin  or  any  other  suitable  medium  where  the  characteristic  of  the  force  of 
resistance  as  a  function  of  instantaneous  velocity  is  known.  We  should  recall  at  this  point  that  the 
coefficients  appearing  the  description  of  this  force  of  retardation  are  influenced  by  the  material 
properties  of  target  and  projectile  as  well  as  by  the  geometrical  shape  of  the  bullet  tip.  In  our 


Figure  22.  Transient  Cavities  Plotted  at  Different  Penetration 
Depths  for  Various  Numbers  of  Particles  (3,  5,  7,  9,  and  11)  in  Non- 
Linear  Spring  Model 
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earlier  analysis,  we  developed  a  sufficiently  accurate  model  that  describes  the  force  of 
retardation  as  a  cubic  polynomial  in  the  instantaneous  velocity. 

For  bullets  impacting  bones,  depending  on  the  entry  velocity  of  the  bullet  prior  to  hitting 
the  bone,  various  modes  of  motion  can  be  classified  as  follows. 

(1)  No  penetration  of  hone.  In  this  case,  the  projectile  velocity  is  not  enough  to  penetrate  so  that 
the  bullet  either  turns  or  rebounds.  This  threshold  velocity  is  approximately  200  ft/sec  for 
osteoporotic  bones  and  700  ft/sec  for  normal  bones. 

(2)  Drill  hole  mode:  This  mode  occurs  for  certain  entry  velocity  of  the  bullet.  Bulle'  penetrates 
like  a  drill  hole  with  clean  entry  and  exit  holes. 

(3)  Brittle  fracture  mode:  For  a  given  bullet  and  certain  range  of  entry  velocity,  the  bone  suffers 
radial  fractures  that  are  sometimes  so  bad  that  the  whole  bone  disintegrates,  the  main 
mechanism  for  these  fractures  is  the  propagation  of  highly  compressive  and  reflected  tensile 
waves  through  the  target. 

We  have  successfully  developed  models  which  use.  an  already  existing  two-dimensional  code 
for  determining  the  wound  track  in  a  homogeneous  medium.  This  code  is  first  used  to  determine 
the  motion  of  the  bullet  just  prior  to  hitting  the  bone.  If  the  bullet  velocity  at  this  instant  is  below 
threshold  velocity,  the  bullet  will  either  rebound  or  turn.  This  produces  some  discontinuities  in 
the  bullet  velocity  (both  translational  and  angular).  These  discontinuities  have  been  calculated 
using  impulse  equations  and  then  the  data  on  linear  and  angular  velocity  of  the  bullet  is  fed  to  the 
above  two-dimensional  code  as  initial  values  for  the  calculation  of  the  subsequent  wound  track 
geometry. 

These  cases  of  the  drill-hole  mode  can  also  be  handled  the  same  way  using  the  same  code  as 
above.  In  this  case,  however,  the  bullet  continues  to  move  in  a  straight  line  after  entering  the 
bone  with  different  retardation  coefficients  than  that  of  a  gelatin  target. 

A  possibility  of  a  critical  shear  force  model  has  also  been  tested  for  applications  in  the  cases  of 
drill-hole  mode  fractures  for  low  velocities.  Experimental  data  shows  that  for  the  osteoporotic 
bones,  the  average  force  of  resistance  is  11.68  lb.  for  0.250-in  diameter  spherical  projectiles 
while  it  is  12  lb.  for  0.406-in  diameter  projectiles  moving  at  200-400  fl/sec.  This  data  justifies  a 
critical  force  criterion  for  low  velocity  penetration  of  bullets  in  bones. 
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2.4  Kinematics  of  Bullet  Penetration  in  a  gelatin  with  embedded  bones 


In  general,  when  a  projectile  enters  a  target,  it  experiences  resistance  to  penetration.  The 
resistance  comes  from  the  contact  area  of  the  projectile  with  the  target  medium.  The  force  of 
resistance  is  distributed  over  this  contact  area.  If  we  know  how  these  forces  are  distributed  along 
with  their  directions  and  magnitudes  we  can  exactly  determine  the  path  of  the  projectile  inside 
the  target.  From  the  analysis  that  has  been  done,  it  has  been  found  that  the  magnitude  of  the  force 
of  resistance  at  any  point  depends  on  the  instantaneous  velocity  of  the  projectile  at  the  point  of 
contact  while  the  direction  is  usually  assumed  to  be  in  the  direction  of  the  normal  component  of 
the  velocity  with  respect  to  the  projectile  surface.  This  is,  however,  not  true  in  cases  where  high 
viscosity  is  present  in  the  target  medium  and  shear  forces  may  contribute  to  the  force  of 
resistance.  The  generation  of  shear  forces  does  not  complicate  the  problem  if  the  nature  and 
magnitude  of  these  forces  are  known  in  terms  of  the  instantaneous  velocity,  v.  Thus,  when  the 
projectile  is  completely  immersed  and  moving  in  the  target,  the  kinematics  of  the  projectile  can 
be  calculated  by  writing  down  the  equations  of  motion  of  the  projectile  for  the  center  of  mass 
and  about  the  center  of  mass  {CM).  For  the  case  of  a  projectile  moving  in  an  ordnance  gelatin 
target,  these  equations  are  given  by  (Figure  12,  page  28) 

i 

mxG  =  -J/?i?(|i^pn|)sgn(i;p„)cosAir  -myi?(|ue„|)sgn(i;e„)sin^ 

~1  (26) 

^ilG  =  +j  pJ^i\vpn\)^Q^(Ppn)^'^0dr -mfR{\vgn\)sgn{Ven)cos0 
-l 

The  equation  of  motion  about  the  CM  can  be  obtained,  as  before,  by  taking  moments  of  the 
equations  of  motion  of  the  ensemble  of  particles  constituting  the  projectile  about  the  CM.  This 
yields 


i 

lGO=-\pR(\Vpn\^Qr\{Vpn)rdr  (27) 

-I 

When  a  bullet  or  a  projectile  is  entering  a  surface  corresponding  to  a  material  discontinuity  e.g., 
in  the  early  phase  of  gelatin  penetration  from  air  or  early  phase  of  bone  penetration  from  the 
gelatin/tissue  side,  the  projectile  kinematics  is  more  complicated  that  of  the  completely 
immersed  case.  In  this  case,  the  part  of  the  projectile  in  contact  with  the  surface  of  material 
discontinuity  suffers  a  different  distribution  of  resistance  forces  than  the  rest  of  the  projectile. 
Thus,  Equations  26  and  27  should  contain  two  integrals;  one  for  the  contact  area  of  the  material 
discontinuity  and  other  for  the  area  that  has  not  yet  penetrated  the  medium.  This  introduces 
discontinuities  in  the  velocity  distribution  for  both  the  linear  and  angular  components  of  the 
projectile.  It  is  necessary  to  analyze  these  issues  in  order  to  describe  projectile  penetration 
through  a  gelatin/soft  tissue  embedded  bone.  Below  we  give  some  analysis  of  this  problem, 
which  may  yield  some  information  about  the  physical  nature  of  the  force  of  resistance  coming 
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from  the  surface  of  discontinuities.  We  also  mention  some  of  the  fracture  modes  observed  during 
bone  penetration. 

Analysis  of  Various  Modes  of  Bone  Fracture 

It  has  not  yet  been  conclusively  established  what  parameters  of  the  dynamic  process  of  projectile 
penetration  are  responsible  for  bones  to  preferentially  fracture  in  a  specific  mode.  Two  important 
parameters  of  penetration  dynamics  are  projectile  size  and  entry  velocity.  The  kinetic  energy  of 
the  projectile  is  directly  proportional  to  both  of  these  quantities.  From  ballistic  impact 
experiments  reported  in  the  literature  involving  embedded  bones  in  ordnance  gelatin  targets  and 
spherical  projectiles  of  various  sizes;  it  is  seen  that  increasing  the  size  of  the  spherical  projectile 
is  effective  in  penetrating  the  bone  with  little  or  no  radial  fractures  on  the  anterior  side  while  a 
decrease  in  the  size  of  the  spherical  projectile  is  less  effective  in  penetration  but  effective  in 
creating  brittle  mode  radial  fractures. 

For  example,  an  osteoporotic  femur  subject  to  impact  by  a  0.406-in.  sphere  at  600  fl/sec 
penetrates  the  bone  with  little  or  no  radial  fractures  while  a  0.250-in.  sphere  at  602  ft/sec  did  not 
penetrate  but  produced  radial  fractures  emanating  from  the  entrance  hole.  If  we  increase  the 
speed  of  the  smaller  size  sphere  to  908  ft/sec,  the  projectile  also  stops  but  the  bone  is  broken  into 
two.^  This  case  has  about  41%  more  kinetic  energy  than  that  of  the  larger  size  sphere  but  no 
penetration  occurs.  This  shows  that  the  distribution  of  energy  for  smaller  size  spheres  is  more 
focused  towards  radiation  of  energy  than  towards  penetration.  Hence  the  increase  in  kinetic 
energy  alone  does  not  necessarily  guarantee  penetration.  We  need  to  know  how  projectile  size 
and  its  velocity  affect  the  distribution  of  energy. 

To  calculate  the  radial  fracture  from  compressive  waves,  MRC  has  developed  a  singular 
displacement  model  where,  to  the  lowest  order  in  singularity,  the  radial  displacement  decays  like 
J/r  where  r  is  the  distance  from  the  center  of  the  projectile  footprint.  For  smaller  footprints,  it 
radiates  more  intense  compressive  waves.  This  model  gives  the  same  qualitative  result  as 
observed  in  the  above  experiments.  For  fractures  in  coated  and  uncoated  polycarbonate 
substrates,  these  models  yield  quantitative  estimates  of  radial  fracture  that  is  in  extremely  good 
agreement  with  experiments.'*  However,  this  mechanism  of  energy  distribution  has  not  yet  been 
completely  understood.  More  systematic  impact  experiments  are  needed  where  one  parameter  of 
impact  should  be  changed  while  others  remain  unchanged  so  that  the  effect  of  a  particular 
parameter  on  the  penetration  mechanism  can  be  understood. 

2.4.1  Characteristics  of  Bone  Impact  and  Penetration 

When  a  projectile  penetrates  a  target  modeled  by  an  embedded  bone  in  a  gelatin/tissue,  host 
matrix,  and  its  penetration  characteristics  after  it  impacts  the  embedded  bone  depends  mostly  on 
the  following  quantities; 

a.  Incident  velocity  of  the  projectile  immediately  before  it  impacts  the  bone 

b.  Projectile  material,  shape  and  size 

c.  Bone  condition  (Normal,  Osteoporotic,  Mildly  Osteroporotic  etc.) 
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Table  2.  Energy  and  Velocity  Data  on  Spherical  Steel  Projectiles  Impacting  Human  Femurs 
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Table  2  shows  data  from  ballistic  tests  performed  on  embalmed  human  femurs  with  spherical 
projectiles.  These  data  were  garnered  from  the  literature  and  represent  unloaded,  normal  femurs. 
The  literature  sources  are  between  14  and  50  years  old  and  therefore  limited  in  the  types  of  data 
that  could  be  extracted  from  these  experiments. 

From  the  data  that  has  already  been  examined  however,  it  appears  that  a  specific  impulse  of  0.4 
Ibf-sec/in^  represents  the  threshold  required  to  cause  a  spiral  fracture  in  the  femur.  This 
corresponds  to  a  dynamic  pressure  of  2.16  ksi.  This  is  similar  to  what  was  recorded  on  ballistic 
experiments  discussed  earlier.  A  spiral  fracture  is  the  fracture  mode  typically  observed  when  the 
projectile  does  not  directly  strike  the  bone  but  the  projectile  passes  close  to  the  bone.  A  point 
load  of  approximately  287  bounds  appears  sufficient  to  cause  a  complete  tensile  failure  of  the 
femur. 


0  200  400  600  800  1000  1200  1400  1600  1800  2000 


Figure  23  through  Figure 
25  reduce  the  data  from 
this  table  in  different 
ways.  Figure  23  shows 
striking  velocity  versus 
exit  velocity  for  two 
different  diameter  spheres. 
Figure  25  shows  energy 
loss  versus  projectile 
average  velocity  during 
penetration  of  the  bone.  In 
these  two  plots  the  results 
show  considerable 
divergence  for  the  two 
different  diameter  spheres. 


Exit  Velocity  (ft/sec) 

Figure  23.  Entry  and  Exit  Velocity  of  Steel  Projectiles  Striking 
Embalmed  Human  Femurs 


Figure  25  shows  energy 
loss  normalized  to  the 
projected  area  of  the 
sphere  versus  average 
projectile  velocity  during 

penetration.  In  this  case  the  results  collapse  onto  each  other.  This  is  consistent  with  the 
mechanics  of  fracture  where  the  “energy  density”  (i.e.,  energy  loss  per  unit  area)  is  a  parameter 
used  to  describe  fracture  initiation  and  propagation.  More  generally,  the  criterion  should  be 
energy  density  per  unit  time.  For  spheres  however  the  contact  area  per  unit  time  is  probably 
similar.  We  would  expect  to  see  profound  differences  however  in  the  time  history  of  the  contact 
area  for  projectiles  with  different  leading  edge  geometries. 


We  prepared  a  DARPA  proposal  with  Dr.  Paul  Dougherty,  Maj.,  USAMC  (orthopedic  surgeon  at 
Fort  Knox)  to  conduct  ballistic  testing  of  human  femurs  in  a  gelatin  matrix  while  subject  to  a 
preload.  This  would  provide  an  opportunity  to  validate  analytical  predictions  and  develop  more 
relevant  data  that  elucidate  the  issues  above. 
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Figure  24,  Energy  Loss  and  Average  Velocity  (Scaled  Display) 
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Figure  25.  Energy  Loss  per  unit  Projected  Area  and  Average  Velocity 
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Figure  26.  Nomenclature  c^sociated  with  pre-  and  post  impact 

kinematics 


If  F  is  of  the  order  O(v^),  then  ^  =  n -9  and  we  have  a  complete  rebound  of  the  bullet.  Thus,  for 
penetration  F  is  at  least  0(v)  that  gives  the  new  angle  (p  as 


— sin6>At 

^  =  - 

1-— cos6>At 
m 

F  =  av 

If  F=0(v),  then  (p=0  so  that  bullet  moves  with  unchanged  direction. 

Since  the  penetration  of  projectiles  in  gelatin  are  well  understood  and  the  associated  algorithms 
at  MRC  give  reasonably  good  results,  we  may  solve  the  penetration  problem  up  to  the  time  it 
makes  contact  with  the  bone.  After  this,  we  need  to  establish  the  nature  of  the  resistance  force  as 
a  fimction  of  velocity  in  order  to  predict  the  projectile  kinematics  resulting  from  bone 
interaction.  These  distributions  can  be  obtained  from  the  geometry  of  the  wound  tract  observed 
ballistic  experiments  involving  targets  with  embedded  bones  and  inversion  of  penetration  depths 
under  various  impact  scenarios. 

2.4.4  “Drill  Hole”  Mode 

This  mode  occurs  when  v^cr  ^2cr  •  mode,  the  bullet  leaves  behind  a  wound  tract 

similar  to  a  drill  hole  with  clean  entry  and  exit  holes. 

2.4.5  Brittle  Fracture  Mode:  Bone  Fragmentation 

When  Vffi  >  V2cr ,  the  bone  manifests  radial  fractures  that  are  so  severe  that  the  whole  bone 
disintegrates.  The  main  mechanism  for  these  fractures  is  the  propagation  of  highly  compressive 
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waves  through  the  target.  Analytical  models  have  been  developed  at  MRC  to  interpret  the 
existence  of  and  to  quantify  such  radial  fractures.  In  the  BONGEL  Code  (see  Section  2.5),  bone 
fragments  and  their  final  resting  locations  are  modeled  by  a  uniform  random  distribution. 

2.4.6  Uniform  Random  Distribution 

Since  the  injury  data  on  bone  fragmentation  is  very  limited,^  it  is  not  possible  to  determine  the 
probabilistic  distribution  of  fragment  size,  number  of  fragments  and  their  final  resting  location  in 
the  body  after  bone  disintegration  takes  place  due  to  high  velocity  of  impact.  For  the  purpose  of 
modeling,  we  have  assumed  that  these  distributions  are  uniformly  random  but  some  physical 
characteristics  have  been  maintained  to  match  experimental  observations. 

The  maximum  number  of  bone  fragments  used  in  BONGEL  Code  is  10.  This  result  is  in 
conformity  with  observations  of  reported  bullet  injuries  that  resulted  in  bone  fragmentation.  The 
energy  partition  and  distribution  of  bone  fragments  is  adjusted  so  that  larger  projectiles  have  a 
greater  proportion  of  kinetic  energy  devoted  to  the  penetration  (as  opposed  to  the  brittle  fracture 
process. 

2.4.7  Final  Location  of  Bone  Fragments 

The  penetration  depth  of  projectile  and  bone  fragments  is  calculated  in  BONGEL  using  the 
previously  described  penetration  algorithm  with  a  known  initial  velocity.  The  circumferential 
and  radial  distribution  along  the  mean  fragment  path  is  calculated  from  a  random  number 
generator  using  the  assumption  of  a  uniform  random  distribution. 


2.5  Development  of  The  BONGEL  Code 


MRC  has  developed  the  BONGEL  code  to  analyze  various  aspects  of  the  kinetics  and  kinematics 
of  projectile  penetration  in  a  tissue  target  with  embedded  bone.  To  retain  simplicity  during  the 
initial  phase  of  the  development,  the  bone  is  modeled  as  a  right  circular  cylinder.  In  subsequent 
development,  described  in  Section  4,  the  bone  is  allowed  to  take  any  other  shape  occupying  a 
given  volume  of  the  target.  In  Section  4,  the  bone  geometry  is  read  from  the  virtual  anatomy 
from  the  MGI  Lower  Extremity  Trauma  Simulator.  Three  types  of  projectiles  can  be  used  in  the 
code;  these  are  (1)  bullets,  (2)  spherical  projectiles,  and  (3)  fragments. 

In  the  code,  these  projectiles  are  identified  by  the  code  parameter  PTYPE  that  takes  the  value  1, 
2,  or  3  for  bullet,  sphere  or  fragment  type  projectile,  respectively.  Since  the  input  requirements 
for  different  types  of  projectiles  may  differ,  for  the  ease  of  the  user,  projectile  properties  can  be 
entered  or  modified  independently.  These  files  are  BULLET.in,  SPHERE.in  and 
FRAGMENT.in,  respectively.  When  the  projectile  type  PTYPE  is  prescribed,  the  BONGEL 
code  generates  its  own  input  file,  BONGEL.in,  by  modifying  the  user  specified  projectile  input 
files  as  necessary.  These  files  also  contain  impact-scenario  information,  e.g.,  velocity,  location 
and  orientation  of  the  projectile  prior  to  impact. 
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Other  than  writing  various  graphics  and  other  information  files,  this  code  is  not  a  stand-alone 
code.  The  development  has  been  done  in  modular  forms  so  that  each  module  can  be  modified 
independent  of  other  modules.  There  are  two  screen  utility  codes  that  need  to  be  present  in  the 
directory  or  file  path.  These  files  are  necessary  only  when  any  modifications  of  this  code  are 
made  and  new  compilation  is  performed.  These  two  screen  utility  modules  are  CLS  (clear 
screen)  and  TLOCATE  (cursor  positioning).  These  two  modules  are  only  provided  in  object- 
module  forms  so  that  only  compilations  can  be  done.  No  source  codes  for  these  two  modules  are 
either  provided  or  deemed  necessary. 

Besides  these  two  screen  utility  codes,  other  stand-alone  modules  are  also  used  by  the  BONGEL 
code  as  necessary.  When  this  is  done,  the  BONGEL  code  creates  the  proper  input  files  prior  to 
running  these  codes.  There  are  two  such  stand-alone  codes  that  BONGEL  uses;  these  are 
BULPENT2.for  and  CAVLOC2.for.  The  BULPENT2  code  is  used  in  the  analysis  of  bullet 
penetration  through  tissue  or  bone  medium  while  the  CAVLOC2  code  is  used  to  generate 
temporary  and  permanent  cavity  locations  as  the  projectile  passes  through  the  target.  MRC  has 
also  developed  a  stand-alone  code  named  MM4  that  uses  a  one-dimensional,  nonlinear  spring 
model  to  calculate  temporary  and  permanent  cavity  location  for  a  chain  of  particles  where  only 
one  end  of  the  chain  is  given  an  initial  velocity.  Since  the  repeated  use  of  the  code  in  various 
directions  at  all  locations  of  the  projectile  is  necessary  to  generate  the  three-dimensional 
temporary  and  permanent  cavity  locations,  and  hence,  is  very  much  time  consuming,  BONGEL 
uses  analytical  fits  obtained  fi’om  MM4  code  to  predict  temporary  and  permanent  cavities. 
BONGEL  prepares  an  input  file  for  CAVLOC2  prior  to  calling  this  module  for  the  calculation  of 
three-dimensional  temporary  and  permanent  cavity  surfaces. 

As  mentioned  earlier,  code  developments  are  made  in  modular  forms  so  that  one  module  may 
execute  another  module  if  necessary.  In  order  to  implement  such  uses  of  other  modules  fi'om  any 
other  module,  BONGEL  uses  a  SYSTEM  command  that  is  provided  in  the  MICROSOFT 
FORTRAN  Version  5  and  above.  Thus,  other  FORTRAN  compilers  lacking  the  SYSTEM 
command  are  not  capable  of  compiling  this  version  of  the  code.  Some  typical  examples  of  such 
uses  are  shown  in  the  following  code  lines 

l=SYSTEM('COPY  SXYD.PLT  SXYD2.PLr) 
l=SYSTEM('COPY  SYYD.PLT  SYYD2.PLT') 
l=SYSTEM('COPY  BULLIMP2.IN  BULPENT2.IN’) 
l=SYSTEM('BULPENT2') 
l=SYSTEM(‘CAVLOC2’) 

In  the  first  three  examples,  the  code  uses  the  DOS  copy  command  to  copy  some  necessary  files 
while  in  the  last  two  examples,  the  code  executes  the  modules  BULPENT2  and  CAVLOC2 
using  DOS  commands. 

A  completely  integrated  version  of  this  code  can  be  generated  easily  if  other  compilers  lacking 
the  SYSTEM  command  are  necessary. 
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The  following  files  are  necessary  to  run  BONGEL: 

BULLET.IN 

SPHERE.IN 

FRAGMENT.IN 

BONGELEXE 

BULPENT2.EXE 

CAVLOC2.EXE 

Other  inputs  files  are  automatically  generated  by  BONGEL. 

The  following  files  are  necessary  to  compile  BONGEL: 

BONGELFOR 

CLS.OBJ 

TLOCATE.OBJ 

along  with  a  FORTRAN  compiler  similar  to  Microsoft  FORTRAN  Version  5  and  above. 
A  basic  flowchart  of  BONGEL  code  is  shown  on  the  next  page  in  Figure  27. 
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Figure  27.  Flowchart  ofMRC BONEGEL  code 
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2.5. 1  Enhanced  Tissue  Damage  due  to  Ejected  Bone  Fragments 

When  bone  fragmentation  occurs  due  to  high  impact  velocity,  these  fragments  move  inside  the 
tissue  medium,  and  create  additional  damage.  This  damage  is  usually  referred  to  as  enhanced 
tissue  damage.  This  part  of  the  damage  calculation  is  not  incorporated  in  the  current  version  of 
the  code  due  to  a  lack  of  program  resources.  As  previously  mentioned  there  are  no  fundamental 
techmcal  obstacles  to  incorporating  this  enhanced  damage  and  this  phenomenon  could  be 
incorporated  at  a  future  date. 

2.5.2  Modifications  Made  to  Incorporate  Arbitrary  Bone  Shape  and  Location 

The  BONEGEL  code  includes  an  interface  to  the  virtual  anatomy  included  in  the  AfG/  Lower 
Extremity  Battlefield  Trauma  Simulator.  This  interface  is  discussed  further  in  Section  4.  The 
virtual  anatomy  is  available  in  a  digitized  form  so  that  when  a  coordinate  is  specified  with 
respect  to  a  user  defined  rectangular  coordinate  system,  the  associated  part  can  be  identified.  We 
therefore  modified  the  BONGEL  code  so  that  for  a  given  impact  scenario,  entry  and  exit  points 
on  the  bone  can  be  identified  in  cases  where  the  projectile  impacts  the  bone  and  a  bone 
penetration  may  occur.  This  has  been  done  as  follows. 

When  the  impact  scenario  is  specified  by  the  INFLICT  module  discussed  in  Section  4,  the 
direction  cosines  (^,m,n)  and  the  coordinates  (a,b,c)  ofthe  point  of  entry,  4  in  the  body  are 
known.  Thus  the  coordinates  {x,  y,  z)  of  any  other  point,  P,  on  the  line  of  penetration  of  the 
projectile  can  be  written  in  terms  of  the  distance  AP  =  r  as 

x  =  a  +  £r 
y  =  b  +  mr 
z  =  c  +  rtr 

By  tracking  the  coordinates  (jc,  z)  as  the  projectile  moves  (gradual  increase  in  r  )  along  its  line 
of  penetration,  we  can  determine  m  advance  if  bone  impact  occurs,  and  if  so,  its  entry  and  exit 
points  in  cases  where  bone  penetrations  occur.  This  algorithm  can  replace  the  relevant  part  of  the 
BONGEL  code  when  the  integration  of  BONGEL  with  the  digitized  database  is  made. 

2.5.3  Axis  Systems  used  in  BONGEL  Code 

Frame  F;  (Oxyz)  Global  or  User  Defined  Coordinate  System 
(Qx’y’z’)  Local  coordinate  system  at  the  bone  impact  point  Q 
{£,  m,  n)  direction  cosines  of  the  original  bullet  shot  line  with  respect  to  the  frame  F. 

{a,  p,  Y )  direction  cosine  of  the  impulsive  force  at  Q  with  respect  to  F.  This  direction  is 

along  the  outward  normal  direction  to  the  bone  surface  at  the  impact  point  Q.  Plane  n  is  the 
plane  containing  the  lines  and  (a.p.r) . 
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(^1,  mj,  Mj)  direction  cosines  of  the  new  direction  of  the  projectile  center  of  mass,  G,  velocity 
right  after  bone  impact  in  the  case  where  rebounding  occurs  after  bone  impact.  This  forms  the 
new  x’  axis  at  Q.  This  is  also  referred  to  as  LP,MP,NP  in  the  BONGELT  code. 

Cross  product  of  the  lines  {£,m,n)  and  {a,p,Y)  forms  the  new  z’  axis  at  Q. 

The  direction  cosines  of  the  Qz’  axis  are  denoted  as  (^2>^2'^2) 

The  cross  product  of  z’-  and  x’-axes  at  Q  forms  the  new  y’  axis  at  Q. 

The  direction  cosines  of  the  Qy’  axis  are  denoted  as  (^3,  m3,  n3) 

2.5.4  File  Structures  in  BONGEL  Code 


The  generation  of  internal  files  in  BONGELT  and  Other  Called  Programs  (RK2S)  and  Codes 

(BULPENT2) 


RK2S  Subroutine  (For  One  dimensional  motion  x  =  x(t) ,  x  -space,  t  -time) 


File 

Unit 

Description 

SXY.PLT 

9 

X  vs.  f 

SXYD.PLT 

10 

dx  . 

—  vs.  t 
dt 

SYYD.PLT 

11 

dx 

xvs.  — 
dt 

STPHS.COD 

15 

Code  Control  Parameter 

BULLPENT2  Code  (Run  by  BONGELT2  Code,  for  Two-dimensional  Motion 

X  =  x{t),  y  -  y{t) ;  Subscript  G  indicates  the  center  of  mass  of  the  bullet,  subscript  1 
front  tip  while  subscript  2  indicates  back  tip  of  the  bullet) 

0  indicates  bullet  orientation  with  the  x-axis. 

File 

Unit 

Description 

BXG.PLT 

2 

Xq  vs.  t 

BXGD.PLT 

3 

^  vs.t 
dt 

BYG.PLT 

4 

Vg  vs.  t 

BYGD.PLT 

5 

— ^  vs.  t 

dt 

BTH.PLT 

6 

0{6eg)  vs.  t 

BTHD.PLT 

7 

d0 

—  vs.  f 

dt 

BULLPENT2.DBG 

8 

Debug  File 

BXY1.PLT 

9 

(x,i/)  for  tip  1 

BXY2.PLT 

10 

{x,y)  for  tip  2 

BXYG.PLT 

11 

{x,y)  for  the  center  of  mass  G 

BYAW.PLT 

12 

Yaw  angle  vs.  time  t 

BULLPEN.ANM 

14 

Graphic  Animation  File 

indicates  the 
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STPHS.COD  16  Code  control  parameter 

BXY1 D. PLT  1 8  Velocity  components  for  the  front  tip 

along  the  x-  and  y-axis. 

Note:  These  file  names  have  been  modified  by  adding  1, 2  or  3  after  their  names  to  indicate  the 
data  obtained  during  phase  1, 2  or  3  of  the  projectile  penetration  where  the  different  phases  of 
motion  are  defined  as  follows: 

Phase  1:  Projectile  motion  prior  to  bone  impact 

Phase  2:  Projectile  rebound  mode  or  penetration-through-bone  mode 

Phase  3:  Projectile  reentry  in  tissue  medium  after  bone  penetration  is  completed. 

Significant  efforts  have  been  made  in  developing  analytical  models  capable  of  calculating  wound 
tract  geometry  for  the  motion  of  a  bullet  or  other  similar  bodies  of  high  aspect  ratio  in  a  bone- 
tissue  environment.  We  assumed  that  the  bullet  suffered  no  rotation  until  it  hits  the  bone  so  that 
its  center  of  mass  moves  in  a  straight  line.  Although  this  is  not  the  most  general  case,  for 
extremities  that  tend  to  have  short  wound  paths,  more  often  than  not,  the  projectile  follows  a 
rectilinear  path  unless  an  interaction  with  bone  tissue  occurs.  This  tends  to  not  be  true  in  other 
regions  of  the  body  where  the  wound  tracts  are  substantially  longer.  We  have  modeled  the  three 
dimensional  curvilinear  path  the  projectile  follows  after  bone  interaction  however  we  elected  not 
to  include  non-rectilinear  paths  in  soft  tissue  prior  to  bone  interaction.  This  scope  was  defined  in 
the  proposal  and  does  not  represent  any  fundamental  technical  obstacles  but  merely  a  decision 
relative  to  managing  resources. 

It  is  possible  to  relax  this  condition  on  the  projectile  kinetics  but  the  resulting  problem  becomes 
three-dimensional.  Besides,  for  most  problems  of  interest,  this  condition  does  not  impose  any 
limitations  in  applications  of  these  models.  The  problem  we  are  currently  working  on  can  be 
formulated  as  follows.  We  consider  a  target  composed  of  tissue  matrix  where  a  cylindrical 
shaped  bone  is  embedded.  Here  the  tissue  material  is  replaced  with  a  20%  gelatin  or  any  other 
suitable  medium  where  the  characteristic  of  the  force  of  resistance  as  a  function  of  instantaneous 
velocity  is  known.  We  should  recall  at  this  point  that  the  coefficients  appearing  in  the  description 
of  this  force  of  retardation  are  influenced  by  the  material  properties  of  target  and  projectile  as 
well  as  by  the  geometrical  shape  of  the  bullet  tip.  In  our  earlier  analysis,  we  have  developed  a 
sufficiently  accurate  model  that  describes  the  force  of  retardation  as  a  cubic  polynomial  in  the 
instantaneous  projectile  velocity,  v .  For  bullets  impacting  bones,  depending  on  the  entry 
velocity,  ,  of  the  bullet  prior  to  hitting  bone,  three  interaction  modes  can  be  classified  as 

follows. 

(1)  No  penetration  of  bone.  In  this  case,  Vf.  is  not  sufficient  to  effect  penetration  so  that  the 
projectile  either  turns  or  rebounds.  This  threshold  velocity,  ,  is  approximately  200  ft/sec 
for  osteoporotic  bones  and  700  ft/sec  for  normal  bones. 

(2)  Drill  hole  mode.  This  mode  occurs  for  certain  entry  velocities  where  the  projectile  penetrates 
like  a  “drill  hole”  with  clean  entry  and  exit  holes. 
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(3)  Brittle  fracture  mode.  For  a  given  bullet  and  certain  range  of  entry  velocity,  the  bone  suffers 
radial  fractures  that  are  sometimes  so  severe  that  the  whole  bone  disintegrates.  The  main 
mechanism  for  these  fractures  is  propagation  of  highly  compressive  and  reflected  tensile 
waves  through  the  target. 

We  have  successfully  developed  models  that  used  an  existing  two-dimensional  code  for 
determining  the  wound  track  in  homogeneous  medium.^  This  code  is  first  used  to  determine  the 
motion  of  the  bullet  just  prior  to  hitting  the  bone.  If  the  bullet  velocity  at  this  instant  is  below, 

Vrf,  the  bullet  will  either  rebound  or  turn.  This  produces  some  discontinuities  in  the  bullet 
velocity  (both  transnational  and  angular).  These  discontinuities  have  been  calculated  using 
impulse  equations  and  then  the  data  on  linear  and  angular  velocity  of  the  bullet  is  fed  to  the 
above  two-dimensional  code  as  initial  values  for  the  calculation  of  the  subsequent  wound  track 
geometry. 

The  drill-hole  mode  can  also  be  handled  the  same  way  using  the  same  code  as  above.  In  this 
case,  however,  the  bullet  continues  to  move  in  a  straight-line  path  after  entering  the  bone  with 
different  retardation  coefficients  than  that  of  a  gelatin  target. 

A  possibility  of  a  critical  shear  force  model  has  also  been  tested  for  applications  in  the  cases  of 
drill-hole  mode  fractures  for  low  velocities.  Experimental  data  shows  that  for  the  osteoporotic 
bones,  the  average  force  of  resistance  is  1 1.68  lb.  for  0.250-in  diameter  spherical  projectiles 
while  it  is  12  lb.  for  0.406-in  diameter  projectiles  moving  at  200-400  ft/sec.  This  data  justifies  a 
critical  force  criterion  for  low  velocity  penetration  of  bullets  in  bones. 

2.5.5  Some  Additional  Notes  on  the  File  Structures  used  in  the  BONGEL  code 

There  are  two  main  modules  of  the  code.  Both  Codes  are  written  in  FORTRAN.  One  is  called 
BONGELT2.FOR  that  is  a  preprocessor  that  processes  the  user  input  and  generates  the  input 
code  for  the  second  module  BULLPENT.FOR.  This  module  actually  calculates  the  velocity¬ 
time,  space-time,  paths  of  the  projectile  for  the  front  tip,  back  tip  and  center  of  mass,  G. 
Currently,  the  code  works  for  a  bullet  or  any  other  slender  object  with  high  aspect  ratio  (e  g.,  a 
flechette)  as  well  as  for  spherical  projectiles  and  fragment  penetrations.  The  code  uses  a 
Microsoft  SYSTEM  command  through  I=SYSTEM(‘dos  command’C)  which  may  not  be 
available  to  other  compilers.  This  code  also  uses  two  FORTRAN  object  modules  called  CLS  and 
TLOC  ATE.  These  subroutines  are  used  to  clear  the  screen  and  to  locate  the  cursor  on  the  screen, 
respectively.  If  necessary,  the  code  can  be  modified  to  bypass  these  statements.  Descriptions  of 
the  internal  files  and  input  files  are  as  follows: 


®  R  D.  Eisler,  A  K.  Chatteijee,  and  G.H.  Burghart,  Algorithm  Development  to  Describe  Retardation  in  Human 
Tissue  of  a  19.6  Grain  Flechette,  Mission  Research  Corporation  Report  MRC-COM-R-93-0372,  US  Army  Natick 
Technical  Report  TR-95/023L,  May  1995  (Unclassified  -  Limited  Distribution). 
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2.5.5. 1  File  Description  for  BONGELT2  Code 

Files  opened  in  Bone-Gelatin  Interaction  Code:  B0NGELT2.F0R 


File 

Unit 

Description 

BONGELT2.IN 

1 

Input  File 

BONGELT2.DBG 

20 

Debug  File,  ifDEBUG.NE.O 

BULLIMP1.IN 

2 

Phase  1  Input  File  for  BULLPENT.for  Code 

BXYG1.PLT 

9 

First  Phase  X-Y  data  for  the  center  of  mass  G 

STPHS.COD 

15 

Projectile  arrest  flags 

BXGD.PLT 

10 

Velocity  vs.  space  data  for  G 

BULLIMP2.IN 

2 

Phase  II  input  file  for  BULLPENT 

BXYG2.PLT 

9 

Second  Phase  X-Y  data  for  G 

STPHS.COD 

15 

Projectile  arrest  flags 

BXGD.PLT 

10 

Described  earlier 

BULLIMP3.IN 

2 

Phase  III  input  file  for  BULLPENT 

BXYG3.PLT 

9 

Third  Phase  X-Y  data  for  G 

BULLIMP2.IN 

2 

Described  above 

BXY12.PLT 

9 

Front  Tip  X-Y  data  in  Phase  II 

BXY22.PLT 

10 

Back  Tip  X-Y  data  in  Phase  II 

BXYG2.PLT 

11 

X-Y  data  for  G  in  Phase  II 

XYZ1.PLT 

12 

Front  Tip  XYZ-data 

XYZ2.PLT 

14 

Back  Tip  XYZ-data 

XYZG.PLT 

16 

XYZ-data  for  G 

XYZALL.PLT 

18 

XYZ-data  for  Tipi  (Front),  Tip  2(Back)  and  G 

TXYZ12.PLT 

19 

Time  and  XYZ-data  for  Tip  1 ,  Tip  2  and  G 

Files  copied  through  SYSTEM  Copy  Command  in  BONGELT2 


l=SYSTEM{'COPY  BULLIMP1.IN  BULLPENT.IN'C) 
l=SYSTEM('COPY  BXY1  .PLT  BXY1 1  PLTC) 
l=SYSTEM('COPY  BXGD.PLT  BXGD1.PLTC) 
l=SYSTEM('COPY  BXY2.PLT  BXY21.PLTC) 
l=SYSTEM('COPY  BXYG.PLT  BXYGI.PLT'C) 

C  l=SYSTEM('COPY  BXG.PLT  BXG1  .PLTC)  commented  out 

l=SYSTEM(’COPY  BULLIMP2.IN  BULLPENT.IN'C) 
l=SYSTEM('COPY  BXY1.PLT  BXY12.PLT'C) 
l=SYSTEM('COPY  BXGD.PLT  BXGD2.PLTC) 
l=SYSTEM('COPY  BXY2.PLT  BXY22.PLTC) 
l=SYSTEM('COPY  BXYG.PLT  BXYG2.PLTC) 

C  l=SYSTEM('COPY  BXG.PLT  BXG2.PLTC)  commented  out 

l=SYSTEM('COPY  BULLIMP3.IN  BULLPENT.IN'C) 
l=SYSTEM(’COPY  BXY1.PLT  BXY13.PLTC) 
l=SYSTEM('COPY  BXGD.PLT  BXGD3.PLTC) 
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l=SYSTEM('COPY  BXY2.PLT  BXY23.PLTC) 
l=SYSTEM(’COPY  BXYG.PLT  BXYGS.PLT'C) 

C  l=SYSTEM(’C0PYBXG.PLTBXG2.PLTC)  commented  out 

l=SYSTEM('COPY  BULLIMP2.IN  BULLPENT.IN’C) 
l=SYSTEM('COPY  BXY1.PLT  BXY12.PLTC) 
l=SYSTEM('COPY  BXGD.PLT  BXGD2.PLTC) 
l=SYSTEM('COPY  BXY2.PLT  BXY22.PLTC) 
l=SYSTEM(’COPY  BXYG.PLT  BXYG2.PLTC) 

C  l=SYSTEM('COPYBXG.PLTBXG2.PLrC)  commented  out 

2.5.52  File  Description  for  Bullet  Penetration  Code  BULLPENT 


File 

Unit 

Description 

BULLPENT.IN 

1 

Input  File 

BULLPENT.DBG 

8 

Debug  File 

BXY1.PLT 

9 

XY-File  for  Front  Tip 

BXY2.PLT 

10 

XY-File  for  Back  Tip 

BXYG.PLT 

11 

XY-File  for  the  center  of  mass  G 

BYAW.PLT 

12 

Yaw  Data 

BULLPEN.ANM 

14 

Animation  Data  File 

BXDX.PLT 

15 

Velocity  vs.  Space  for  G 

STPHS.COD 

16 

Data  File  for  Arrest  Flags 

BXG.PLT 

2 

X-Time  File  for  G 

BXGD.PLT 

3 

X-Velocity-Time  File  for  G 

BYG.PLT 

4 

Y-Time  File  for  G 

BYGD.PLT 

5 

Y-Velocity-Time  File  for  G 

BTH.PLT 

6 

Yaw-Time  File 

BTHD.PLT 

7 

Yaw  Rate-Time  File 

2.5. 5. 2. 1  Coordinate  System  Used  in  BONGEL  Code 


With  reference  to  Figure  28,  a  description  of  the  two  reference  frames  that  have  been  used  in  the 
development  of  BONGEL  Code.  All  coordinate  input  to  the  BONGEL.IN  file  must  be  made  in 
the  user  or  global  frame.  F.  But  the  code  uses  some  internal  frame  so  that  other  modules  of  the 
code  can  be  used.  This  is  referred  to  as  the  local  frame  (Qx’y’z’)  with  origin  at  the  bone  impact 
point  Q.  A  description  of  the  shotline  in  the  local  coordinate  frame  and  the  impact  direction  for 
the  rebound  mode  is  shown  in  Figure  29. 

Frame  F:  (Oxyz)  Global  or  User  Defined  Coordinate  System 

(Qx’y’z’)  Local  coordinate  system  at  the  bone  impact  point  Q 

(£,  m,  n)  direction  cosines  of  the  original  bullet  shot  line  with  respect  to  the  frame  F. 

{a,p,y)  direction  cosine  of  the  impulsive  force  at  Q  with  respect  to  F.  This  direction  is  along  the 

outward  normal  direction  to  the  bone  surface  at  the  impact  point  Q. 

Plane  ;r  is  the  plane  containing  the  lines  {l,m,n)  and  {a,p,Y). 
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>  X 

usi:r  frame: 

Figure  28.  Coordinate  Systems  Used  in  BONEGEL  Code 
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SHOT  LINE 

OXYZ  (USER  FRAME) 

Figure  29.  Shotline  and  Impact  Direction  in  the  Local 
Frame  System 

(^i,  mj,  Mj)  direction  cosines  of  the  new  direction  of  velocity  of  the  center  of  mass  G  right  after 
bone  impact  in  the  case  where  rebounding  occurs  after  bone  impact.  This  forms  the  new  x’  axis 
at  Q.  This  is  also  referred  to  as  LP,MP,NP  in  the  BONGELT  code. 

Cross  product  of  the  lines  and  {a,/3,/)  forms  the  new  z’  axis  at  Q. 

The  direction  cosines  of  the  Qz’  axis  are  denoted  as 

Cross  product  of  z’-  and  x’-axes  at  Q  forms  the  new  y’  axis  at  Q. 

The  direction  cosines  of  the  Qy’  axis  are  denoted  as  (^3,  m3,  ^3) 


2.5.5.2.2  Generation  of  Internal  Files  in  BONGELT  and  Other  Called  Programs  (RK2S)  and  Codes 
(BULPENT2) 


RK2S  Subroutine  (For  One 

dimensional  motion 

X  =  x(t) ,  X  -space,  t  -time) 

File 

Unit 

Description 

SXY.PLT 

9 

xvs.  t 

SXYD.PLT 

10 

dx 

—  versus  t 
dt 

SYYD.PLT 

11 

dx 

X  versus — 
dt 

STPHS.COD 

15 

Code  Control  Parameter 
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BULLPENT2  Code  (Run  by  BONGELT2  Code,  for  Two-dimensional  Motion 
X  =  x(t),y  =  yit) ;  Subscript  G  indicates  the  center  of  mass  of  the  bullet,  subscript  1 

indicates  the  front  tip  while  subscript  2  indicates  back  tip  of  the  bullet) 

0  indicates  bullet  orientation  with  the  x-axis. 


File 

Unit 

Description 

BXG.PLT 

2 

Xq  vs.  t 

BXGD.PLT 

3 

*"3  VS.t 

dt 

BYG.PLT 

4 

Vg  vs.  t 

BYGD.PLT 

5 

vs.( 

dt 

BTH.PLT 

6 

^(deg)  vs.t 

BTHD.PLT 

7 

de 

—  vs.  t 
dt 

BULLPENT2.DBG 

8 

Debug  File 

BXY1.PLT 

9 

ix,y)  for  tip  1 

BXY2.PLT 

10 

{x,y)  for  tip  2 

BXYG.PLT 

11 

(x,  y)  for  the  center  of  mass  G 

BYAW.PLT 

12 

Yaw  angle  vs.  time  t 

BULLPEN.ANM 

14 

Graphic  Animation  File 

STPHS.COD 

16 

Code  control  parameter 

BXY1D.PLT 

18 

Velocity  components  for  the  front  tip 
along  the  x-  and  y-axis. 

Note:  These  file  named  has  been  modified  by  adding  1, 2  or  3  after  their  names  to  indicate  the 
data  obtained  during  phase  1, 2  or  3  of  the  projectile  penetration  where  the  different  phase  of 
motion  are  defined  as  follows: 

Phase  1:  Projectile  motion  prior  to  bone  impact 

Phase  2:  Projectile  rebound  mode  or  penetration-through-bone  mode 

Phase  3:  Projectile  reentry  in  tissue  medium  after  bone  penetration  is  completed. 

2.5.6  Units  Used  in  the  BONGEL  Code 

Unless  otherwise  specified,  the  following  units  are  used  for  all  input  and  output  data: 

Length  (L)  in  centimeters  (cm) 

Time  (T)  in  seconds  (s) 

Mass  (M)  in  grams  (gm) 
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All  other  units  are  consistent  with  the  above  units.  For  example,  velocity 


^  —  I  is  measured  in 
T 


cm 

s 

etc. 


,  acceleration 


..  cm  j,  „(ML^.  .  gm.cm  .  ,  -  gm.cm 

IS  m  — r- ,  force  F  is  m  - — r — ,  impulse  (F.r )  is  m  - 

yr  J  s 
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3  BALLISTIC  TESTING 

Table  3  shows  a  complete  log  of  ballistic  test  implemented  under  this  program.  The  purpose  of 
the  ballistic  testing  is  fourfold:  (1)  determine  for  a  spherical  projectile  the  threshold  velocity 
required  to  penetrate  -  i.e.,  inelastically  deform  -  various  formulations  of  gelatin  targets  with 
parametrically  varied  mechanical  properties,  (2)  collect  data  relative  to  the  effect  of  leading  edge 
projectile  geometry  on  target  transient  response,  (3)  collect  data  relative  to  target  boundary 
effects;  and,  (4)  segregate  pressure  drag  versus  lateral  surface  viscous  effects. 

Item  1  was  investigated  by  using  spherical  projectiles  with  two  diameters  and  the  same  material 
density.  The  striking  velocity  required  to  enter  10, 15,  and  20%  gelatin  targets  was 
experimentally  determined.  Since  concentrations  less  than  20%  are  sensitive  to  temperature 
variations,  testing  was  accomplished  at  controlled  temperatures. 

Items  2,  3,  and  4  involved  investigating  the  transient  response  of  the  gelatin  for  two  different 
gelatin  formulations  and  variable  projectile  and  target  geometries  (gelatin  blocks  of  different 
sizes  and  with  different  support  conditions).  This  phase  of  the  experimental  effort  involved 
collecting  high  frequency  pressure  measurements  supplemented  with  temporally  correlated  flash 
photography. 

An  analytical  model  of  dynamic  pressure  was  developed  to  support  investigation  of  issues  2, 3, 
and  4.  The  dynamic  pressure  model  describes  the  early  time  distribution  of  pressure  promoted  by 
interaction  of  the  projectile  and  target  medium.  By  modeling  attenuation  characteristics  of  the 
target  material  and  experinientally  determining  pressures  at  some  distance  away  from  the  wound 
tract,  the  dynamic  pressure  model  was  used  to  determine  time  resolved  pressure  distributions 
local  to  the  projectile.  This  enabled  segregation  of  base  pressure  drag  from  lateral  surface 
viscous  effects  in  the  retardation  algorithm  and  contributed  to  development  of  a  non-linear  spring 
lattice  that  simulates  projectile-soft  tissue  interaction  and  response.  The  data  also  enable  a 
physical  interpretation  of  parameters  governing  retardation  of  the  projectile.  This  physical 
interpretation  of  model  parameters  is  necessary  to  extend  models  initially  being  developed  to 
describe  projectile  interaction  with  biosimulants,  to  human  tissue. 

From  previous  work  (U.  S.  Army  Natick  Research,  Development,  and  Engineering  Center 
Contracts  DAAK60-92-C-0003  and  DAAK60-92-C-0008)  an  algorithm  describing  the  retarding 
force  per  unit  mass,  R(v),  as  a  function  of  projectile  velocity  for  a  sphere  in  20%  ordnance 
gelatin  was  developed.  This  algorithm  had  the  form  of  an  n-dimensional  polynomial  in  terms  of 
the  instantaneous  velocity  of  a  spherical  projectile;  i.e.,  R(v)  =  a  +  6v  +  cv^  +  ... .  The  constant 
coefficients,  a,  h,  c, ... ,  were  empirically  determined  and  were  specific  to  the  projectile 
geometry-target  medium  combination. 

In  this  effort,  a  second  order  polynomial  was  used  to  describe  penetration  of  a  sphere  into  20% 
ordnance  gelatin.  In  this  polynomial,  the  a  and  cv^  terms  can  be  shown  to  be  of  similar 
magnitudes  and  identically  equal  for  a  projectile  velocity  of  873  fps.  If  the  retarding  force  is 
divided  by  the  presented  area  of  the  projectile.  A,  in  the  direction  of  motion  to  obtain  pressure 
immediately  in  front  of  the  projectile,  the  implausible  result  that  pressure  is  of  similar 
magnitudes  at  very  high  and  very  low  projectile  velocities,  is  obtained.  This  problem  was 
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Table  3.  Summary  Descriptions  of  Ballistic  Tests  Implemented 


INSrRUM&JTATiON 
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fzyh^i'g'ipais" 


LOCATION  RANGE  PRESSURE 
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. 
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_ 

L 
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2.4 
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~ 
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■ 

1 

II 

■  ■ 

*1 

i 

II 

II 

s 

& 
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1 
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370* 
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_ 

790* 

_ 

20% 

20% 

20%  (COOL) 
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r 
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f 
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corrected  by  segregating  dynamic  and  static  components  of  pressure.  That  is,  total  pressure  was 
scaled  by  the  relative  magnitude  of  the  dynamic  pressure  that  must  be  a  maximum  at  the  time  of 
projectile  entry  into  the  target.  This  is  shown  below 

(  ^(^)Y  ^  Y JL Y 

I,  A  \pvj  A  Yvo; 

where  vo  is  the  projectile  striking  velocity  and  p  is  the  target  media  density. 

Another  enhancement  included  addition  of  the  following  term 


to  describe  pressure  attenuation  with  distance  and  dispersion  or  frequency  dependence  of 
attenuation.  In  this  expression  a  is  the  attenuation  coefficient,  oR  describes  attenuation  due  to 
geometric  scattering  where  is  radial  distance,  and  Q(co)  describes  dispersion  of  pressure  as  a 
function  of  frequency,  co.  In  the  dynamic  pressure  model  the  resulting  expression  is  integrated 
over  time  for  a  moving  line  source.  The  waves  arriving  from  the  various  sources  are  then 
accumulated  at  field  points  to  describe  time-resolved  pressure  at  a  particular  location  away  from 
the  wound  tract. 

The  approach  described  above  was  used  to  predict  pressure  away  from  the  projectile  trajectory 
for  the  experimental  configuration  shown  in  Figure  30.  Calibration  tests  were  conducted  using, 
pressure  transducers  that  were  not  cast  in  the  gelatin  (as  they  were  in  subsequent  tests)  but  were 
bonded  to  plates  that  rested  on  the  free  surfaces  of  the  gelatin.  (The  transducers  were  therefore 
held  in  compression  and  wave  reflection  between  the  free  surface  and  the  transducer 
configuration  promoted  occurrence  of  high  frequencies  in  the  traces  shown  in  Figure  31  and 
Figure  32). 


Figure  30.  Ballistic  Test  Configuration 

The  calibration  experiment  involved  ballistic  penetration  of  a  1.04  gram,  0.250-inch  diameter 
steel  sphere  striking  20%  ordnance  gelatin  at  1370  feet  per  second.  The  total  penetration  of  the 
sphere  was  7.0  inches.  Pressure  transducers  were  located  at  three  and  six  inches  from  the  struck 
edge  of  the  4.5  X  5  X  9.75-inch  gelatin  target.  The  xi  and  xa  dimensions  were  1.5  and  3.2  inches. 


68 


Mission  Research  Corporation 


DAMD17-94-C-4099 


- 

■ 

TR1  = 
1 

/ 

1 

0.  2\ 

S/W  /2 

■■O.Zr 

€ 

ns 

3" 

-  }20mr 

(Kfii 

1  1  t  < 

M  1  t 

i 

-j- 

- 

,TRS 

/ 

/ 

—  I 

>2  0.2 

/S2 

f 

/ 

1 

f 

/ 

1 

Figure  31.  PCS  Transducers  1.5  inches  from  Trajectory  at  48  jjsec 


The  purpose  of  this  calibration  experiment  was  to  define  time  synchronization  and  ifequency 
response  requirements  for  subsequent  experiments  since  instrumentation  circuitry  had  to  be  built 
and  the  frequency  response  (and  filter  strategy)  of  the  pressure  transducers  had  to  be  determined. 
Since  data  was  generated  however,  albeit  difficult  to  interpret,  the  data  was  used  to  determine  if 
order-of-magnitude  agreement  could  be  obtained  by  exercising  preliminary  models. 
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Figure  32.  PCB  Transducers  3.2  inches  from  trajectory  at  40  fjsec 
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Traces  from  the  four  transducers  are  shown  in  Figure  31  and  Figure  32.  Figure  31  shows  traces 
of  pressure  versus  time  from  sensors  located  on  a  plane  1.5  inches  from  the  trajectory  at  3  and  6 
inches  from  the  impact  surface  of  the  gelatin  target.  The  peak  pressures  associated  with  the 
initial  expansion  of  the  temporary  cavity  are  45  psi  at  3  inches  and  30  psi  at  6  inches.  These 
peaks  occur  at  0.6  and  1.8  msec,  respectively.  Similar  results  for  a  plane  3.2  inches  from  the 
trajectory  is  shown  in  Figure  32  where  peak  pressures  of  41  and  30  psi  is  developed  at  0.5  and 
0.6  msec. 

A  comparison  of  Figure  31  and  Figure  32  show  in  both  cases  an  initial  high  frequency  peak 
which  is  not  very  well  resolved  and  low  frequency  response  where  associated  amplitudes  are 
attenuated  between  10  and  30%.  The  analytical  model  predicts  peak  dynamic  pressures  at  these 
locations  between  40  and  80  psi.  The  prediction  assumed  that  the  attenuation  coefficient  was 
identically  equal  to  zero  and  the  velocity  of  the  displacement  field  normal  to  the  projectile 
motion  was  between  0. 1  and  0.2  of  the  projectile  velocity.  This  estimate  was  based  on  empirical 
observation  of  high-speed  movies  of  other  experiments.  As  explained  below  the  displacement 
field  normal  to  the  projectile  motion  is  critical  for  determining  the  initial  conditions  for  the  non¬ 
linear  spring  lattice.  Figure  33  and  Figure  35  show  predicted  pressure  versus  change  in  velocity 
and  projectile  velocity  versus  time. 


Figure  33.  Predicted  Pressure  versus  Change  in  Velocity  (vo  is  the 
striking  velocity  and  v  is  the  instantaneous  projectile  velocity) 


A  lattice  of  non-linear  springs  is  envisioned  to  describe  the  mechanical  response  of  soft-tissue. 
This  is  depicted  in  Figure  34.  The  high  frequency  portion  of  the  power  spectrum  promotes 
damage  and  the  energy  content  of  these  frequencies  dissipate  rapidly  due  to  conversion  of  energy 
into  mechanical  work.  The  initial  model  was  one-dimensional  where  spring  constants  were 
modeled  as  a  quadratic  equation:  kix  +  k2^  and  assumed  elastic  unloading.  This  model  was 
expanded  to  a  two  and  then  three-dimensional  lattice.  A  more  general  formulation  of  the  spring 
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Figure  34.  Schematic  of  Non-Linear  Spring  Model 
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constants  and  unloading  function  was  also  developed  based  on  the  mechanical  stress-strain  curve 
of  relevant  tissue.  Development  of  the  model  depicted  in  Figure  34  highlights  four  theoretical 
issues  discussed  in  items  /  through  iv  below. 

(i)  Development  of  the  Projectile  Retardation  Algorithm.  The  first  problem  concerns 
development  of  the  projectile  retardation  algorithm.  Two  insights  from  previous  work  at  MRC 
(DAAK60-92-C-0003  and  D  AAK60-92-C-0008)  demonstrated  that:  (a)  Close  approximations  of 
the  retarding  force  on  spheres  of  different  dimensions  can  be  obtained  by  knowing  the  retarding 
force  on  a  prescribed  sphere  and  scaling  the  force  by  the  ballistic  coefficient  of  the  unknovm 
sphere,  and,  (b)  For  certain  velocity  regimes,  ,  lift  and  drag  forces  on  a  projectile  penetrating  a 
tissue  simulant  can  be  determined  by  representing  a  projectile  as  an  ensemble  of  spheres  and 
integrating  the  retarding  force  on  each  sphere  over  the  geometry  of  the  projectile.  This  last 
problem  is  a  complicated  highly  non-linear  integration  since  the  center  of  projectile  rotation  and 
integration  limits  change  with  time.  These  two  insights  however  reduce  the  original  problem  of 
determining  the  detailed  forces  acting  on  a  projectile  during  penetration  (which  is  probably 
intractable)  to  the  problem  of  determining  the  retarding  force  per  unit  mass  on  a  sphere. 

Currently,  the  model  used  to  determine  the  retarding  force  on  a  sphere  is  semi-empirical  and  is 
described  in  detail  in  Section  2.  A  goal  of  the  current  effort  is  to  develop  the  retardation  model  in 
terms  of  physical  properties  of  the  projectile  and  target  media  so  that  the  model  can  be  used  to 
describe  what  occurs  in  a  living  person  with  the  appropriate  substitution  of  tissue  properties.  A 
problem  that  fhistrates  this  strategy  is  the  fact  that  target  properties  change  with  projectile 
penetration  velocity.  To  circumvent  this  difficulty,  properties  of  the  target  medium  at  limiting 
physical  states  (i.e.,  a  viscous  fluid  at  high  penetration  velocities  and  a  viscoplastic  solid  at  low 
penetration  velocities)  where  target  properties  are  well  defined  (and  a  database  exists  for  human 
tissue);  are  used  to  construct  a  relation  between  retarding  force  and  penetration  velocity. 

It  was  shown  in  Section  2  that  the  low  and  high  velocity  asymptotes  that  describe  projectile 
retardation  can  be  represented  as  an  infinite  series  where  the  first  term  in  each  series  dominates. 

It  was  further  shown  that  the  first  term  in  the  high  velocity  as)rmptote  could  be  related  to  the 
physical  density  of  the  target  medium  and  the  drag  coefficient  of  the  projectile.  Further,  the 
ultimate  mechanical  static  can  be  related  to  the  first  term  of  the  low  velocity  asymptote.  However, 
this  property(s)  was  not  well  defined  and  remained  to  be  established.  Without  this  set  of 
properties  being  rigorously  defined,  the  models  being  developed  in  this  effort  can  not  be 
confidently  extrapolated  to  human  tissue. 

(ii)  Determination,  as  a  Function  cf  Projectile  Velocity,  Initial  Particle  Velocity  Transverse  to 
Projectile  Motion.  This  problem  concerns  determination  of  the  initial  velocities  on  the  various 
chains  of  the  non-linear  spring  lattice  intercepted  by  the  projectile  (see  Figure  34).  A  videotape 
acquired  fromDr.  Ronald  Bellamy,  Col.,  USAMC  shows  an  AK74  projectile  penetrating  10% 
ordnance  gelatin  photographed  at  20,000  frames  a  second.  Unfortunately,  no  timing  information 
was  conveyed  on  the  tape  and  the  framing  rate  is  non-uniform  (since  it  ta^es  time  for  the  camera 
to  accelerate  to  the  prescribed  rate  of  20K  frames  per  second).  However,  it  is  estimated  that  the 
displacement  velocity  normal  to  the  projectile  trajectory  is  approximately  10  to  20%  of  the 
projectile  penetration  velocity. 
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Figure  37.  Schematic  and  Nomenclature  of  Boundary  Layer  Separation  from 
Sphere 


The  problem  of  developing  a  model  to  more  generally  describe  the  transverse  displacement  field 
is  alluded  to  in  Figure  37  where  a  schematic  of  a  sphere  penetrating  a  viscous  medium  is  shown. 
The  dashed  triangular  region  with  subtended  angles  y  and  6  represent  the  boundary  layer.  If  the 
transverse  particle  velocity  of  the  target  medium  can  be  represented  diS  y  =  Pv  where  v  is  the 

projectile  velocity  and  P  is  an  undetermined  coefficient  and  the  point  at  which  the  boundary 
layer  separates  from  the  surface  of  the  sphere.  The  boundary  layer  separation  can  be  further 
represented  by  the  subtended  angle  a  that  is  a  fi'ee  parameter.  The  a  and  P  parameters  can  be 
shown  to  be  related  as  shown  in  Figure  36.  Analytically  Predicted  P  =p  (a).  In  Figure  36. 
Analytically  Predicted  P  =P  (a),  for  P  =  0.2  (i.e.,  transverse  particle  velocity  of  target  media  is 
20%  of  projectile  velocity)  the  point  of  boundary  layer  separation  corresponds  to  subtended 
angle  of  a  =  0.02. 


ro 


P(a)  := 


1-x 

fufr  -  X 


dx 


Figure  36.  Analytically  Predicted  p=p  (a) 
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The  problem  then  reduces  to  two  issues.  The  first  issue  is  how  to  predict  a  for  different  projectile 
geometries,  striking  conditions,  projectile  velocities,  and  target  media.  Second,  knowing  the 
distribution  of  ct  —  ct  (y)  and  v  as  a  function  of  penetration  depth,  transverse  velocities  must  be 
mapped  onto  the  correct  chain  of  springs.  This  issue  is  related  to  the  viscosity  of  the  tissue  and 
how  much  the  boundary  region  stretches  and  sticks  to  the  surfaced  of  the  projectile  before  the 
boundary  layer  sheds.  Ballistic  experiments  were  designed  to  gain  insight  into  both  of  these 
issues  and  are  explained  below. 

Projectile  geometries  included  a  sphere,  conically  sharp  cylinder  that  was  confined  to  prevent 
yaw,  and  a  cylinder  launched  at  a  prescribed  yaw.  Each  of  these  geometries  was  expected  to 
produce  very  different  trends  in  transverse  target  particle  velocity.  For  example,  the  boundary 
layer  should  remain  attached  to  the  conical  surface  of  conically  sharp  cylinder.  In  the  expression 
y  =  ^  described  earlier  where  is  x  the  projectile  velocity  in  the  direction  of  the  projectile 

motion  and  yis  the  transverse  particle  velocity  promoted  in  the  target  by  the  projectile  motion,  P 
should  be  a  constant.  The  maximum  extent  of  the  temporary  cavity  and  the  extent  of  the  residual 
cavity  should  therefore  be  linearly  proportional  to  the  velocity  history  of  the  projectile.  In  the 
case  of  a  sphere  the  separation  point  will  be  a  fiinction  of  velocity  and  therefore  P  =  P(v).  The 
cylinder  will  be  the  most  challenging  geometry  in  that  the  rotational  kinematics  of  the  cylinder 
will  result  in  a  time  dependent  leading  edge  geometry. 

From  ballistic  experiments  using  military  bullets,  the  permanent  wound  tract  is  seen  to  be  highly 
non-uniform.  This  is  due  to  the  rotational  kinematics  of  the  projectile  changing  the  leading  edge 
geometry  and  presented  area  of  the  projectile  in  the  direction  of  motion  as  well  as  the  variation  in 
projectile  velocity.  Second,  there  are  a  lot  of  radial  fi-actures  associated  with  passage  of  the 
projectile.  The  presence  of  these  radial  fractures  highlights  two  additional  issues.  First,  how  does 
this  type  of  mechanical  damage  translate  to  the  human  body?  Second,  how  is  the  size  of  the 
temporary  cavity  measured.  Traditionally,  this  has  been  based  on  photographically  observing  a 
discontinuity  in  the  optical  properties  of  the  target  or  change  in  the  target  index  of 
refraction  over  time. 

This  may  be  very  misleading  however  since  the  fi-acturcs  form  at  about  the  bulk  sound  speed  of 
the  material  (about  5,000  ft/sec).  This  is  on  the  order  of  tens  of  microseconds.  The  temporary 
cavity  however  forms  on  the  order  of  milliseconds.  The  observed  changes  in  target  optical 
properties  are  therefore  very  likely  due  to  diffraction  patterns  caused  by  pre-existing  fracture 

planes  opening  and  closing.  These  diffraction  patterns  may  completely  obscure  observation  of 
the  temporary  cavity. 

These  difficulties  were  addressed  by  construction  of  an  array  of  light  tubes  that  will  back  light 
the  targrt  and  are  loca^ted  perpendicular  to  the  film  plane.  The  tubes  were  made  from  ABS  plastic 
with  stainless  steel  shim  stock  inserts.  The  tubes  were  3.5  inches  in  diameter  and  13  inches  long 
with  xenon  flash  lamps  (1  microsecond  on/off).  The  purpose  of  these  tubes  was  to  produce  an 
intense  highly  diffuse  light  source  that  will  tend  to  demphasize  diffraction  effects  when 
photographed.  Additionally,  gelatin  targets  were  constructed  with  specific  types  of  defects 
fabricated  into  the  targets.  Ifre  pre-damaged  blocks  were  photographed  with  the  constructed  light 
source  so  that  more  reliable  inferences  can  be  made  regarding  the  interpretation  of  photographed 
features. 
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(iii)  Tissue  Damage  Criteria.  The  residual  cavity  from  a  penetrating  injury  results  from  initial 
damage  promoted  by  the  distribution  of  pressure  associated  with  the  projectile  motion.  The 
damage  can  also  be  enhanced  after  passage  of  the  projectile.  This  later  damage  is  associated  with 
transient  displacement  of  the  target  medium  and  elastic  recovery  of  the  material  behind  the 
projectile.  In  both  cases  a  stress  field  is  promoted  which  can  strain  tissue  beyond  the  elastic  limit 
of  the  respective  tissue.  When  tissue  damage  occurs,  elastic  properties  of  the  tissue  can  be 
degraded  and/or  stress  concentrations  can  be  promoted.  This  alters  the  time-resolved  description 
of  the  stress  and  strain  fields  in  the  neighborhood  of  the  wound  tract.  The  accumulation  of 
inelastic  strain  over  time  however  governs  the  geometry  of  the  wound  tract  and  adjacent  tissue 
damage.  Therefore,  damage  criteria  are  required  for  the  various  tissues  affected  which  describe 
the  nature  of  the  damage  and  at  what  levels  of  stress/ strain  the  damage  occurs. 

This  issue  also  relates  to  the  initial  velocity  promoted  by  the  projectile  motion  in  the  transverse 
direction.  When  the  projectile  interacts  with  the  target,  material  stretches  around  the  projectile 
and  eventually  breaks.  The  tensile  stress  at  which  the  target  material  around  the  projectile 
eventually  breaks  can  be  analytically  related  to  the  transverse  particle  velocity  of  the  target 
material,  the  subtended  angle  on  the  projectile  surface  at  which  the  boundary  layer  separates,  and 
initial  conditions  on  the  various  chains  of  springs  in  the  lattice. 

(iv)  Tissue  Prestress.  In  vivo  tissues  are  generally  stressed.  If  an  artery  is  cut,  it  will  shrink  away 
from  the  cut.  A  broken  tendon  retracts  and  the  lung  is  in  tension  at  all  times.  Pre-stress  at  the 
time  of  wounding  may  significantly  influence  the  geometry  of  the  wound.  This  is  particularly 
true  since  biological  materials  deviate  from  Hookes  law  (i.e.,  the  relation  between  stress  and 
strain  is  non-linear).  At  a  minimum,  some  approximation  of  the  pre-existing  stress  state  of  the 
affected  tissues  may  have  to  be  superimposed  on  the  stress  field  resulting  from  the  projectile 
interaction  with  the  tissue.  This  pre-existing  stress  state  may  vary  with  activity  as  different 
muscles  and  tendons  react  during  task  performance.  At  a  maximum,  the  pre-stress  may  influence 
the  mode  of  damage.  This  is  likely  to  be  the  situation  with  bone  tissue  where  the  fracture  mode 
can  be  affected.  Unfortunately,  virtually  all  of  the  ballistic  testing  reported  in  the  literature  has 
been  done  on  unstressed  tissue.  In  the  case  of  autopsies,  the  activity,  and  therefore  kinematic 
state  of  the  various  muscles  and  tendons  at  the  time  of  wounding  is  not  known. 


3.1  Ballistic  Impact  of  Spherical  Projectiles  into  1 0%  and  1 5% 
Ordnance  Gelatin 

Tests  were  conducted  at  a  nominal  striking  velocity  of  1,400  feet  per  second,  using  1/4-inch 
spherical  steel  projectiles.  Ordnance  gelatin  targets  were  fabricated  in  the  shape  of  rectangular 
prisms,  measuring  nominally  5x5x10  inches.  All  but  two  blocks  were  cast  with  two  pressure 
probes  installed,  on  centerline  and  nominally  0.5  inches  above  the  projectile  trajectory.  Spacing 
from  the  striking  surface  was  1.50  and  4.20  inches  respectively.  The  extent  of  the  temporary 
cavity  was  assumed  to  be  given  by  the  dimensions  of  the  diffraction  pattern  recorded  on  film, 
using  backlighted  flash  photography.  Testing  started  testing  with  10%  gelatin  and  then 
progressed  to  15%  and  20%  gelatin. 
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Figure  38.  Shadowgram  to  show  scale  of  subsequent 
photos  (caliper  open  1  inch) 


Figure  38  was  taken  to  establish  proper 
scaling  of  detail  in  the  photographs.  It 
shows  a  backlighted  image  of  calipers, 
opened  to  1 ,000  inches,  placed  in  the 
position  of  the  projectile  trajectory.  The 
rectangular  patterns  visible  in  the 
photograph  are  present  on  the  optical 
random  diffuser  used  in  the  illumination 
beam  and  are  not  used  in  these 
experiments. 

The  photographs  in  Figure  39  through 
Figure  41  are  from  different  shots  but 
may  be  viewed  as  the  progression  of  the 
development  of  the  temporary  cavity, 
which  has  been  defined  here  as  the 
boundaries  of  the  recorded  optical 
diffraction  patterns  in  10%  ordnance 
gelatin.  Figure  39  shows  the  image  of 
the  two  pressure  probes  in  the  gelatin 
target,  before  impact.  Figure  40  was 
recorded  106  microseconds  after 
impact,  with  the  projectile  clearly 
visible  directly  below  the  first 
transducer.  The  first  transducer  showed 
a  pressure  pulse  arriving  35 
microseconds  after  impact  and 
transducer  two  showed  a  pressure  rise 
55  microseconds  later,  or  16 
microseconds  before  the  flash  recorded 
the  picture  of  Figure  40.  The  arrival 
times  of  these  pressure  pulses  were 
consistent  with  sound  speeds  in  the 
target  material.  Figure  41  was  recorded 

287  microseconds  after  impact;  the  first  transducer  showed  a  rise  47  microseconds  after  impact 
and  transducer  two  responded  66  microseconds  later,  or  somewhat  slower  than  the  sound  speed. 
Based  on  the  photo  in  Figure  41,  the  temporary  cavity  has  nearly  engulfed  the  first  transducer  at 
that  time  and  the  transducer  shows  near  zero  pressure. 


Figure  39.  Shadowgram  at  0  microseconds  in  10% 
ordnance  gelatin 


We  then  progressed  to  15%  gelatin.  It  had  become  clear  that  the  temporary  cavity  would 
continue  to  expand,  well  after  the  times  that  were  documented  and  therefore  was  not  of  interest. 
The  photo  of  Figure  42  was  recorded  286  microseconds  after  impact  on  a  15%  gelatin  block  at  a 
somewhat  lower  velocity  (nominally  1,200  Ft/sec).  Pressure  pulse  arrival  times  were  29  and  71 
microseconds  later.  It  was  noted  with  interest  that  the  presence  of  the  pressure  probe  influences 
the  growth  of  the  temporary  cavity.  To  address  this  issue,  we  attempted  to  record  a  late  time 
photo  of  the  temporary  cavity  that  would  engulf  the  probes.  The  photograph  of  Figure  43 
resulted  when  the  flash  was  fired  2,095  microseconds  after  impact. 
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Figure  40.  Shadowgram  at  106  micorseconds  in 
10%  ordnance  gelatin 


Figure  41.  Shadowgram  at  287  microseconds  in 
10%  ordnance  gelatin 


The  influence  of  the  probes  on  the 
geometry  of  the  temporary  cavity  is 
spectacular.  Also  it  is  clear  that  contact 
between  the  gelatin  and  probes  cannot  be 
assured  for  late  time  events.  Pressure 
traces  recorded  after  the  arrival  of  the 
temporary  cavity  are  clearly  not  useful. 

To  document  the  growth  of  the  temporary 
cavity  in  the  absence  of  probes  we 
recorded  the  shadowgram  of  Figure  44, 
also  2,095  microseconds  after  impact,  into 
an  uninstrumented  gelatin  block.  It  seems 
that  the  diameter  would  continue  to  grow; 
we  estimate  to  2.0  inches  maximum  or  8 
projectile  diameters. 

Measurements  of  the  permanent  wound 
track,  unfortunately,  are  difficult  because 
the  target  material  easily  tears  and  also 
changes  shape  due  to  moisture  loss.  Figure 
45  contains  back  lighted  photographs  of 
the  permanent  wound  tracks,  one  with 
pressure  probes  and  the  other  without. 
There  appear  to  be  no  significant 
differences  between  the  two  tracks.  It  can 
then  be  concluded  that  the  effect  of  the 
presence  of  the  pressure  probes  is  limited 
to  the  formation  of  the  temporary  cavity. 


Figure  46  was  obtained  during  a  check-out  shot.  Here  we  used  a  previously  tested  gelatin  block 
with  a  permanent  wound  track  and  a  spherical  projectile  imbedded  in  it,  turned  on  is  side.  The 
shadowgram  was  recorded  2,000  microseconds  after  impact.  The  temporary  cavity  looks  quite 
similar  to  the  late  time  cavity  of  Figure  42,  except  near  the  rear  (exit)  surface  of  the  block.  It 
seems  that  the  free  surface  effects  have  created  the  necked  down  temporary  cavity  structure 
visible  and  are  also  responsible  for  the  appearance  of  the  dark  ring,  concentric  with  projectile 
trajectory  and  visible  on  the  rear  surface.  The  cause  of  the  ring’s  appearance  is  currently  not 
definitively  known. 


It  is  also  interesting  to  observe  that  the  previously  straight,  pre-existing  wound  track  was 
displaced  on  the  order  of  1/4  inch  in  the  direction  of  the  bullet  travel.  Thus,  the  gelatin  is  not 
only  moving  radially  outward  to  form  the  temporary  cavity,  but  also  longitudinally,  in  the 
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Figure  42.  Shadowgram  in  15%  gelatin  recorded 
286  micorseconds  after  impact 


direction  of  projectile  travel,  as  discussed  in 
Section  2  relative  to  the  “sticking 
hypothesis”  of  the  target  medium  on  the 
projectile. 

The  absolute  magnitudes  of  the  pressures 
recorded  are  probably  not  meaningful  except 
for  the  very  first  response  of  the  pressure 
probes.  After  that,  the  large-scale 
displacement  of  the  material  around  the 
probes  and  probe  mounts  will  cause  local 
and  temporary  separation  of  the  material 
from  the  transducers,  followed  by  the 
unpredictable  collapse  of  these  cavities.  We 
are  comfortable  only  reporting  the 
magnitude  of  the  first  pressure  spike  that 
arrives  at  each  probe.  The  rest  of  the  traces 
may  not  be  useful.  Unfortunately  we  have 
also  consistently  underestimated  the 
magnitudes  of  the  early  time  spikes  so  that 
the  traces  initially  exceeded  the  range 
available  in  the  storage  medium.  From  these 
traces  we  have  only  retrieved  arrival  times 
of  the  pressure  pulses  at  the  traducers,  as 
correlated  with  the  entry  event  (ballistic 
screen)  and  firing  of  the  Xenon  flash 
producing  the  shadowgrams. 


For  run  5-22-01,  which  produced  the  late 
time  image  of  the  temporary  cavity,  as 
shown  in  Figure  42,  we  reduced  the 
sensitivity  of  the  storage  scope  to  reveal  the 
peaks  as  seen  by  the  two  pressure  probes. 

Figure  47  shows  the  screen  event  of  the 
projectile  entering  the  block  on  the  lower  trace  and  the  response  of  the  first  transducer.  The 
elapsed  time  from  screen  to  transducer  is  29  microseconds  which  is  consistent  with  the  sound 
speed  of  the  material.  The  amplitude  of  the  first  peak  is  nominally  900  psi  and  the  second  peak  is 
1,200  psi.  The  early  arrival  time  of  pressure  pulses  of  this  magnitude  is  not  consistent  with  our 
expectations  for  striking  velocities  well  below  the  sound  speed  in  the  target  material.  The  total 
duration  of  the  pressure  pulse  is  on  the  order  of  140  microseconds  which  is  probably  dominated 
by  reflected  waves  from  the  boundaries  of  the  block. 


Figure  43.  Shadowgram  in  15%  gelatin  with 
pressure  transducers  and  recorded  2.095  msec 
after  impact 


Figure  48  shows  data  from  the  same  test;  the  lower  trace  shows  the  time  firing  of  the  flash 
correlated  with  the  response  of  the  second  transducer.  The  pressure  pulse  shape  is  quite  different 
from  the  first  pulse  and  the  maximum  pressure  reached  is  1,350  psi.  that  is  probably  to  be  taken 
as  equal  to  the  first  pulse.  The  difference  in  the  profile  shape  is  probably  attributable  to  the 
influence  of  reflected  waves  from  the  sidewalls.  The  flash,  which  yielded  the  photo  of  Figure  42 
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Figure  44.  Shadowgram  recorded  2. 095  msec 
after  impact  in  15%  gelatin  without  pressure 
transducers 


uncertainties  from  events  that  do  not  have  a 
very  crisp  rise  and  are  recorded  at  rather  slow 
writing  speeds. 

It  seems  quite  clear  from  the  above 
considerations  that  high-pressure  pulses  are 
observed  that  are  transmitted  at  or  near  sound 
speed  and  are  not  significantly  attenuated  over 
distances  on  the  order  of  inches.  The  physical 
significance  of  pressure  indications  after  this 
postulated  shock  passage  is  unclear.  The 
arrival  of  the  boundary  of  the  temporary  cavity 
at  the  transducer  face  does  not  appear  to  be 
accompanied  by  a  significant  pressure  rise. 
Pressures  measured  at  late  time,  typical  of 
maximum  temporary  cavity  extent  are  near 
zero  and  slightly  negative;  i.e.,  sub- 
atmospheric.  As  stated  previously,  at  the 
present  time  we  cannot  be  sure  that  these 
measurements  have  physical  significance. 

It  is  expected  that  the  response  of  the 
instrumented  gelatin  targets  may  be  analogous 
to  what  can  be  expected  when  bone  is  present, 
surrounded  by  large  amounts  of  soft  tissue,  like 
a  thigh.  If  the  projectile  does  not  impact  the 
femur,  then  initial  damage  to  the  femur  must 
be  the  result  of  the  pressure  loading  from  the 
striking  event.  If  the  earlier  claim  of  minimal 
attenuation  over  short  distances  is  valid,  then 


is  seen  to  have  occurred  much  later  than  the 
pressure  pulses,  2.020  milliseconds  after  pulse 
arrival  at  the  second  transducer.  Thus  the 
growth  of  the  temporary  cavity  is  seen  to  lag  far 
behind  the  pressure  pulse. 

The  data  presented  in  Figure  47  and  Figure  48 
is  temporally  correlated  by  an  independent 
measurement  (not  shown)  of  the  elapsed  time 
between  ballistic  screen  and  flash;  in  this  case, 
2.095  milliseconds.  The  time  between  pressure 
pulse  arrival  at  probes  one  and  two  is  calculated 
by  subtracting  two  large  numbers,  yielding  an 
interval  of  41  microseconds.  This  would 
suggest  faster  than  sound  speed  transmission 
but  is  more  likely  attributable  to  measurement 


Figure  45.  Residual  cavities  in  15%  ordnance 
gelatin.  (Top)  in  presence  of  pressure 
transducer,  and  (Bottom)  without  probes. 
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the  probability  of  a  femur  fracture  is  not  strongly  dependent  on  close  proximity  of  the  femur  to 
the  wound  tract.  Rather,  the  striking  conditions  seem  to  dominate.  Based  on  very  limited  data  at 
this  point,  it  seems  that  a  pressure  pulse  on  the  order  of  1,000  psi  with  tens  of  microseconds 
duration  might  cause  a  bone  fracture.  Further,  analysis  and  testing  is  required  to  substantiate  this 
however. 


Figure  46.  Late-time  cavity  in  15%  gelatin  target 
with  pre-existing  residual  cavity  orthogonal  to 
trajectory 
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Figure  48.  Traces  from  Flash  Detector  and  Pressure  Transducer  2 


Figure  49.  Spherical  projectile  with 


impact  velocity  of 490 fps  at  1.1  msec  Figure  50.  Spherical  projectile  with  impact  velocity  of 
after  impact  1,870 ^s  107  microsec  after  impact 
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3.2  BALLisTic  Impact  of  Spherical  Projectiles  into  20%  Ordnance 
Gelatin 


Figure  52.  Spherical  projectile  with  impact 
velocity  of 1,300 jps  220  microsec  after  impact 


Figure  52.  Spherical  projectile  with 
striking  velocity  of  1,870 jps  107 
microsec  after  impact 


Testing  was  initiated  on  20%  gelatin  blocks, 
some  with  and  some  without  pressure 
probes.  The  target  geometry  remained  the 
same,  namely  a  rectangular  block,  5  by  5  by 
10  inches,  with  the  first  and  second  pressure 
probes  located  2.4  and  5.0  inches 
downstream  of  the  impact  surface.  The 
photographic  scale  remained  identical  to  the 
fiducial  shown  in  Figure  38.  Some 
photographs  include  a  shadow  image  of  a 
vertical  needle.  This  fiducial  is  located  2.00 
inches  downstream  from  the  impact  surface. 

The  shadowgraph  of  Figure  49  shows  a 
typical  diffraction  pattern  resulting  from  a 
low  impact  velocity  in  20%  gelatin.  The 
appearance  of  the  temporary  cavity  is  quite 
similar  to  typical,  early  time,  cavities 
produced  in  10%  and  15%  gelatin  with  little 
detail  visible.  A  more  revealing  picture  is 
shown  in  Figure  50  for  a  significantly  higher 
impact.  The  similarity  with  the  late  time 
pictures  of  Figure  43  and  Figure  44  for  15% 
gelatin  is  apparent.  Basically,  the  extent  of 
the  cavity  is  a  little  smaller  in  20%  gelatin  as 
expected  but  the  scale  of  "turbulence",  i.e., 
the  size  of  the  diffraction  cells  is  similar.  A 
post  test  photograph  (Figure  51)  recorded  at 
the  same  location  reveals  the  permanent 
cavity  which  is  also  similar  to  the  15%  data 
of  Figure  45. 

It  was  pointed  out  in  that  pressure  response, 
as  measured  by  the  probes  installed  in  the 
gelatin  targets  occurs  quite  early  in  the 
penetration  process.  It  was  therefore 
reasonable  to  attempt  to  record  the  early  time 
response  of  the  temporary  cavity  coinciding 
with  the  pressure  rise  and  fall  indicated  by 
the  pressure  transducers.  The  temporary 
cavity,  recorded  0.22  msec,  after  impact  is 
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Figure  54.  Spherical  projectile  with  striking 
velocity  of 1,940 jps  133  microsec  after 
impact  (small  probe) 


i 
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Figure  55.  Time  resolved  pressure  pulse  for 
Figure  7  test 

the  different  amplification  and  magnification 
of  the  two  figures  but  show  good  agreement. 
From  these  observations  it  was  concluded 
that  the  size  of  the  probe  mounting  structure 
does  not  alter  the  measured  response  of  the 
transducer.  Testing  was  therefore  continued 
using  the  rigid  structures. 

To  complete  the  pictorial  presentation  of  the 
growth  and  collapse  of  the  temporary  cavity 
we  conducted  tests  at  very  early  and  very  late 
times.  Temporary  cavities  in  20%  gel  were 
recorded  without  pressure  probes,  at  a 
nominal  striking  velocity  of  1,000  ft/sec. 


nicely  depicted  in  Figure  52;  unfortunately  the 
pressure  transducer  failed  to  trigger  the  recording 
medium.  The  test  was  then  repeated  with  a  shorter 
delay  to  attempt  to  capture  the  shadowgram  near 
the  peak  pressure.  A  predictable  temporary  cavity 
was  recorded  in  Figure  53  and  the  recorded 
pressure  response  is  shown  in  Figure  55.  The 
transducer  started  to  respond  about  60  psec  before 
flash  and  reached  a  peak  value  of  940  psi,  about 
40  psec.  after  flash.  The  onset  of  the  pressure  rise 
is  consistent  with  shock  arrival  times,  based  on 
wave  speeds.  The  peak  pressure  value  is  typically 
observed  when  the  projectile  is  closest  to  the 
probe. 

The  mounting  structure  designed  and  built  for  the 
purpose  of  installing  and  protecting  the 
piezoelectric  transducers  from  pressure  pulses  that 
might  exceed  several  thousand  psi  or  direct 
particle  impacts  were  relatively  rigid  and  massive. 
To  investigate  possible  importance  of  these 
structures  in  influencing  the  gelatin  response,  a 
test  was  performed  using  a  minimal  supporting 
structure  on  one  of  the  transducers.  The  resultant 
shadowgram  of  Figure  54  was  recorded  133  psec 
after  impact,  and  shows  the  extra  light  structure  of 
the  transducer  mount.  The  temporary  cavity  in 
Figure  54  is  very  similar  to  those  shown  in  Figure 
52  and  Figure  53  with  large  probes.  The  pressure 
measurement  obtained  with  the  "light"  probe  is 
shown  in  Figure  56.  A  comparison  of  this  trace 
with  that  of  Figure  55  is  a  little  complicated  by 


Figure  56.  Time  resolved  pressure  pulse  for 
experiment  in  Figure  8 
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Figure  58.  Spherical  projectile  with  striking  velocity  of 
750 jps  88  microsec  after  impact 


Figure  58  shows  the  temporary  cavity 
recorded  88  psec  after  impact  and 
Figure  57  shows  a  cavity  recorded 
nearly  2.73  msec  after  impact.  The  early 
cavity  is  much  as  expected.  The  partial 
mirror  image  of  the  cavity  seen  on  the 
impact  surface  is  not  fully  understood 
however.  That  surface  also  seems  to  be 
showing  an  upheaval,  typical  of  a 
surface  wave  or  crater  formation. 

The  late  time  photo  shows  a 
substantially  collapsed  cavity,  which 
seems  to  have  closed  at  the  entrance 
site.  The  evident  fringe-like  rings  on  the 
impact  surface  would  also  seem  to  be 
the  result  of  a  surface  wave,  radiating 
outward  from  the  impact  site.  Further 
studies  would  seem  to  be  of  interest,  in 
particular,  since  impacts  on  live  tissue 
are  known  to  produce  significant  surface 
upheavals  and  ejecta. 


To  further  address  issues  of  wave 
propagation,  transit  times  and  arrival 
times,  a  single  test  was  conducted  with 
one  pressure  probe  in  the  forward  (2.4") 
in  depth  installation  and  a  second  probe 
mounted  flush  on  the  rear  surface  of  the 
gelatin  block.  The  response  of  the  first 
transducer  (standard  location)  is  shown 
in  Figure  59  to  be  similar  to  the  pulses 
from  similar  experiments.  The  trace 
shown  in  Figure  60  was  recorded  with 
the  transducer  on  the  rear  surface,  some 

5  inches  from  the  impact  site.  The  two  traces  are  temporally  correlated  by  the  flash  event.  The 
two  peaks  marked  in  Figure  60  are  consistent  with  the  arrival  times  of  the  first  compression  wave 
from  the  impact  and  the  reflection  of  that  wave  off  the  striking  surface.  The  significance  of  the 
amplitudes  of  anything  but  the  first  pulse  may  be  questioned.  A  schematic  of  the  temporal 
correlation  of  the  various  events  derived  from  these  pressure  traces  and  other  fiducials  used  in 
the  recording  is  shown  in  Figure  61. 


Figure  57.  Spherical  projectile  with  impact  velocity  of 
1, 150  jps  2. 730  msec  after  impact 


3.2. 1  Pressure  Measurements  from  Ballistic  Testing 


The  pressure  data  presented  here  was  recording  using  a  Kistler  piezoelectric  transducer.  Two 
pressure  ranges  were  used;  a  high  sensitivity  transducer,  calibrated  to  500  psi  and  a  high  range 
transducer,  calibrated  to  5,000  psi.  Typically,  these  types  of  transducers  can  be  over-ranged  by 
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Figure  59.  Time  resolved  pressure  pulse 
from  first  transducer 


factors  of  2  and  3.  All  pressure  data  presented 
here  was  recorded  1/2  inch  above  the  projectile 
trajectory,  unless  noted  otherwise,  with  the 
transducer 

installed  2.40  inches  downstream  of  the  impact 
site.  The  recorded  pressure  traces  are  basically 
similar  in  shape.  They  exhibit  an  initial  fast  rise 
to  some  level,  typically  less  than  1/2  of  the  later 
peak  value.  A  drop  follows  the  first  peak  and  then 
a  gradual  rise  to  a  second  peak  that  is  also  the 
maximum  pressure  recorded.  This  second  peak  is 
then  followed  by  an  expansion  to  zero.  It  is  this 
second  peak,  or  the  absolute  peak  value  recorded 
that  will  be  used  in  the  subsequent  presentation.  It 
can  also  be  said  that  the  generalized  description 
above,  applies  both  to  the  20%  and  10%  gelatin 
data. 

It  should  be  pointed  out  here  that  the  absolute 
accuracy  of  these  pressure  measurements  is  very 
difficult  to  quantify.  The  greatest  complication 
arises  from  the  fact  that  ordnance  gelatin  is  not  a 
fluid  under  normal  conditions.  Assuming  that 
even  if  a  well  coupled  interface  between 
transducer  and  gelatin  exists  at  the  onset  of 
testing,  it  is  unclear  that  this  interface  will  remain 
frilly  in  tact  during  initial,  and,  in  particular  later 
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Figure  6L  Temporal  correlation  of  pressure  histories 
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Figure  63.  Peak  pressure  is  proportional  to  striking  velocity 

squared 


Figure  62.  Time  resolved  pressure  for  spherical 
projectile  striking  10%  gelatin  at  1,570 fps 


shock  passages.  We  have  reasoned  that 
during  periods  of  increasing  pressure  it 
should  be  possible  to  sustain  a  well-coupled 
interface.  It  is  therefore  suggested  that 
rising  pressures  be  reasonably  well 
represented.  Rarefactions  on  the  other  hand 
are  probably  not  well  tracked.  A  totally 
confusing  picture  results  when  gross  motion 
of  the  gelatin  material,  adjacent  to  the 
transducer,  has  initiated.  The  material  now 
has  to  flow  past  or  around  the  obstacle 
caused  by  the  transducer  and  that  becomes 
a  troublesome  event  especially  for  a  semi¬ 
solid  material. 

The  pressure  profiles  recorded  often  exhibit 
rather  sharp  spikes.  The  transducers  are 
capable  of  responding  in  the  Megahertz 
regime  but  the  recording/storage  medium  is 
limited  in  its  storage  capacity.  A  limited 
number  of  data  points  can  be  stored  at  any 
one  time.  The  resultant  temporal  resolution 
of  the  recorded  data  is  therefore  dependent 
on  the  recording  rate.  Our  storage 
oscilloscopes  are  10-Megahertz  devices 
however  at  the  recording  speeds  typical  of 
-this  data,  much  lower  sampling  rates  are 
used  because  of  the  limited  memory 
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Figure  65.  Time  resolved  pressure  for  spherical  projectile  striking  20% 

gelatin  at  1, 760 fps 

capacity.  It  must  therefore  be  expected  that  certain  short  duration  peaks  from  the  transducers  will 
not  be  accurately  replicated  in  the  data  as  it  relates  to  peak  values.  For  the  interpretation  of  the 
present  data  this  does  not  appear  to  present  a  significant  problem.  The  pressure  peaks  used  here 
are  typically  of  rather  long  duration  and  profound  accuracy  is  neither  claimed  nor  required. 
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Peak  pressures  recorded  are  plotted  as 
a  function  of  striking  velocity  in  Figure 
63.  The  data  identified  by  circles  is  for 
20%  gelatin  and  triangular  symbols 
represent  10%  gelatin  data.  All 
pressure  data  was  recorded  1/2  inch 
removed  from  the  trajectory  of  the  1/4 
inch  spherical  projectile  except  for  the 
data  point  represented  by  the  triangle; 
this  data  was  recorded  closer  in, 
nominally  1/4  inch  removed  from  the 
trajectory.  This  point  falls  well  above 
all  other  data  that  we  have  attributed  to 
the  close  proximity  to  the  trajectory  for 
this  point.  If  this  trend  is  confirmed,  it 
will  be  very  desirable  to  systematically 
vary  the  distance  of  the  transducer 
from  the  trajectory,  for  constant 
striking  velocities.  Peak  pressure  data  recorded  for  10%  and  20%  gelatin  appears  to  coincide. 

The  data  scatter  in  data  is  large,  however  no  trends  can  be  detected  that  would  distinguish 
between  10%  and  20%  data.  A  series  of  tests  at  low  striking  velocities  were  also  performed  but 
their  inclusion  in  the  presentation  of  Figure  63  is  not  practical.  For  striking  velocities  up  to  440 
fl/sec.  we  recorded  pressures  less  than  20  psi  in  20%  gelatin  and  that  result  is  indicated. 

Since  the  recorded  peak  pressures  presented 
are  not  the  result  of  the  material  being  shock- 
loaded,  but  what  seems  to  be  related  to  the 
passage  of  the  projectile,  i.e.,  a  fluid- 
mechanical  pressure,  it  seemed  reasonable  to 
check  for  a  velocity  squared  dependence  of 
the  peak  pressure.  In  fluid  mechanics  this 
would  be  equivalent  to  a  dynamic  pressure, 
which  includes  a  velocity  squared  and  density 
term.  Data  scatter  would  obscure  any  possible 
density  differences  for  10%  and  20%  gelatin 
(specific  gravity  of  1.02  versus  1.06).  This 
forced  a  velocity-squared  curve  through  a 
well-anchored  point  in  the  middle  of  the 
densest  data  scatter  (no  other  justification)  and 
calculated  the  remainder  of  the  parabola.  This 
second  order  dependence  on  velocity  fits  the 
data  probably  no  better  than  an  arbitrary  linear  fit,  but  it  seems  more  gratifying  from  a  physics 
stand  point  for  the  reason  given  above  and  because  it  will  go  to  zero  in  a  consistent  fashion  and 
also  seems  to  lend  itself  to  extrapolation  to  higher  values  of  velocity. 


Figure  67.  Spherical  projectile  incident  at  1, 750 
jps  332  microseconds  after  impact 


Figure  66.  Spherical  projectile  incident  at  1,370 jps, 
320  microsec  after  impact 
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Figure  68.  Time  Resolved  Pressure  for  Glancing 
Impact  on  Instrumented  “Femur”  at  1750 jps 


Figure  69.  Spherical  projectile  incident  at  2,160  Jps 
178  microsec  after  impact 


It  does  not  seem  desirable  to  include  all 
recorded  pressure  traces  but  a 
comparison  of  typical  pressure  traces  is 
informative.  The  pressure  trace 
recorded  in  10%  gelatin,  from  Run  8- 
10-03,  for  a  striking  velocity  of  1,570 
ft/sec.  is  shown  in  Figure  62.  A  similar 
test,  at  a  slightly  higher  striking 
velocity  of  1,760  ft/sec.,  is  shown  in 
Figure  65  for  20%  gelatin  from  Run  7- 
3 1-03.  A  comparison  of  the  two  traces, 
which  is  somewhat  complicated  by  the 
fact  that  different  recording  speeds 
were  used,  shows  that  the  traces  are 
quite  similar  in  shape  and  duration, 
with  a  marginal  match  in  peak 
pressure. 

3.2.2  Ballistic  Test  Data  For 
Simulated  Femur 

A  preliminary  effort  was  made  to 
investigate  the  pressure  field 
experienced  by  a  simulated  bone, 
embedded  in  the  ordnance  gelatin.  We 
nominally  attempted  to  approximate 
femur  in  size,  with  the  gelatin  matrix 
being  representative  of  soft  tissue 
surrounding  the  femur.  The  femur  was 
simulated  using  a  1.50-inch  diameter 
aluminum  cylinder,  embedded  in 
gelatin  as  shown  in  Figure  64.  The 
cylinder  was  instrumented  with  a  single, 
flush  mounted,  piezoelectric  transducer 
on  the  curved  surface.  The  axis  of  the 
transducer  was  at  45  degrees  to  the 
projectile  trajectory  in  all  but  one  test, 
where  it  was  parallel. 

For  the  first  test  on  a  simulated  femur, 
the  trajectory  of  the  projectile  was  3/4 
inches  removed  from  the  femur  surface. 
It  yielded  the  temporary  cavity 
photograph  of  Figure  66,  recorded  320 
psec.  after  impact.  Unfortunately,  the 
pressure  level  was  not  sufficient  to 
trigger  the  recording  medium  so  that  no 
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pressure  was  recorded.  The  second 
attempt  represented  an  over¬ 
correcting  and  resulted  in  a 
glancing  impact  on  the  femur,  as 
seen  in  Figure  67.  The  pressure 
trace  recorded  for  this  test  is  shown 
in  Figure  68  and  seems  to  be 
revealing  a  structural  shock  from 
the  solid  impact,  superimposed  on 
the  pressure  trace.  The  response 
went  off-scale  at  an  indicated 
pressure  level  of  1,240  psi. 

The  test  was  then  repeated  two 
times  with  only  partial  success. 
Figure  69  shows  the  temporary 
cavity,  recorded  178  psec.  after 
impact,  for  a  striking  velocity  of 
2,160  ^sec.  Unfortunately,  the 
pressure  trace  presented  in  Figure 
71  shows  that  the  transducer  went 
into  saturation  at  an  indicated 
1,410  psi.  The  third  attempt 
yielded  a  good  photograph  of  the 
temporary  cavity,  but  once  again, 
the  transducer  failed  to  trigger  the 
storage  scope  so  that  the  pressure 
data  was  lost. 


Since  these  tests  were  considered 
exploratory,  we  felt  justified  in 
changing  a  basic  condition,  before 
we  ever  had  a  complete  success 
with  the  original  configuration. 

For  this  test  we  rotated  the  femur 
so  that  the  axis  of  the  transducer 
would  be  parallel  to  the  projectile  trajectory,  i.e.,  so  that  the  diaphragm  would  be  facing  the 
projectile  (perpendicular).  The  shadowgram  of  Figure  70  shows  a  typical  temporary  cavity  and 
the  pressure  trace,  shown  in  Figure  72,  depicts  a  rather  different  response  than  was  seen  from  the 
pressure  probes  tested.  A  preliminary  observation  is  that  the  pressure  experienced  at  this 
location,  with  an  indicated  peak  of  nearly  2,100  psi  and  a  very  rapid  rise  to  this  peak,  followed 
by  an  equally  sharp  drop,  is  typical  of  shock  loading  and  reflection.  The  large  femur,  compared 
to  the  transducer  or  pressure  probes  is  a  much  more  efficient  surface  for  shock  reflection,  at  this 
orientation.  This  may  explain  the  high  spike.  The  second  peak,  seen  on  the  pressure  probes,  is 
not  expected  to  be  dominant  here  because  of  the  much  greater  distance  from  the  trajectory  and 
the  90-degree  rotation  of  the  transducer. 


Impact 


Photograph 


Figure  71.  Time  Resolved  Pressure  for  2,160  jps  Striking 

Velocity 
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Figure  72.  Time  Resolved  pressure  for  Ballistic  Experiment  in 

Previous  Figure. 


3.3  Ballistic  Impact  of  Projectiles  with  Various  Geometries  and 
Constituent  Materials 


If  projectile  retardation  within  ordnance  gelatin  targets  can  be  viewed  as  a  fluid  mechanical 
event,  basically  dependent  on  a  drag  coefficient,  an  area,  and  a  dynamic  pressure  (proportional  to 
density  and  velocity  squared),  then  the  ballistic  coefficient,  3,  which  is  the  weight  of  the 


W 

projectile,  W,  divided  by  its  drag  coefficient,  Co,  and  representative  area,  A{P  = - ),  should 

Cj^A 

correlate  projectiles  that  differ  in  mass  and  geometry.  In  the  case  of  spherical  projectiles,  the 
drag  coefficient,  Cd,  should  be  constant  over  a  reasonable  size  range  so  that  the  scaling 
parameter  is  reduced  to  W/d,  where  d  is  the  diameter  of  the  projectile. 


This  3-scaling  approach  appeared  to  be  reasonably  successful  for  addressing  projectile  geometry 
related  issues.  However,  when  current  data,  produced  with  1/4-inch  diameter  steel  spheres  was 
scaled  for  comparison  with  17  caliber  BB  data,  it  was  found  that  about  a  10%  discrepancy  was 
present.  The  1/4-inch  spheres  penetrated  10%  deeper.  This  observation  prompted  us  to  carefully 
reexamine  our  present  data  as  well  as  older  BB  data. 
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Scale 


1  inch 


Figure  76.  Corrfiguration  of  In-Situ  Experiments  and 
Variance  between  Maximum  Penetration  Depth  and  Final 
Resting  Location 


The  results  of  this  investigation, 
shown  in  Figure  73  through 
Figure  75,  confirmed  that  no 
differences  could  be  attributed 
to  differences  in  preparation. 
Figure  73  shows  the  penetration 
data  for  BB's  into  20%  gelatin, 
all  at  temperatures  between  37 
and  44  degrees  F.  Data  from 
previous  tests  as  well  as  14  tests 
performed  under  this  program 
are  shown.  Old  and  new  data  are 
indistinguishable.  The 
penetration  depths  for  1/4-inch 
steel  spheres,  all  generated 
under  this  program,  are  shown 
in  Figure  74  and  are  basically 
unremarkable.  The  relatively 
large  data  scatter  near  zero 
penetration  depth  is  not 
unexpected. 

This  data  was  then  scaled  by 
ballistic  coefficient  to 
correspond  to  BB-conditions. 
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Table  4.  Depth  of  Penetration  Data for  V4  ”  Steel  Sphere 
in  10, 15,  and  20%  Gelatin 


V-STRIKING 

DEPTH 

V-RESIDUALl 

<FPS) 

m 

12/14/01 

20% 

MEm 

12/14/02 

20% 

1600 

8.0 

20% 

1370 

3/8/02 

20% 

1320 

6.8/ 7.0 

3/9/01 

20% 

1390 

7.1  /  7.2 

3/13/01 

15% 

1370 

9.4 

3/13/02 

15% 

8.8 

3/21/01 

■HtUI 

1440 

7.3 

3/21/02 

20% 

1340 

5/17/03 

15% 

1200 

8.8 

7/18/01 

20% 

1370 

7.2 

m-tnbyM 

20% 

160 

0.35 

20% 

205 

0 

8/8/02 

20% 

320 

8/8/03 

20% 

170 

0 

8/8/04 

20% 

440 

1.53 

8/8/05 

20% 

250 

0.54 

8/10/02 

10% 

1610 

4.8 

890 

1570 

4.8 

770 

1  8/17/01 

20% 

1940 

9.8  /  9.6 

20% 

930 

4.7 

20% 

0.8 

8/17/04 

20% 

750 

3.3 

8/17/05 

1150 

5.6 

1060 

4.8 

470 

370 

4.8 

<50 

8/21/04 

1230 

4.8 

8/21/05 

WBSSM 

840 

4.8 

370 

2690 

4.8 

790 

■Ri'gai 

1780 

8.7 

1020 

5.0 

isa 

1330 

6.6/6.4 

1  8/31/03 

mwm 

740 

N/A 

9.0 

1700 

8.8 

9/6/01 

20% 

1570 

8.3/8.1 

9/6/04 

20% 

840 

3.6 

1170 

1520 

1  9/12/03 

20% 

1550 

7.7/7.5 

470 

1.6 

9/18/02 

20% 

1350 

6.9/6.7 

9/18/03 

20% 

1510 

7.6/7.3 

This  scaled  data  is  superimposed  on  the  penetration  data  of  Figure  73  and  presented  as  Figure  75. 
The  BB-data  is  presented  as  solid  symbols  and  falls  about  10  per  cent  below  the  penetration  data 
for  scaled  spheres.  This  discrepancy  is  real  and  mustbe  taken  to  mean  that  P-scaling  is  not 
entirely  valid,  or  that  the  drag  coefficients  are  not  constant.  One  factor  contributing  to  this 
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Figure  77.  Experiment  9-12-03 


discrepancy  is  illustrated  in  Figure  76 
where  it  was  observed  that  the  final 
resting  location  of  the  projectile  does 
not  correspond  to  the  maximum  depth 
of  penetration.  The  reason  for  this  is 
that  the  target  material  recoils  as  the 
projectile  slows  down  and  causes  the 
projectile  to  rebound.  This  appears  to 
be  about  a  3%  effect  in  20%  gelatin 
penetration  depth  measurements  and 
substantially  less  than  this  in  10  and 
15%  formulations  where  penetration 
depth  is  greater  and  material  elasticity 
less. 


Figure  78.  Experiment  9-18-02 


Impacts  on  lower  concentration  gelatin 
produce  far  greater  penetration  depths 
than  impacts  in  20%  gelatin.  A 
limited  number  of  tests  performed  on 
15%  gelatin  came  to  rest  within  the 
10-inch  long  model  for  velocities 
below  1,400  ft/sec.  For  10%  gelatin 
the  problem  was  greater.  We  therefore 
decided  to  address  the  penetration 
issue  by  measuring  the  residual 
projectile  velocity  after  penetration  of 
a  block  of  gelatin  in  the  transverse 
direction  (4.80  inches).  The  data  from 
these  tests  is  tabulated  in  the  summary 
tabulation  (Table  4)  of  all  1 /4-inch 
steel  sphere  penetration  data.  Two 
tests  performed  using  1/4-inch  glass 
sphere  produced  data  that  should  not 
be  used  here;  the  high  velocity  test 
resulted  in  a  shattered  projectile  and 
the  low  velocity  test  falls  in  a  region 
of  large  data  scatter. 


3.3. 1  In-Situ  Velocity  of  the  Penetrating  Projectile. 

The  mathematical  model  to  describe  the  retardation  of  the  spherical  projectile  in  gelatin  has  used 
final  penetration  depth  data,  for  a  given  striking  velocity,  to  describe  the  retardation  process. 

This  approach  must  assume  that  the  instantaneous  velocity/penetration  depth  is  given  by  the  final 
penetration  depth  versus  striking  velocity.  To  test  the  validity  of  this  approach  we  recorded 
double  flash  photographs  of  the  projectile  during  the  penetration  process.  The  distance  between 
the  two  photographic  images  of  the  projectile,  properly  scaled,  and  the  known  temporal 
separation  of  the  two  flashes  yield  a  mean  velocity  over  the  distance  between  the  two  images. 
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Figure  79.  Experiment  9-18-03 


The  results  from  three  tests  that 
yielded  usable  images  are  shown 
in  Figure  78  through  Figure  79 
and  in  the  summary  of  runs 
performed  (Table  3  and  Table 
4),  as  runs  9-12-3  and  9-18-2 
and  -3.  The  instantaneous 
velocity/penetration  data  shows 
that  the  velocity  loss  per 
distance  of  gelatin  penetrated  is 
not  constant,  as  might  be 
deduced  from  an  examination  of 
the  final  penetration  versus 
velocity  plot  of  Figure  74.  At 
higher  velocities  the  velocity 
loss  per  distance  penetrated  is 
greater  than  at  low  velocities. 

3.3.2  Tests  Using  Conical 
Projectiles 


Figure  80.  Pressure  Pulse  from  Experiment  8-23- 
02 


Testing  using  conically  sharp,  cylindrical 
projectiles  was  initiated  on  20%  gelatin, 
instrumented  with  high-speed  pressure 
transducers.  One  of  the  obvious  problems 
encountered  here  was  the  expected 
pitching/yawing  motion  of  this  type  of 
projectile,  especially  after  entering  the 
gelatin  block.  To  alleviate  this  problem  we 
launched  projectiles  that  were  over  4  inches 
long  and  placed  the  muzzle  close  to  the 
gelatin  target  so  that  yawing  motion  would 
be  minimal  during  the  early  phases  of 
penetration.  The  first  test  (8-23-01)  yawed 
excessively  and  curved  away  from  the 
transducer.  Shot  8-23-02  produced  the 
pressure  trace  of  Figure  80  but  failed  to 
yield  a  photograph.  Also,  the  rod  nearly 
impacted  the  transducer  that  accounts  for 
the  high  peak  pressure  shown.  For  8-23-03 
the  projectile  pitched  up  and  struck  the 
transducer  mount,  as  seen  in  the  rather 
spectacular  picture  of  Figure  82. 


It  became  clear  that  an  even  longer  projectile,  with  very  little  free-flight  distance  would  be 
required  to  assure  minimal  pitch  and  yaw.  Run  9-26-01  produced  the  pressure  trace  of  Figure  81, 
but  failed  to  yield  a  photograph  or  a  striking  velocity.  A  peak  pressure  of  1,290  psi  was  recorded. 
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This  test  was  then  repeated  as  run  10- 
05-01,  using  identical  gun  conditions. 
This  test  yielded  a  velocity  of  1,710 
ft/sec,  and  the  photograph  of  Figure 
83.  The  flow  field  visible  about  the 
projectile  is  typical  of  a  sharp  cone  in 
a  fluid  medium.  The  temporary  cavity 
evolves  more  slowly  than  for  spheres 
which  is  also  consistent  with  travel 
through  a  fluid.  The  faint  image  to  the 
left  of  the  projectile  is  the  result  of  a 
late  time,  inadvertent  second  flash  and 
should  be  ignored. 


Figure  82.  Shadowgram  from  Experiment  8-23-03 


A 


>  zyo 


The  pressure  trace  recorded  for  this 
sharp  cone  test  is  shown  in  Figure  84 
and  again  in  Figure  86  after 
amplification.  A  peak  pressure  of 
1,010  psi  is  indicated.  The  shape  of  the 
pressure  profile  repeats  quite  well  the 
profile  shown  in  Figure  81  for  a 
similar  projectile.  Comparison  with 
sphere  data  shows  significant 
similarity,  as  expected,  if  the  response 
is  fluid  mechanical  in  nature.  The 
sphere-generated  pressure  profile  of 
Figure  85  is  somewhat  difficult  to 
compare  to  that  of  the  cone,  in  Figure 
84,  because  of  the  different  recording 
amplitudes  and  writing  speeds.  The 
amplified  cone  pressure  profile  of 
Figure  86  is  better  suited  for  this 
comparison.  Clearly,  the  peak  pressure 
agrees  very  well  with  the  peak 
pressure  recorded  for  a  sphere,  at 
nearly  identical  striking 


velocity.  The  pressure  profiles  are  also 
similar;  pulse  durations  are  nearly 
identical  as  are  the  times  to  reach  the 
peak  value.  The  primary  difference  is 
in  the  early  time  pressure  peaks  that 
are  typical  of  the  arrival  of  the  shock 
from  the  impact  of  the  projectile  on 
the  gel  block.  Higher,  more 

pronounced  peaks  are  seen  for  the  blunt  spherical  projectile  than  for  the  sharp  cone,  which  is 
consistent  with  wave  propagation/expansion  considerations. 


Figure  81.  Pressure  Pulse  from  Experiment  9-26-01 


Figure  83.  Shadowgram  from  Experiment  10-05-01 
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Figure  85.  Pressure  Pulse  from  Experiment  7-31-03 
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4  GRAPHICAL  INTERFACE 

4.1  SOFTWARE  INTERFACE 


The  injury  modeling  software  is  organized  into  four  separate  programs,  or  modules,  as  shown  in 
Figure  87  below.  The  programs  use  a  consistent  input/output  file  format  so  that  the  output  of  one 
module  can  be  loaded  as  input  to  the  next  module.  This  enabled  a  script  to  be  used  to  automate 
the  execution  of  the  four  modules. 


Figure  87.  Injury  Modeling  Software 


GUI 

This  graphical  interface,  written  in  €++,  lets  the  user  aim  a  generic  weapon  at  the  limbs  (upper 
arm  or  upper  leg).  It  uses  loadmodelQ  to  load  polygonal  descriptions  of  the  tissue  surfaces,  and 
uses  writebulletQ  to  write  out  the  projectile  parameters  to  a  file  (bullet  in)  that  bongel  can  read. 

A  description  of  the  graphical  user  interface  of  this  module  is  contained  in  the  next  section. 

Bongel 

This  module  contains  the  three  programs  {bongel,  bulpent,  and  cavloc)  for  calculating  the  shapes 
of  the  temporary  and  permanent  wound  cavities.  It  queries  the  healthy  model  for  material 
properties  of  the  tissue  along  the  projectile  path,  and  writes  out  a  polygonal  description  of  the 
surface  of  each  cavity.  A  complete  description  of  its  design  and  implementation  is  given  below. 
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Boolean 

This  Boolean  geometry  software,  written  in  C,  loads  descriptions  of  the  temporary  and 
permanent  cavities  (using  loadcavQ),  as  well  as  the  healthy  anatomy  models,  and  uses  a  Boolean 
subtraction  algorithm  to  remove  the  tissue  within  the  permanent  cavity.  It  also  marks  the  tissue  in 
the  temporary  cavity  as  devitalized  using  a  Boolean  intersection  algorithm.  It  writes  .mgi  files 
(MGI  Inventor-like  format)  that  can  be  read  by  mgiview  for  display  purposes. 

Mgiview 

Mgiview  is  a  graphical  display  utility  that  had  been  previously  developed  by  MusculoGraphics. 
For  this  project,  it  was  enhanced  to  display  two  windows,  one  containing  the  3D  surface  models 
of  the  temporary  and  permanent  cavities  (created  by  module  2),  and  one  containing  the  injured 
limb  model  (created  by  module  3).  It  reads  one  or  more  ".mgi"  files,  which  can  contain  polygons, 
colors,  and  textures. 


4.2  GRAPHICAL  INTERFACE 

The  graphical  user  interface  is  divided  into  three  windows,  as  shown  in  Figure  88  below.  The 
large  window  on  the  left  contains  a  3D  surface  model  of  a  human  figure.  The  model  is  low- 
resolution  so  that  it  can  be  manipulated  quickly  and  easily  (the  high-resolution  model  is  used  for 
the  resulting  injury  calculations).  Also,  in  this  window,  a  virtual  rifle,  represented  by  the  3D 
arrow  at  the  right  edge  of  the  window  in  Figure  88.  A  dotted  line  extends  from  the  tip  of  the 
arrow  to  the  human  figure.  A  cross-hair  target  indicates  the  point  of  contact  with  the  body  (in  the 
picture  below  the  line  contacts  the  middle  of  the  right  thigh).  The  dotted  line  represents  the  bullet 
path  preceding  impact  with  the  body,  and  the  cross  hairs  represent  the  impact  location  of  the 
bullet  on  the  surface  of  the  body.  The  user  can  change  the  view  of  this  model  interactively,  but 
the  rifle  and  body  move  together  during  the  view  manipulations. 

To  aim  the  rifle  at  a  different  portion  of  the  body,  the  user  moves  the  body  within  the  window  in 
the  lower  right  comer  of  the  screen.  This  window  displays  a  view  from  the  rifle  barrel  and  down 
the  bullet  path  towards  the  body.  The  cross  hairs  remain  in  the  center  of  the  window,  and  the  user 
can  translate  the  body  in  3  dimensions  to  change  the  impact  point  of  the  bullet.  As  the  body  is 
moved,  the  coordinates  of  the  gun  direction,  impact  location,  and  list  of  tissues  in  the  line  of  fire 
are  updated  in  the  upper  right  of  the  screen.  If  the  user  wants  to  model  an  injury  involving  bullet- 
bone  contact,  he  can  move  the  body  until  the  label  bone  appears  in  the  tissue  list.  The 
information  window  in  the  upper  right  also  contains  two  interface  elements  that  allows  the  user 
to  change  the  properties  of  the  bullet.  The  user  can  choose  from  one  of  four  bullet  types  (AK47, 
AK74,  Ml 6,  hollow  point  -  currently  only  the  AK47  is  implemented),  and  the  user  can  set  the 
speed  of  the  bullet  at  the  time  of  impact  with  the  body  (from  200  -  2000  feet/sec).  When  satisfied 
with  the  bullet  properties  and  direction,  the  user  can  press  the  shoot  button  to  write  the 
parameters  to  a  file  and  then  mn  a  script  that  will  execute  the  other  three  modules  of  the 
program.  After  the  three  modules  are  finished,  two  windows  will  be  displayed  on  the  screen,  one 
containing  the  3D  surface  models  of  the  temporary  and  permanent  cavities,  and  one  containing  a 
3D  model  of  the  injured  limb. 
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Figure  88.  User  Interface  for  the  INFLICT  module  of  Injury  code 


4.3  BOOLEAN  CODE  FOR  INJURY  CREATION 

As  described  above,  module  3  is  Boolean  software  that  uses  the  injury  cavities  and  healthy  limb 
model  to  create  an  injured  limb  model.  In  the  first  step.  The  volume  of  tissue  fi'om  the  healthy 
model  that  lies  within  the  permanent  cavity  is  removed  (subtracted)  from  the  model  (call  the 
result,  model  M).  The  second  step  is  to  determine  the  region  of  tissue  in  model  M  that  lies  within 
the  temporary  cavity,  and  to  label  it  as  devitalized.  This  is  done  with  a  Boolean  intersection 
algorithm.  Devitalized  tissue  is  usually  darker  and  stiffer  than  live  tissue,  so  its  color  and  texture 
are  changed  to  reflect  these  changes.  The  third  step  is  to  subtract  the  temporary  cavity  Ifom 
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model  M,  and  then  add  the  devitalized  tissue  model  to  it.  Boolean  addition  is  not  needed  for  this 
step  because  there  is  no  overlap  between  the  two  model  portions.  The  result  is  an  injured  limb 
model  containing  a  wound  cavity  surrounded  by  devitalized  tissue.  The  equations  defining  this 
Boolean  process  are  shown  below: 

=  one  tissue  (e.g.,  skin,  femur)  from  healthy  limb  model  containing  n  tissues 
H  =  healthy  limb  model  = 

0<j<n 


P  =  permanent  wound  cavity 
T  =  temporary  wound  cavity 

M.^H.-P 

injured  limb  model  =  ^  (A/j  -T)  +  (Mj  n  T) 

0<i<n 


To  perform  these  operations,  we  implemented  Boolean  subtraction  and  intersection  algorithms. 
The  algorithms  operate  on  two  closed  polyhedra  (PI  and  P2)  of  arbitrary  shape.  The  polyhedra 
can  contain  any  number  of  vertices  and  planar  polygons,  and  can  be  convex  or  concave.  Both 
algorithms  begin  by  calculating  the  intersections  of  each  polygon  in  the  first  polyhedron  (PI) 
with  the  polygons  in  the  second  polyhedron  (P2).  In  the  general  case,  a  polygon  either  lies 
completely  inside  the  other  polyhedron,  completely  outside  the  other  polyhedron,  or  it  intersects 
it  and  is  divided  into  two  or  more  pieces.  For  Boolean  subtraction,  the  set  of  polygons 
representing  PI  -  P2  are  collected  and  formed  into  a  polyhedron.  For  intersection,  the  polygons 
representing  PI  &  P2  are  collected  in  a  similar  fashion. 

Much  of  the  design  and  implementation  of  these  algorithms  was  performed  by  MGI  for  another 
project,  but  several  enhancements  were  necessary  in  order  to  use  the  code  in  this  project. 
Specifically,  the  numerical  robustness  of  the  geometric  intersection  code  was  increased  so  that  it 
would  work  on  models  with  vertex  coordinates  much  less  than  one  (the  models  are  in  units  of 
meters)  without  losing  accuracy.  Also,  the  error  checking  and  error  recovery  routines  were 
improved  so  that  the  code  better  recognized  ill-formed  polyhedra,  attempted  to  correct  the 
defects,  and  printed  informative  error  messages  when  it  was  unable  to  recover.  Finally,  several 
bugs  were  fixed  so  that  once  a  polyhedron  was  loaded,  several  Boolean  operations  could  be 
performed  on  it  in  sequence  without  having  to  write  it  to  a  temporary  data  file  and  then  reload  it 
between  operations. 


4.4  CODE  TO  ASSEMBLE  INJURED  LIMB  MODEL  WITH  TEXTURE  MAPS 


Much  of  the  work  for  this  task  had  to  be  skipped  for  budgetary  as  well  as  scientific  reasons.  With 
a  reduced  scope,  we  did  not  create  complete  3D  models  of  each  tissue  (skin,  muscle,  fat,  fascia, 
and  bone)  in  both  limb  models,  nor  were  realistic  texture  maps  created  for  each  tissue  type  for 
healthy  and  devitalized  conditions.  As  a  result,  the  limb  models  were  simplified  to  contain  just 
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models  of  the  soft  tissue  and  bone.  All  of  the  material  between  skin  and  bone  was  assumed  to 
consist  of  a  generic  soft  tissue.  We  believe  that  although  this  simplification  is  a  large  one,  it  does 
not  prevent  the  software  from  producing  realistic  injury  models  in  some  cases.  Much  of  the 
tissue  in  the  upper  arm  and  upper  leg  is  muscle,  whose  retardation  coefficients  are  similar  to 
those  of  generic  soft  tissue  model  employed. 


4.5  MODEL  QUERY  CODE 

The  primary  interface  between  the  bongel  code  in  module  2  and  the  healthy  limb  model  is  the 
material jxt()  function.  This  function,  written  in  C++,  is  used  by  the  bongel  code  to  determine 
the  material  properties  of  the  tissue  at  the  tip  of  the  bullet  as  it  travels  through  the  limb.  The 
prototype  for  this  function  is; 

int  material_at(double  modefpt[],  double  tissue_coeff[],  char  *tissue_type); 

This  function  takes  as  input  the  X,  Y,  Z  coordinates  of  the  point  at  which  you  want  to  know  the 
material  properties  (modelpt[]),  and  returns  as  output  the  retardation  coefficients  of  the  material 
(tissue_coeff[])  and  the  type  of  tissue  (tissue  type)  at  the  given  point.  Material  jxtQ  keeps  track 
of  the  previous  point  queried,  and  when  it  is  called  again,  it  forms  a  line  segment  from  the 
previous  point  to  the  current  point.  It  then  checks  to  see  if  this  line  segment  crosses  any  tissue 
boundaries.  If  it  does  not,  then  the  current  point  is  in  the  same  tissue  that  it  was  in  last  time.  If  the 
line  does  cross  one  or  more  boundaries,  the  sequence  of  boundary  crossings  is  used  to  determine 
the  current  tissue.  The  retardation  coefficients  and  type  of  this  tissue  are  then  returned  to  the 
bongel  code. 


4.6  MATERIAL  PARAMETERS 

The  material  properties  that  MRC  calculated  from  the  gelatin  experiments  were  the  same  as  the 
properties  that  bongel  needs  to  calculate  the  wound  cavities.  Thus  we  did  not  have  to  do  any 
further  processing  of  the  values.  We  inserted  them  into  the  healthy  model  file,  assigning  them  to 
the  skin  and  muscle  tissues.  When  bongel  queries  the  model  for  the  current  tissue  and  its  material 
properties,  these  values  are  returned  without  modification. 


4.7  CURVATURES  AND  SURFACE  NORMALS 

The  normal  vector  at  each  point  on  the  surface  of  the  skin  and  the  bone  is  needed  in  order  to 
calculate  the  bullet  trajectory  through  the  limb.  The  healthy  model  file  contains  normal  vectors 
for  each  polygon  in  each  3D-tissue  model.  These  normals  were  calculated  using  a  polyhedron 
processing  utility  that  had  been  developed  at  MusculoGraphics.  Surface  curvature  values  are  not 
currently  used  in  the  injury  modeling  software. 

Module  1  uses  the  surface  normal  at  the  point  of  impact  between  the  bullet  and  the  skin  to 
calculate  several  parameters  that  are  written  to  BULLET.IN,  and  subsequently  used  by  module  2. 
The  code  that  calculates  the  intersection  between  the  line  of  fire  and  the  skin  model  returns  the 
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impact  point  as  well  as  the  polygon  in  which  the  impact  occurred.  The  normal  vector  for  this 
polygon  can  then  be  looked  up  from  the  list  stored  in  the  healthy  model  file.  This  normal  vector, 
the  line  of  fire,  and  a  vector  in  the  plane  of  the  impacted  polygon  are  used  to  calculate  THD  and 
PHD,  which  describe  the  angles  between  the  reference  frame  of  the  bullet  and  the  reference 
frame  of  the  skin  and  bone  model. 


4.8  GRAPHICAL  DISPLAY  OF  INJURED  LIMB 

The  work  for  this  task  centered  on  mgiview,  a  graphical  display  utility  that  had  been  previously 
developed  by  MusculoGraphics.  For  this  project,  it  was  enhanced  to  display  two  windows,  one 
containing  the  3D  surface  models  of  the  temporary  and  permanent  cavities  (created  by  module 
2),  and  one  containing  the  injured  limb  model  (created  by  module  3).  It  reads  one  or  more  "  mgi" 
files,  which  can  contain  polygons,  colors,  and  textures.  Due  to  funding  limitations,  the  tissue 
models  and  texture  maps  that  are  displayed,  while  accurately  representing  the  injury  shape,  do 
not  present  a  realistic  image  of  a  penetrating  injury. 


4.9  Integration  of  healthy  limb  models  into  boolean  code 

3D  geometric  models  of  the  skin  and  bone  of  the  upper  arm  and  leg  were  obtained  from  the 
Visible  Human  male.  These  models  were  processed  with  a  commercially  available  polygon 
editor  to  make  sure  that  they  had  no  holes  or  seams,  and  that  the  polygons  were  all  triangles  with 
counter-clockwise  vertex  ordering  and  outward-facing  normals.  They  were  then  put  into  an 
Inventor-like  format  so  that  they  could  be  loaded  into  each  of  the  first  three  modules  of  the  injury 
modeling  software.  Shown  below  is  an  abridged  version  of  the  model  file. 

#Inventor  V2 . 1  ascii 
Separator  { 

Separator  { 

Label  { 

label  "Body" 

} 

DEF  Tissue  Separator  { 

DEF  Bone  Group  { 

Label  { 

label  "Humerus" 

} 

MaterialProperties  { 

coefficients  [  0.00329  0.0  5660000.0  ] 

} 

Material  { 

ambientColor  111 
diffuseColor  111 

} 

Separator  { 

Coordinates  { 

point  [  “0.08689  0.14564  0.09006,  ...  “0.18034  0.11287 

0.95479  ] 

} 

ShapeHints  { 
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vertexOrdering  COUNTERCLOCKWISE 

) 

Normal  { 

vector  [  0.337718  -0.821392  -0.459632,  ...  0.255227  - 
0.962267  -0.0943473  ] 

} 

NormalBinding  { 

value  PER_VERTEX_INDEXED 

} 

IndexedFaceSet  { 

coordindex  [  0,  1,  2,  -1,  . . .  31,  30,  25,  -1  ] 
normalindex  [  0,  1,  2,  -1,  . . .  62,  52,  51,  -1  ] 


} 

} 

DEF  Tissue  Separator  { 

DEF  Muscle  Group  { 

Label  { 

label  "Upper  Arm" 

} 

MaterialProperties  { 

coefficients  [  0.284  -11300.0  148000000.0  ] 

} 

Material  { 

ambientColor  1  0.5  0.5 
diffuseColor  1  0.5  0.5 

} 

Separator  { 

Coordinates  { 

point  [  0.17224  0.11403  0.13076,  ...  0.13139  0.14543 

0.14347  ] 

} 

ShapeHints  { 

vertexOrdering  COUNTERCLOCKWISE 

} 

Normal  { 

vector  [  -0.660695  0.739354  -0.129756,  ...  -0.456261 
0.838466  -0.297993  ] 

} 

NormalBinding  { 

value  PER__VERTEX_INDEXED 

} 

IndexedFaceSet  { 

coordindex  [  0,  1,  2,  -1,  ...  35,  34,  29,  -1  ] 
normalindex  [  0,  1,  2,  -1,  ...  65,  55,  54,  -1  ] 

} 

} 

} 

} 

DEF  Tissue  Separator  { 

DEF  Bone  Group  { 

Label  { 

label  "Femur" 

} 

MaterialProperties  { 

coefficients  [  0.00329  0.0  5660000.0  ] 
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Material  { 

ainbientColor  111 
diffuseColor  111 

} 

Separator  { 

Coordinates  { 

point  [  -0.17939  -0.59382  0.00347,  ...  -0.16927  -0.63917 

0.01092  ] 

} 

ShapeHints  { 

vertexOrdering  COUNTERCLOCKWISE 

} 

Normal  { 

vector  [  0.337718  -0.821392  -0.459632,  ...  0.255227  - 
0.962267  -0.0943473  ] 

} 

NormalBinding  { 

value  PER_VERTEX_INDEXED 

} 

IndexedFaceSet  { 

coordindex  [  0,  1,  2,  -1,  ...  31,  30,  25,  -1  ] 
normalindex  [  0,  1,  2,  -1,  ...  62,  52,  51,  -1  ] 


} 


} 

DEF  Tissue  Separator  { 

DEF  Muscle  Group  { 

Label  { 

label  , "Thigh" 

} 

MaterialProperties  { 

coefficients  [  0.284  -11300.0  148000000.0  ] 

} 

Material  { 

ambientColor  1  0.5  0.5 
diffuseColor  1  0.5  0.5 

} 

Separator  { 

Coordinates  { 

point  [  -0.17472  -0.58372  0.00852,  ...  -0.16927  -0.60936 

0.01638  ] 

} 

ShapeHints  { 

vertexOrdering  COUNTERCLOCKWISE 

} 

Normal  { 

vector  [  -0.660695  0.739354  -0.129756,  ...  -0.456261 
0.838466  -0.297993  ] 

} 

NormalBinding  { 

value  PER_VERTEX_INDEXED 

} 

IndexedFaceSet  { 

coordindex  [  0,  1,  2,  -1,  ...  35,  34,  29,  -1  ] 
normalindex  [  0,  1,  2,  -1,  ...  65,  55,  54,  -1  ] 
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} 


} 


} 

DEF  Tissue  Separator  { 

DEF  Skin  Group  { 

Label  { 

label  "Whole  Body  Skin" 

} 

MaterialProperties  { 

coefficients  [  0.284  -11300.0  148000000.0  ] 

} 

Material  { 

ambientColor  0.7  0.4  0.2 
diffuseColor  0.7  0.4  0.2 

} 

Separator  { 

Coordinates  { 

point  [  -0.508422  0.11458  0.55278,  ...  -0.47785  0.14844 

0.53186  ] 

} 

ShapeHints  { 

vertexOrdering  COUNTERCLOCKWISE 

} 

Normal  { 

vector  [  0.849744  -0.104314  -0.516773,  ...  0.326578  - 
0.133913  0.935636  ] 

} 

NormalBinding  { 

value  PER_VERTEX_INDEXED 

} 

IndexedFaceSet  { 

coordindex  [  0,  1,  2,  -1,  ...  193,  192,  184,  -1  ] 
normalindex  [  0,  1,  2,  -1,  ...  288,  287,  266,  -1  ] 


} 


} 


4.10  Port  MRC  code  to  SGI  and  integrate  into  MGI  code 

The  injury  modeling  code  developed  by  MRC  works  as  three  separate  programs  {bongel, 
hulpent,  and  cavloc).  The  Fortran  source  code  for  each  of  these  programs  was  ported  to  the  MGI 
SGI  workstation  (Indigo2  Max  Impact  running  ERIX  6.2).  The  port  required  minimal  changes  to 
the  source  code  and  it  was  a  straightforward  process  to  compile  and  run  the  programs  on  the  SGI. 

The  next  step  was  to  make  changes  to  the  code  so  that  it  worked  in  conjunction  with  modules  1 
and  3  as  outlined  in  the  software  design.  The  interactive,  command-line  interface  that  prompted 
the  user  for  bullet  characteristics  was  removed  and  instead  all  of  the  necessary  input  variables 
were  defined  in  a  single  data  file  (which  is  generated'by  module  1).  The  assumptions  in  bulpent 
were  then  addressed,  i.e.,  that  the  limb  model  is  a  right,  circular  cylinder  of  soft  tissue  with  a 
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right,  circular  cylinder  of  bone  in  the  middle  of  it.  Instead  of  making  these  assumptions,  the  code 
was  changed  to  load  a  user-defined  3D-limb  model,  and  query  it  for  tissue  shapes  and  material 
properties.  This  query  code  is  also  used  to  determine  whether  or  not  the  bullet  hits  the  bone.  A 
description  of  this  query  code  is  included  above. 

The  final  change  that  was  made  was  to  alter  the  format  of  the  injury  cavity  files  that  are  output 
by  cavloc  and  loaded  into  module  3,  the  Boolean  software.  Cavloc  was  enhanced  to  create  closed 
polyhedra  containing  counter-clockwise  triangles  for  both  the  permanent  and  temporary  cavities. 
These  files  could  then  be  loaded  directly  into  the  Boolean  software  without  further  processing. 

The  injury  depicted  in  Figure  89  and  Figure  90was  generated  using  a  bullet  with  cross-sectional 
radius  of  0.254cm  and  an  impact  velocity  of  2000  m/s.  The  bullet  was  aimed  such  that  it  did  not 
impact  the  femur.  It  had  enough  speed  to  exit  the  thigh;  hence  there  is  an  exit  wound  as  well  as 
entrance  wound. 

The  wire-fi'ame  picture  (Figure  89)  shows  the  permanent  and  temporary  cavities  clipped  to  the 
surface  of  the  thigh.  The  wider  section  of  the  temporary  cavity  is  at  the  entrance  wound.  Because 
the  bullet  did  not  have  time  to  tumble  within  the  soft  tissue  of  the  thigh,  the  permanent  cavity 
remains  small.  Figure  90  shows  the  thigh  with  tissue  within  the  permanent  cavity  removed,  and 
tissue  within  the  temporary  cavity  colored  darker  than  the  rest.  The  view  is  of  the  entrance 
wound. 


Figure  89.  Wireframe  model  of  predicted  temporary  and  permanent 
cavities 
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5  CONCLUSIONS  AND  RECOMMENDATIONS 

This  effort  used  different  configurations  of  ordnance  gelatin  as  surrogate  materials  to  develop  an 
analytical  and  experimental  methodology  that  describes  mechanical  damage  to  human  tissue 
from  penetrating  wounds.  The  assumptions  employed  were  specific  to  the  lower  extremities  and 
battlefield  threats  but  could  be  relaxed  to  generalize  the  methodology  to  other  body  regions, 
other  threats,  and  other  types  of  wounds.  A  software  package  based  on  this  methodology  was 
developed  to  predict  wound  geometries.  This  software  was  interfaced  with  a  virtual  anatomy 
based  on  the  Visible  Human  database  to  describe  tissue  damage  relative  to  the  anatomy  of  the 
lower  extremities.  The  software  predicts  the  extent  of  soft  tissue  damage  as  well  as  bone 
interaction,  including  fracture  mode  and  ejected  bone  fragments. 

Recommendations  for  future  work  are  sixfold  and  do  not  encompass  major  technical  obstacles. 
First,  the  existing  soft  tissue-bone  interaction  models  should  be  enhanced  to  include  the  effect  of 
tissue  tearing  from  the  presence  of  ejected  bone  fragments  in  the  projectile  wake.  The 
importance  of  this  effect  is  significant.  See  for  example  the  discussion  of  the  casualty  described 
in  Figures  1  and  2  of  Appendix  A. 

Second,  the  analytical  models  developed  to  describe  wound  geometry  in  soft  tissue  are 
sufficiently  detailed  that  they  could  be  easily  extended  to  characterize  the  nature  of  the  soft  tissue 
damage  adjacent  to  the  wound  tract.  For  example,  since  the  models  describe  the  stress/strain  field 
promoted  by  the  projectile  interaction  with  soft  tissue,  the  predicted  stresses  and  strains  could  be 
used  to  describe  the  radius  over  which  blood  vessels  of  different  diameters  rupture,  the  extent  of 
tissue  concussion,  etc. 

Third,  the  mechanical  property  database  for  the  constituent  tissues  could  be  more  highly  resolved 
in  terms  of  number  of  tissue  categories  considered.  The  various  models  developed,  in  most  cases, 
are  functions  of  independently  measured  target  properties  where  the  data  is  already  available  but 
merely  has  to  be  reduced  properly.  In  particular,  it  would  be  desirable  to  have  a  model  specific  to 
skin. 

Fourth,  in  order  to  make  the  wound  visualization  more  compelling,  improved  texture  maps  of 
damaged  tissue  and  bleeding  need  to  be  developed.  The  bleeding  should  not  only  consider  free- 
flowing  blood  but  also  blood  that  is  not  free-flowing  and  has  altered  the  appearance  and 
mechanical  properties  of  other  tissue. 

Fifth,  bone  trauma  from  penetrating  wounds  is  a  major  source  of  battlefield  morbidity  and  long 
hospital  stays,  yet  very  little  is  known  about  the  mechanics  of  projectile  bone  interaction. 
Presumably,  just  as  the  pathogenesis  of  any  disease  is  studied  to  develop  remedial  treatment,  the 
same  logic  would  be  applicable  to  bone  trauma.  An  effort  in  this  area  would  presumably  include 
ballistic  experiments  on  pre-loaded  human  bones  in  various  gelatin  matrices  with  in-situ 
instrumentation. 

Finally,  there  is  a  need  to  more  extensively  validate  the  models  that  have  been  developed,  in 
particular,  by  comparison  with  human  autopsy  data.  - 
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1.  BACKGROUND 

The  overwhelming  source  of  morbidity  on  the  battlefield  is  penetrating  injuries  to  the 
musculoskeletal  system.  Wounds  to  lower  extremities,  in  particular,  accoimt  for  the 
majority  of  total  and  nonfatal  military  casualties  in  both  conventional  European-style 
campaigns  and  jungle  warfare.  Physical  characteristics  of  these  woxmds,  however,  are 
highly  variable  and  thwart  assessment  in  terms  of  their  medical  consequences  and 
operational  effects  on  soldier  performance.  This  is  illustrated  Ity  the  two  combat 
casualties  shown  in  Figures  1  and  2. 

Both  casualties  were  struck  in  the  forearm  Ity  bullets  fired  from  an  AK47  assault  rifle  at 
ranges  of  about  20  meters.  For  both  casualties  the  bullets  were  probably  of  identical 
design  and  construction  and  had  the  same  striking  velocity  and  kinetic  energy.  The 
casualty  shown  in  Figure  1  however  sustained  a  relatively  superficial  soft  tissue  woimd 
while  the  casualty  shown  in  Figure  2  almost  sustained  a  traumatic  amputation  of  the 
arm. 

This  subject  effort  is  providing  an  integrated  software  package,  applicable  to  the 
extremities,  that:  (1)  Simulates  tissue  damage  from  small  arm  projectiles  and  fiagments 
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Figure  1.  Casualty  with  perforating  gunshot  wound  of  the  forearm,  Fig,  4S2,  Ref  1, 


Fig.  2a,  Casualty  with  perforating  gunshot  wound  of  the  forearm,  Fig.  4-53,  Ref  L 

penetrating  the  musculoskeletal  system;  (2)  Performs  a  detailed  biomechanical 
assessment  of  residual  physical  capability;  and,  (3)  Provides  a  high  resolution  graphical 
interface  for  medical  simulation  and  training. 

The  program  started  in  September  1994  and  has  a  36  month  period  of  performance  and 
consists  of  five  tasks.  The  first  task  develops  models  describing  soft  tissue  damage  from 
a  penetrating  wound.  This  description  includes  the  permanent  wound  cavity,  tissue 
transient  response,  and  hemorrhagic  zone  of  extravasation.  Models  developed  in  the 
second  task  describe  interaction  of  bone  tissue  with  the  missile  including  mode  of 
fracture  produced. 


A-2 


Figure  2b.  Exit  wound  for  Casualty  shown  in  Figure  2a,  Fig.  4-54,  Ref  1. 


Figure  2c.  Roentgenogram  of  casualty  shown  in  Figure  2 
shattered  elbow,  Fig.  4-55,  Ref  1. 


Task  3  consists  of  three  subtasks.  The  first  subtask  develops 
models  to  describe  coupling  between  the  ballistic  response 
of  the  bone  and  the  surrounding  soft  tissue.  This  coupling 
occurs  due  to  ejection  of  secondary  bone  fragments  (Figure 
2c),  fragmentation  of  the  projectile,  and  enhanced  projectile 
yaw. 


In  the  second  subtask,  an  interface  will  be  developed  which  will  allow  visualization  of 
the  wound.  This  visualization  will  be  used  in  the  final  subtask  which  involves  selecting 
casualties  from  an  autopsy  data  base  for  simulating  specific  wounds.  This  will  be  done  to 
refine  models  and  certify  model  fidelity. 


Task  4  integrates  biomechanical  software  which  describes  the  moment  generating 
capability  of  joints  given  the  initial  and  deformed  geometry  of  the  bones,  muscles  and 
tendons  resulting  from  the  projectile  interaction.  The  final  task.  Task  5,  is  an 
administrative  task  related  to  preparation  of  deliverables. 


The  models  from  this  effort  will  eventually  be  imported  into  a  lower  extremity  battlefield 
trauma  simulator  under  development  by  Musculographics,  Inc.  for  ARPA.  The  simulator 
will  be  used  to  train  combat  medics  in  the  management  of  lower  extremity  trauma. 
During  1996-97,  these  simulators  will  begin  replacing  the  use  of  animals  in  trauma 
management  training  of  combat  medics. 

The  current  paper  describes  progress  during  the  first  year  which  focuses  on  the  first  two 
tasks,  simulation  of  soft  tissue  response  and  bone  trauma. 


2.  PHENOMENOLOGY 

The  gross  morphology  of  a  typical  penetrating  wound  fi’om  a  projectile  into  skeletal 
muscle  tissue  can  be  thought  of  as  divided  into  three  zones.  First  there  is  the  permanent 
wound  cavity  which  is  a  hole  containing  blood  clots,  detached  chunks  of  tissue,  and 
foreign  material  that  may  have  been  aspirated  or  carried  into  the  wound.  Smrounding 
the  permanent  wound  cavity  is  a  zone  of  extravasation  which  is  grossly  homorrhagic, 
shredded,  pulped  muscle  around  the  permanent  cavity.  Farther  away  is  the  zone  of 
concussion  which  is  grossly  normal  muscle  but  with  histiological  abnormalities  such  as 
interstitial  hemorrhage  and  abnormal  myocytes. 

Traditional  soft-tissue  wound  surgeiy  excises  the  hemorrhagic,  shredded,  pulped  tissue 
of  the  zone  of  extravasation  and  some  of  the  neighboring  zone  of  concussion.  Skeletal 
muscle  tissue  is  surgicaUy  excised  if  its  color  is  abnormal,  if  it  fails  to  contract  when 
simulated,  if  its  circulation  is  impaired  (i.e.,  not  bleed  when  cut),  and  if  it  has  abnormal 
consistency.’ 

Immediately  upon  projectile  impact  with  tissue,  high  amplitude  stress  waves  radiate 
fi’om  the  impact  region.  Figure  3  shows  a  schematic  of  a  ballistic  test  series  recently 
completed  on  ordnance  gelatin  targets.  Ordnance  gelatin  of  various  formulations  is 
conventionally  used  to  simulate  skeletal  muscle  tissue  during  ballistic  testing. 

In  our  experiments,  gelatin  formulation  (10, 15  and  20  percent),  projectile  geometry 
(spheres  of  various  diameters,  cones,  and  parallelepipeds  of  various  constituent 
materials)  and  striking  velocity  were  parametrically  varied.  Target  geometry  was  also 
varied  (homogeneous  parallelepipeds  of  different  lengths  and  widths,  with  and  without 
various  types  of  inclusions,  e.g.,  simulated  femm).  Striking  velocity,  temporally  resolved 
pressure  histories  at  particular  insitu  target  locations,  and  transient  and  permanent 
gelatin  displacements  were  recorded  for  the  various  experiments. 


Col.,  R.F.  Bellamy.  M.D,  FACS,  and  BG.  R.  Zajtchuk,  M.D.,  FACS,  "The 
Physics  and  Biophysics  of  Wound  Ballistics,”  Chapter  4  of  Textbook  of  Military  Medicine: 
Conventional  Warfare,  Office  of  the  Surgeon  General,  Department  of  the  Army,  1991 , 
page  142. 
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Figure  3,  Schematic  of  Ballistic  Test 

The  shadowgraph  (back  lighted  photograph  from  diffuse  light  source)  shown  in  Figure  4 
shows  a  typical  diffraction  pattern  (longitudinal  section)  from  a  low  velocity  impact  of  a 
sphere  on  a  20%  gelatin  target.  Figure  5  shows  a  typical  diffraction  pattern  at  a 
significantly  higher  velocity  recorded  1.33  msec  after  impact.  These  diffraction  patterns 
record  the  transient  displacement  of  the  gelatin  target  after  passage  of  the  projectile. 
These  displacements  eventually  collapse  due  to  the  elastic  recoil  of  the  tissue  resulting  in 
a  much  smaller  permanent  wound  tract  similar  to  that  shown  in  Figure  6. 


Figure  4.  Shadowgraph  of  Low  Velocity  Impact 


Preceding  passage  of  the  projectile  and  coinciding 
with  target  impact  are  large  pressure  transients. 
Figure  7  shows  a  pressure  time  history  from  a 
0.25  inch  steel  sphere  with  a  striking  velocity  of  1,570  ft/sec  on  a  10%  ordnance  gelatin 
target.  This  data  was  recorded  0.5  inches  above  the  projectile  trajectory  2.40  inches 
downstream  of  the  impact  site.  The  measured  peak  pressure  was  1,349  psi. 


■ 
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Figure  5.  Shadowgraph 
of  High  Velocity  Impact 


Recorded  pressure 
versus  time  traces  are 
similar  in  shape  and 
exhibit  an  initial  steep 
rise  to  a  local  maxima 
about  1/2  of  the  peak 
value  that  is  manifested  later  in  the  trace.  This  first  peak  is  followed  by  a  reduced 
pressure  and  gradual  pressure  elevation  to  a  second  peak  which  is  the  maximum 
pressure  recorded.  The  second  peak  is  then  followed  by  an  exponential  decay  to  ambient 
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pressure.  The  pressure 
transients  evolve  over  about 
a  hundred  p,sec  in  these 
experiments  reaching  a 
peak  when  the  projectile  is 
nearest  the  pressure 
transducer.  The  gelatin 
displacement  evolves  over 
hundreds  of  psec  in  the 
wake  of  the  projectile.  Finally,  intense  dynamic  oscillations  affecting  the  entire  target 
are  evident  in  late  time.  These  structural  oscillations  occur  on  the  order  of  milliseconds. 


Peak  pressures  are 
plotted  in  Figure  8  as  a 
function  of  striking 
velocity  for  0.250  inch 
diameter  spherical 
projectiles.  The  data 
identified  by  circles  is 
for  20%  gelatin  targets 
and  the  sqaure  symbols 
represent  10%  gelatin 
targets.  The  data  point 
represented  by  the 
triangle  was  recorded 
closer  than  the  nominal 
0.5  inch  lateral  distance 
from  the  trajectory.  For 
this  data  point  the 
lateral  distance  was 
estimated  to  be  0.25 
inches.  This  point  falls 
well  above  the  other 
data  which  we  have 
attributed  to  the  close 
proximity  of  the 
trajectory. 

Figures  9  and  10  show  a  shadowgraph  and  pressure  trace  from  a  ballistic  test  on  a 
simulated  femur  imbedded  in  a  gelatin  matrix.  The  pressure  experienced  at  the  location 
of  the  pressure  transducer  on  the  surface  of  the  femur  (45  degrees  off-axis)  is  nearly 
2,100  psi  with  a  very  steep  rise  and  sharp  drop.  This  is  typical  of  shock  loading  and 
reflection  and  is  very  different  from  the  pressure  traces  seen  in  homogeneous  targets. 


Figure  7.  Measured  240  microsecond  Pressure  Pulse  from 
Impact 


Figure  6,  Longitudnal  Cross-Section  of  Permenant  Cavity 
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Figure  8.  Peak  Pressures  are  Proportional  to  Striking  Velocity  Squared. 


3000 


Figure  9.  Shadowgraph  of  Ballistic  Experiment  with  Simulated  Femur. 


3.  MODELING 

The  modeling  effort  is  proceeding  in  two  phases.  In  both  phases,  models  describing 
temporally  and  spatially  evolving  stress  and  strain  fields  in  the  neighborhood  of  the 
projectile  trajectory  are  being  developed. 
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Figure  10.  Shock  Front  and  98 
microsecond  Trace  from  Simulated  Femur 
Insitu  Pressure  Transducer 


Models  in  the  first  phase  are  being 
developed  with  reference  to  ballistic 
experiments  on  biosimulants  where 
various  parameters  in  the  experiments 
are  being  systematically  varied.  At  this 
stage  of  the  model  development,  it  is 
assumed  that  the  biosimulants  manifest 
some  phenomenology  in  common  with 
the  corresponding  projectile  interaction 
in  human  tissue.  The  goal  of  the  model 
development  activities  is  to  obtain 
algorithms  with  a  fimctional  form  that  describes  important  phenomenological 
characteristics  of  the  projectile  target  interaction.  In  addition  to  identifying  the 
fimctional  form  of  model  algorithms,  model  parameters  can  be  physically  interpreted. 
This  physical  interpretation  arises  from  seeing  the  effect  on  projectile  terminal  ballistics 
of  parametric  variations  in  ballistic  test  conditions  such  as  target  mechanical  properties, 
projectile  geometry  and  velocity. 

Once  the  fimctional  forms  have  been  posited  and  model  parameters  interpreted,  the 
second  phase  of  the  modeling  activities  begins.  This  second  phase  which  begins  later 
this  year,  concerns  extrapolation  of  the  first  phase  algorithms  to  human  tissue.  This 
extrapolation  is  accomplished  by  substituting  tissue  properties  from  fresh  unembalmed 
cadavers  into  model  parameters  developed  in  the  first  phase  of  the  effort.  The  resulting 
second  phase  models  are  then  compared  with  wound  tract  and  associated  tissue  damage 
documented  in  wound  data  autopsies  (e.g..  Wound  Data  and  Munitions  Effectiveness 
Teams,  WDMET,  from  Vietnam).^ 

In  the  first  year  of  this  project,  algorithms  have  been  developed  which  describe  the 
retardation  of  the  projectile  as  a  function  of  penetration  depth  into  a  multilayer, 
multiphase  viscoelastic  medium.^  Models  have  also  been  developed  which  describe 


^  The  WDMET  data  base,  maintained  at  the  Uniformed  Services  University  of 
Health  Sciences  in  Bethesda,  Maryland  contains  records  of  more  than  8,000  autopsies 
from  battlefield  wounds  in  Vietnam.  These  records  include  over  120,000  slides  along  the 
wound  tract,  original  x-rays,  histiological  sections,  clothing,  projectile  fragments,  and  in 
some  cases  interviews  with  witnesses. 

^  A.K  Chatterjee,  R.D.  Eisler,  and  G.H.  Burghart,  “Penetration  of  Projectiles  in 

Gelatin,”  Werner  Goldsmith  Symposium  sponsored  by  the  Joint  Applied  Mechanics  and 
FOOTNOTES,  Continued 
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target  transient  displacement  based  on  the  quasistatic  generalized  stress-strain  curve  of 
the  target  material.  The  residual  cavity  is  obtained  by  determining  accumulation  of 
inelastic  strain  in  the  target.  Damage  in  affected  areas  is  determined  comparing 
evolved  stress  and  strain  with  empirically  determined  ultimate  values  for  specific  tissues. 

Lateral  displacement  of  target  material  is  determined  by  using  the  velocity  profile  of  the 
projectile  during  penetration  and  projectile  leading  edge  geometry  to  predict  local 
motion  of  target  material  around  toe  projectile.  In  the  event  that  the  projectile  manifests 
unstable  rotational  motion,  the  leading  edge  geometry  is  changed  as  a  function  of  time 
to  mimic  projectile  kinematics.  The  projectile  rotational  kinematics  are  predictable  since 
the  projectile  retardation  model  entails  modeling  the  distribution  of  retarding  forces 
around  the  projectile.''  Key  elements  of  this  model  are  algorithms  which  describe  the 
location  on  the  projectile  where  boundary  layers  shed  and  how  target  material  sticks  to 
the  surface  of  the  projectile  during  penetration.  Currently,  these  algorithms  have 
developed  to  show  parametric  sensitivities  in  terms  of  the  variables  above  and  how 
target  displacement  velocities  vary.  Ballistic  experiments  will  be  used  to  establish  the 
correct  values  of  these  parameters. 


4.  SIMULATION 

The  algorithms  discussed  above  are  used  in  the  formulation  of  a  three  dimensional 
nonlinear  spring  lattice  which  is  used  to  simulate  soft  tissue  response  local  to  the  region 
of  ballistic  insult  (see  Figure  1 1).  When  a  ballistic  injury  occurs,  this  lattice  is 
superimposed  on  the  polygonal  element  representation  of  the  virtual  anatomy. 

The  tissue  stress-strain  curve  is  measured  from  fi'esh  unembalmed  cadavers  and  is  used 
to  calibrate  the  springs  in  the  lattice.  The  velocity  profile  of  the  projectile  along  the 
trajectory  is  used  to  determine  the  locations  and  time  where  the  lattice  is  cut.  The 
algorithms  which  describe  transienf  displacement  of  the  target  material  are  used  to 
determine  initial  velocities  imposed  on  the  spring  network  at  the  location  of  the  lattice 
defects. 


Materials  Summer  Meeting  of  the  American  Society  of  Mechanical  Engineers  and  the 
University  of  California,  Los  Angeles.  28-30  June  1995. 

''  A  detailed  explanation  of  the  methodology  and  comparison  of  experimental 

results  with  model  predictions  is  given  for  a  fin  stabilized  1 9.6  grain  flechette  in:  R.D. 
Eisler,  AK  Chatterjee,  S.F  Stone,  and  G.H.  Burghart,  Algorithm  Development  to 
Describe  Retardation  in  Human  Tissue  of  a  19.6  Grain  Flechette,  Mission  Research 
Corporation  report  MRC-COM-93-372(R1),  U.S.  Army  Natick  Research,  Development, 
and  Engineering  Center  contract  DAAK60-920C-0003,  15  May  1994  (To  be  published  as 
a  Natick  Technical  Report). 
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Figure  11.  Schematic  of  Non-linear  Spring  Lattice 

The  spring  lattice  enables  simulation  of  the  soft  tissue  dynamic  response  from  which  the 
a  description  of  the  residual  wound  tract  is  described  based  on  accumulation  of  inelastic 
strain.  Damage  adjacent  to  the  wound  tract  is  also  described  based  on  the  stress-strain 
field  developed  the  tissue,  e.g.,  the  region  in  which  blood  vessels  rupture. 

Relative  to  bone  trauma;  penetration,  ricochet,  and  fracture  mode  resulting  from  the 
projectile  interaction  is  described  as  well  as  perturbation  of  the  projectile  kinematics. 
Ejection  of  secondary  bone  fragments  and  their  influence  on  soft  tissue  response  is  also 
being  developed.  Convergence  studies  on  how  tissue  mass  can  be  discretized  and  the 
influence  of  terms  in  the  springs  stiffiiess  have  already  been  studied. 
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Figure  12.  Tissue  Displacement  with  Different  Mass  Discretization  Schemes 


Figure  12  shows  maximum  transient  cavity  radius  as  a  function  of  penetration  depth  for 
different  target  mass  discretization  schemes.  Figure  13  shows  scaled  displacement  as  a 
function  of  penetration  depth  for  a  0. 177  cm  diameter,  43 1  mg  sphere,  striking  20% 
gelatin  at  1,780  ft/sec.  The  diameter  of  the  permanent  cavity,  maximmn  extent  of  the 
transient  displacement,  and  affected  region  of  tissue  is  superimposed  on  this  plot.  We 
are  currently  in  the  process  of  correlating  these  results  with  recently  completed  ballistic 
experiments. 
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